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Abstract 

The  research  focused  on  the  development  of  a  new  method  to  identify  damaged  structural 
elements  from  a  large  flexible  space  structure  on-orhit,  using  limited  measured  modal  data.  Limited 
measured  modal  data  is  loosely  defined  as  measured  data  containing  only  a  few  modal  frequencies 
and  less  than  10%  of  the  total  structural  degrees-of-freedom.  This  effort  was  decomposed  into  four 
specific  fasks.  The  flrsf  is  fhe  idenfiflcafion  of  parfial  modal  properfies  from  measured  dafa  of  fhe 
nominal  space  sfrucfure.  Second,  fhe  flnife  elemenf  model  musf  he  adjusfed  fo  mafch  fhe  measured 
nominal  parfial  dafa.  The  fhird  fask  is  an  analysis  of  fhe  exfenf  fo  which  sfrucfural  damage  can  he 
localized  fo  individual  sfrucfural  elemenfs  using  fhe  measured  dafa.  In  conjuncfion  wifh  fhis  fask  is 
fhe  deferminafion  of  where  fo  hesf  place  fhe  limifed  number  of  sensors  on  fhe  sfrucfure.  Lasfly,  fhe 
idenfiflcafion  of  sfrucfural  damage  musf  be  performed  using  fhe  limited  measured  modal  dafa  from  a 
damaged  space  sfrucfure. 

Idenfiflcafion  of  fhe  modal  parameters  was  accomplished  using  fhe  Eigensysfem  Realizafion 
Algorifhm,  a  lime  domain  based  melhod,  adopled  for  use  wifh  averaged  measured  frequency  response 
funclions.  Model  funing  was  performed  using  fhe  Automated  Sfrucfural  Oplimizafion  Soflware  pack¬ 
age,  adapted  for  model  funing.  The  melhod  minimizes  a  cosl  funclion  based  on  fhe  mismalch  belween 
fhe  measured  and  analylical  eigenslruclure.  The  minimizalion  is  solved  using  fhe  eigenvalue  and  eigen¬ 
vector  sensilivilies  al  each  ileralion  step.  The  deferminafion  of  priorilized  sensor  localions  and  damage 
localizalion  is  performed  using  fhe  eigenvalue  and  eigenvector  sensilivilies.  Damage  idenfiflcafion  is 
performed  using  a  newly  developed  assigned  parfial  eigenslruclure  melhod,  which  determines  required 
sliffness  changes,  consislenl  wifh  fhe  finite  elemenf  formulalion  to  achieve  fhe  measured  dafa. 

The  Iheory  for  each  fask  is  presented  and  illuslraled  on  an  analylical  example  of  a  41 -elemenf 
free-free  planar  Iruss.  Two  experimenlal  demonslralions  were  performed  and  fhe  resulls  reported.  The 
flrsf  was  a  canlilevered  Iruss  modeled  wifh  104  rod  elemenfs  wifh  a  lolal  of  96  degrees-of-freedom.  The 
measured  dafa  consisted  of  fhe  flrsf  five  flexible  modal  frequencies,  and  only  eighl  componenls  of  fhe 
five  corresponding  eigenvectors.  The  second  experimenl  was  a  canlilevered  frame  assembly  modeled 
wifh  96  beam  elemenfs  wifh  a  lolal  of  192  degrees-of-freedom.  The  measured  dafa  consisted  of  fhe  flrsf 
eighl  flexible  modal  frequencies,  and  only  eighl  componenls  of  fhe  eighl  corresponding  eigenvectors. 


XVI 


The  research  showed  that  in  each  test  case,  the  structural  damage  could  successfully  he  localized  to  a 
small  portion  of  the  structure.  The  extent  to  which  damage  can  he  localized  was  limited  hy  hoth  model 
fidelity  and  accuracy  of  the  measured  modes. 


STRUCTURAL  DAMAGE  IDENTIFICATION  FROM  EIMITED  MEASUREMENT 


DATA 


/.  Introduction 

Over  the  past  decade,  both  the  size  and  complexity  of  military  and  civilian  spacecraft  have 
increased  considerably.  Future  Air  Force  space  systems  will  continue  to  increase  in  size  due  to  mission 
requirements,  while  minimizing  weight  to  remain  within  launch  constraints  of  the  booster  inventory. 
In  addition  to  the  size  and  complexity  increase,  the  on-orbit  lifetime  has  increased  as  a  result  of 
new  advances  in  solar  power  generation  and  battery  storage  devices.  These  advances  coupled  with 
NASA’s  demonstrated  ability  to  perform  on-orbit  repair  and  replacement  of  flight  critical  items,  further 
extend  a  satellite’s  usable  lifetime.  As  a  result,  future  large  flexible  space  structures  will  have  an 
unprecedented  requirement  for  verifying  the  structural  integrity  of  such  space  structures  on  orbit,  on  a 
periodic  basis  over  the  lifetime  of  the  space  system.  The  potential  degradation  of  structural  components 
from  long  term  exposure  to  the  space  environment  was  documented  with  the  retrieval  of  NASA’s  Long 
Duration  Exposure  Facility. ''''''  Additional  damage  is  also  possible  due  to  loads  introduced  during 
orbital  maneuvers,  spacecraft  docking  operations,  and  from  collisions  with  space  debris.  NORAD 
currently  tracks  over  5000  objects  10  centimeters  in  size  or  larger  in  low  Earth  orbit. Information 
on  both  the  location  and  extent  of  structural  damage  will  be  critical  in  assessing  required  in-space 
repair  missions  and/or  deviations  from  the  planned  mission  profiles.  The  current  configuration  of  the 
International  Space  Station  Alpha  is  shown  in  Eigure  1.1.  This  research  focuses  on  identification  of 
damage  to  structural  sub-assemblies  typical  of  large  orbiting  space  platforms.  The  truss-like  structure 
of  International  Space  Station  Alpha’s  solar  array  sub-assembly  is  depicted  in  Eigure  1.2.'^*' 

One  method  of  structural  verification  is  visual  inspection,  however  this  method  may  be  impracti¬ 
cal  due  to  the  extra- vehicular  activity  man-hours  required,  or  impossible  as  in  the  case  of  unmanned  and 
high  altitude  missions.  Thus  an  alternative  solution  to  visual  inspection  is  desired,  and  is  the  main  focus 
of  this  research.  Two  current  disciplines,  closely  related  to  this  topic,  are  system  identification  theory 
and  analytical  model  tuning  from  experimental  data.  These  two  disciplines  will  be  briefly  reviewed, 
and  their  correlation  to  structural  damage  identification  will  be  highlighted. 
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Figure  1 .2  Solar  array  sub-assembly  from  International  Space  Station  Alpha. 
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Vibration  testing  has  been  used  extensively  in  industry  to  monitor  rotating  machinery.  From 
observations  of  input/output  relationships,  the  frequency  signature  of  a  system  can  be  used  to  detect 
failures.  Identifying  specific  failure  modes  can  then  be  accomplished  by  a  comparison  of  the  frequency 
signature  to  that  of  a  predetermined  set  for  each  failure  mode.  For  a  large  system,  this  failure  mode 
set  may  be  impractical  to  produce  or  store  and  interpolation  is  not  necessarily  possible  between  failure 
modes.  Thus  an  alternative  to  searching  a  frequency  signature  database  is  sought.  One  alternative  is 
the  use  of  a  system  identihcation  algorithm. 

System  identihcation  techniques  are  based  on  determining  the  underlying  physical  system  from 
a  given  set  of  input/output  relationships.  For  non-parametric  system  identihcation,  the  physical  system 
can  be  viewed  as  a  ‘black  box’  and  thus  no  inherent  knowledge  of  the  structure  of  the  system  is 
required.  The  primary  measure  of  the  effectiveness  of  the  system  identihcation  is  in  how  well  the 
identihed  mathematical  model  produces  an  output  which  matches  the  measured  output  for  a  given 
input  signal.  Since  there  is  a  direct  relationship  between  the  time  and  frequency  domain  through 
the  Fourier  transform,  the  identihcation  can  be  accomplished  in  either  domain.  System  identihcation 
methods  have  been  shown  to  be  very  effective  in  producing  models  which  exactly  (theoretically)  or 
closely  match  (to  within  the  experimental  error)  the  true  system;  however,  they  typically  do  not  directly 
give  information  about  the  physical  structure  of  the  system.  For  structural  damage  identihcation,  the 
system  identihcation  techniques  are  useful  in  obtaining  the  eigenstructure  of  the  physical  system,  which 
will  be  shown  to  be  sufficient  to  determine  structural  damage.  The  feasibility  of  identifying  modal 
frequencies,  damping  ratios  and  shapes  from  on-orbit  testing  has  been  addressed  in  previous  studies. 

Unlike  the  non-parametric  identihcation,  model  tuning  is  a  parametric  approach  to  system 
identihcation  which  includes  a  detailed  physical  model  with  well  dehned  parameters  of  uncertain 
values.  Model  tuning  attempts  to  match  the  input/output  relationships  of  an  analytical  model  and 
the  physical  system  by  varying  parameters  in  the  mathematical  model.  Differences  in  model  tuning 
methods  depend  on  the  assumptions  made  on  how  to  vary  these  parameters.  For  physical  structures, 
one  approach  is  to  vary  the  material  properties  of  the  elements  in  the  model.  This  is  typically  done 
by  adjusting  the  stiffness  and  mass  and/or  damping  matrices  of  the  model.  This  method  assumes  the 
existence  of  a  hnite  element  method  (FEM)  model  which  produces  simulated  output  which  is  reasonably 
close  to  experimental  measurements.  To  exactly  match  experimental  data,  only  small  perturbations  to 
the  model  matrices  are  required.  Model  tuning  can  then  be  approached  as  an  optimization  problem. 
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For  a  given  input,  a  minimization  on  the  difference  between  the  model  output  and  the  measured  output 
can  be  performed  through  the  use  of  an  appropriate  metric,  while  iterating  on  the  perturbation  of  the 
matrices.  The  key  distinction  between  different  model  tuning  methods,  important  to  structural  damage 
identification,  is  whether  or  not  the  physical  connectivity  of  the  structure  is  retained.  Clearly,  to 
correctly  determine  a  change  in  a  structural  element,  the  tuned  model  must  not  remove  load  paths  or 
introduce  load  paths  not  present  in  the  physical  structure. 

An  additional  complexity  with  model  tuning  algorithms  is  in  obtaining  complete  experimental 
information.  A  typical  FEM  model  of  a  large  flexible  space  structure  will  include  hundreds  or  thousands 
of  nodes  with  as  many  as  six  degrees-of-freedom  at  each  node.  In  contrast,  typical  experimental  data 
will  include  accurate  information  on  only  the  lower  frequency  modes  of  vibrations,  taken  at  only  a  small 
subset  of  the  nodes  in  the  FEM  model.  Eurthermore,  typically  only  translational  degrees-of-freedom  are 
measured  which  further  reduces  the  available  data.  Thus  an  additional  distinguishing  feature  between 
model  tuning  algorithms  is  in  how  to  incorporate  the  reduced  experimental  data  sets  into  the  model. 
Several  model  reduction/expansion  algorithms  have  been  used  to  correlate  the  model  to  experimental 
data.  A  common  attribute  of  these  methods  is  the  use  of  the  nominal  EEM  model  to  obtain  the 

unmeasured  degrees-of-freedom.  Eor  damage  detection  however,  the  EEM  model  is  the  unknown  and 
hence  cannot  be  used  as  the  basis  for  the  reduction/expansion. 

Having  briefly  motivated  the  problem  of  on-orbit  damage  identification,  along  with  introducing 
the  concepts  of  system  identification  and  model  tuning,  a  statement  of  the  research  objective  can  now 
be  given.  The  research  focused  on  the  development  of  a  new  method  to  identify  damaged  structural 
elements  from  a  large  flexible  space  structure  on-orbit,  using  limited  measured  modal  data.  This  effort 
was  decomposed  into  four  specific  fasks.  The  flrsf  is  fhe  idenfiflcafion  of  parfial  modal  properfies  from 
measured  dafa  of  fhe  nominal  (i.e.  undamaged)  space  sfrucfure.  Second,  fhe  EEM  model  musf  be 
adjusted  fo  mafch  fhe  measured  nominal  parfial  dafa.  The  fhird  fask  is  an  analysis  of  fhe  exfenl  fo 
which  sfrucfural  damage  can  be  localized  fo  individual  sfrucfural  elemenfs  using  fhe  measured  dafa.  In 
conjuncfion  wifh  Ibis  fask  is  fhe  deferminafion  of  where  fo  besf  place  fhe  limited  number  of  available 
sensors  on  fhe  sfrucfure.  Easfly,  fhe  idenfiflcafion  of  sfrucfural  damage  musf  be  performed  using  fhe 
measured  dafa  from  a  damaged  space  sfrucfure. 

This  chapter  briefly  infroduced  fhe  research  efforf  and  ouflined  some  of  fhe  related  work  in  Ibis 
area.  Chapter  II  presenfs  a  background  of  related  work  in  system  idenfiflcafion  and  model  Inning  and 
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discusses  the  relevance  to  damage  identification.  The  subsequent  chapters  develop  the  four  tasks  of  the 
on-orhit  damage  identification  problem  in  detail,  beginning  with  an  overview  of  the  methodology  for 
each  task  in  Chapter  III.  The  four  tasks  were  integrated  into  a  single  software  package  which  is  presented 
Appendix  B.  Lastly,  validation  of  the  research  effort  was  performed  using  laboratory  experiments  which 
exhibit  the  same  dynamic  properties  as  large  flexible  space  structures,  and  is  presented  in  Chapter  VIII. 
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//.  Background 


2. 1  Motivation 

The  problem  of  producing  analytical  models  capable  of  predicting  dynamic  responses  has  been 
widely  studied.  While  these  solution  methods  may  result  in  analytical  models  which  match  (typically 
in  a  weighted  least  squares  sense)  experimental  data,  they  most  often  do  not  directly  address  the 
problem  of  relating  perturbations  of  the  analytical  model  to  changes  in  the  physical  parameters  of  the 
structure.  Direct  identification  of  failed  structural  components  is  exactly  the  information  required  to 
adjust  mission  profiles  fo  minimize  sfrucfural  dynamic  loading  and  fo  enable  developmenf  of  repair 
missions  where  possible.  Thus  a  mefhod  which  direcfly  idenlifies  damaged  sfrucfural  members  from 
experimenfal  dafa  is  highly  desirable.  Due  fo  fhe  similarifies  among  damage  idenlificalion,  system 
idenlificalion  and  model  funing,  an  overview  of  exisfing  mefhods  for  system  idenlificalion  and  model 
Inning  will  be  presented  as  well  as  currenl  mefhods  of  damage  idenlificalion. 

2.2  System  Identification  Techniques 

The  problem  associaled  wilh  system  idenlificalion  is:  given  fhe  measured  response  fo  a  known 
inpul,  determine  a  malhemalical  represenlalion  of  fhe  system  which  reproduces  fhe  oulpul  sequence 
given  fhe  inpul  sequence.  Differences  in  algorilhms  are  based  on  fhe  assumplions  of  fhe  underlying 
system,  which  Ihen  eslablishes  fhe  slruclure  of  fhe  analylical  model.  For  large  flexible  space  slruclures, 
a  finite-dimensional  linear  lime-invarianl  model  can  adequately  represenl  fhe  dynamics  of  fhe  slruclure 
and  hence  will  be  Ihe  subjecl  of  Ihis  invesligalion.  In  general  system  idenlificalion  Iheory,  such  a 
reslriclive  assumplion  need  nol  be  imposed. 

To  illuslrale  Ihe  general  approach  lo  system  idenlificalion,  a  sampled  dala  single-inpul  single- 
oulpul  system  will  be  considered.  This  presenlalion  is  based  on  work  by  Ljung.  Given  an  inpul 
sequence  u{t)  and  an  oulpul  sequence  y{t)  where  t  =  (0, 1, . . . ,  N  —  1),  Ihe  discrele-lime  system  can 
be  represented  as: 

y{t)  =  G{z)u{t)  +  v(t)  (2.1) 
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with  G{z)u{t)  representing  the  convolution  summation  of  the  system’s  impulse  response  sequence 
g{t)  with  the  input  sequence  u{t),  and  is  given  hy: 

OO 

G{z)u{t)  =  ^  g{k)u{t  -  k)  (2.2) 

k  =  0 

The  function  v{t)  represents  the  measurement  noise  (uncertainty)  associated  with  obtaining  the  output 
sequence.  The  term  G{z)  represents  the  system’s  transfer  function  and  is  defined  as: 

OO 

G{z)  =  Y,9itK'  (2-3) 

t  =  0 

where  ^  represents  the  time-shift  operator  such  that: 

z~^y{k)  =  y{i  -  1)  (2-4) 

A  parametric  approach  to  system  identification  is  then  confined  fo  fhe  deferminafion  of  a  sef  of 
paramefers  which  fully  describe  G{z).  To  illnsfrafe  fhe  differenl  approaches,  a  more  general  form  of 
Equafion  (2.1)  is  used  fo  form  a  general  paramefric  model  given  as  fhe  following  convolufion: 

,  ,  ,  ,  B{z)  .  ,  C{z)  .  , 

Mz)y{t)  =  ^^“(0  +  (2.5) 

where  A{z),  B(z),  C(z),  D(z),  and  F(z)  are  polynomials  in  fhe  lime  shifl  operator  2;  as  follows: 

A(z')  =  1  +  aiZ  ^  +  .  .  .  +  Clna^  (2.6) 

and  similarly  for  fhe  remaining  four  polynomials.  The  order  of  fhe  polynomials  is  given  by  na,nb,nc,  nd 
and  nf  respeclively.  The  sequence  e{t)  is  assumed  to  be  a  while  noise  sequence  and  is  shaped  by 
A{z)D{z)  produce  fhe  measuremenl  noise  v{t)  of  Equafion  (2.1).  The  commonly  used  linear  models 
are  all  special  cases  of  fhe  general  form  of  Equafion  (2.5).  The  auloregressor  wilh  exogenous  inpul 
model  is  oblained  by  selling  nc  =  nd  =  nf  =  0.  Similarly,  fhe  auloregressor  wilh  moving  average 
and  exogenous  inpul  is  oblained  by  selling  nf  =  nd  =  0.  The  oulpul  error  model  corresponds  to 
na  =  nc  =  nd  =  0,  while  Ihe  Box-Jenkins  model  is  oblained  by  selling  na  =  0.  A  large  sel  of 
models  can  Ihus  be  oblained  from  Ihe  general  form  given  in  Equafion  (2.5). 
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Although  the  above  development  was  in  transfer  function  form,  an  equivalent  state-space  repre¬ 
sentation  is  also  possible.  The  state-space  representation  of  Equation  (2.5)  is  given  as: 

x{t  +l)  =  Ax(t)  +  Bu{t)  (2.7) 

y{t)  =  Cx(t)  +  Du(t)  +  v(t)  (2.8) 

where  A,  B,  C,  and  D  are  redefined  as  constant  matrices.  Equation  (2.7)  can  be  rewritten  in  terms  of 
the  time  shift  operator  as: 

zx(t)  =  Ax{t)  -\-  Bu(t)  (2.9) 

Solving  for  x(t)  yields: 

x(t)  =  (zl  —  A)~^  Bu(t)  (2.10) 

Substituting  this  result  into  Equation  (2.8)  yields  the  convolution: 

yit)  =  [C  (zl  -  A)-'  B  +  d]  u{t)  +  v{t)  (2.11) 

Thus  the  transfer  function  is  directly  related  to  the  state-space  quadruple  (  A,  B,C,  D)by. 

G(z)  =  C[zl  -  A]-^B  +  D  (2.12) 

Hence  the  system  identification  problem  is  equivalent  to  determining  the  constant  matrices  of  the 
quadruple  (A,B,C,D)  which  accurately  reproduces  a  given  measured  response  from  a  known  input. 
A  further  distinction  between  techniques  involves  whether  the  input  u(t)  and  output  y  (f )  are  represented 
in  the  time  or  frequency  domain.  The  foundation  of  these  two  methods  will  be  briefly  discussed. 

2.2.1  Time  Domain.  System  identification  in  the  time  domain  seeks  to  reconstruct  the 
transfer  function  based  on  identification  of  the  system’s  impulse  response  function.  Eor  the  scalar  case, 
if  the  input  is  a  pulse  defined  as: 

(a  t  =  0 

u(t)  =  {  (2.13) 

[  0  f  /  0 

fhen  fhe  oufpuf  y(t)  becomes: 

y(t)  =  ag(t)  +  v(t)  (2.14) 
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where  g(t)  is  the  unit-impulse  response  sequence  and  v(t)  is  the  measurement  noise.  If  the  signal  to 
noise  ratio  is  high,  then  v(t)  <<  ag(t)  and  can  he  neglected.  An  estimate  of  the  impulse  response  is 
then  simply: 


5(0  = 


5(0 


(2.15) 


with  an  error  of  .  An  estimate  of  the  transfer  function  can  then  he  determined  directly  from  the 
estimated  impulse  response. 

OO 

Giz)  =  J2m^~'  (2-16) 

t=o 

A  disadvantage  with  the  above  method  is  the  requirement  to  excite  the  system  with  an  impulse. 
This  can  he  approximated  with  a  finite  pulse  input,  hut  may  he  impractical  to  implement  on  certain 
systems,  and  it  may  he  difficult  to  obtain  response  levels  which  are  significantly  above  the  measurement 
noise  levels.  An  alternative  procedure  known  as  correlation  analysis  overcomes  this  shortcoming  by 
requiring  the  input  is  a  zero-mean  white  noise  sequence,  i.e.  it’s  auto-correlation  function  is  given  as: 


A  r  =  0 


(0  r  /  0 

Then,  from  Equation  (2.1)  the  cross-correlation  between  the  input  and  output  is  given  as: 


(2.17) 


Rynir)  =  XgiT) 


(2.18) 


where  g(T)  is  the  desired  impulse  response  sequence.  By  definition, the  cross-correlation  between 
two  zero-mean  sequences  is  given  as: 


1 


N-l 


Rynir)  =  ^  2^  vitHt  -  r) 


(2.19) 


where  it  is  assumed  that  the  limit  exists.  An  estimate  of  the  correlation  is  then  obtained  by  selecting  N 
sufficiently  large.  Thus  the  correlation  estimate  is  defined  as: 


1 


N-l 


^ynir)  =  ^^5(0“(^-^) 


(2.20) 


t  =  0 
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Using  the  estimated  correlation,  an  estimate  of  gir)  can  then  he  obtained  from  the  input  and  output 
sequences  using  Equations  (2.18)  and  (2.20). 

-y  N-1 

=  Xn  ^  -  t)  (2-21) 

Although  the  assumption  was  made  that  the  input  u{t)  was  white  noise,  this  is  not  a  restrictive 
assumption  since  a  shaping  filter  can  always  he  constructed  to  produce  the  actual  sequence 
The  output  sequence  can  then  he  filtered  through  this  same  filter  and  the  estimate  is  computed  using 
Equation  (2.21)  with  filtered  data.  Although  presented  here  for  the  scalar  case,  the  method  is  applicable 
to  multiple  input/output  combinations. 

Based  on  these  concepts,  a  vast  array  of  system  identification  methods  have  been  studied  and 
presented  in  the  literature.  Erom  these  methods,  the  eigensystem  realization  algorithm  (ERA)  was 
singled  out  as  one  which  has  demonstrated  it’s  ability  to  accurately  produce  a  minimum  system 
realization  from  multi-input/multi-output  (MIMO)  experimental  data  typical  of  large  flexible  space 
structures.  The  ERA  method  is  presented  in  Chapter  IV. 

2.2.2  Frequency  Domain.  The  aforementioned  methods  were  based  on  reconstruction  of  the 
impulse  response  in  the  time  domain.  An  alternative  method  is  identification  of  the  system’s  transfer 
function  in  the  frequency  domain.  This  method  is  commonly  referred  to  as  spectral  analysis  and  is 
included  here  for  completeness. 

The  frequency  response  function  of  the  system  in  Equation  (2.1)  is  found  by  evaluating  the 
transfer  function  G{z)  on  the  unit  circle  2;  =  ,  Eurthermore,  the  noise  sequence  v{t)  can  be 

described  in  terms  of  its  power  spectral  density  <i>t,  1,(02),  where  the  power  spectral  density  function  is 
defined  as  fhe  Eourier  fransform  of  fhe  aufocorrelafion  funcfion.  Thus  fhe  inpuf  oufpuf  relafionship  of 
Equafion  (2.1)  is  completely  specified  from  knowledge  of  G{ijj)  and  <i>t,t,(o2).  If  inifially  we  assume 
fhe  measuremenf  noise  is  negligible,  fhe  frequency  response  funcfion  can  be  found  from  sine-dwell 
fesfing.  In  fhis  mefhod,  fhe  system  is  given  a  known  sinusoidal  inpuf. 

u{t)  =  a  sin(o2f )  (2.22) 
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For  a  linear  time-invariant  system,  it  is  well  known  that  the  steady-state  output  response  after  the 
transients  die  out  is  given  hy: 


Vssit)  =  a  |G'(c^)|  siniut  +  (j))  +  v{t)  (2.23) 

where  cj)  is  given  hy  the  phase  angle  of  G{ijj).  In  this  way,  an  estimate  of  the  transfer  function  can 
he  found  hy  repeated  application  of  the  input  sinusoid  at  frequencies  of  interest  and  measuring  the 
magnitude  and  phase  of  the  response.  A  reconstruction  of  the  transfer  function  from  the  data  can  then 
he  performed. 

In  the  preceding  analysis  it  was  assumed  that  the  measurement  noise  was  small  and  hence  was 
neglected  from  the  analysis.  A  less  restrictive  assumption  is  to  assume  u{t)  and  v{t)  are  independent, 
hut  v{t)  not  necessarily  negligible.  The  relationship  between  input  and  output  given  in  Equation  (2.1) 
corresponds  to  a  relationship  between  the  power  spectra  of: 

(2.24) 


and 


(2.25) 


Estimates  of  the  frequency  response  function  and  the  noise  power  spectrum  are  obtained  by  computing 
the  appropriate  estimates  of  the  correlation  functions  using  Equation  (2.20)  and  their  corresponding 
power  spectra  (by  Eourier  transformation). 


G{u) 


(2.26a) 


(jJ 


(2.26b) 


The  above  highlights  the  general  aspects  of  spectral  analysis.  A  more  complete  derivation  of  spectral 
analysis  is  contained  in  Newland. The  preceding  material  was  based  on  a  discrete  time  analysis, 
since  this  is  the  form  of  measured  data.  The  material  presented  in  the  remainder  of  this  chapter  and  all 
subsequent  chapters  is  based  on  a  continuous  time  analysis  unless  otherwise  stated. 
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2.3  Model  Tuning  Techniques 


As  opposed  to  system  identification  techniques  where  no  a  priori  model  is  required  to  perform 
the  identification,  model  tuning  assumes  the  existence  of  an  analytical  model.  For  flexible  structures, 
the  model  can  he  written  as: 

Mx  +  Cx  +  Kx  =  Ff  (2.27) 


where  M,  C,  K  represent  the  mass,  damping,  and  stiffness  matrices  respectively  for  the  degree-of- 
freedom  vector  x.  The  matrix  F  represents  the  influence  matrix  to  a  given  input  vector  /.  This 
analytical  model  is  the  result  of  an  application  of  Newton’s  second  law,  and  most  likely  produced 
as  the  output  of  a  FEM  program  and  used  to  simulate  the  dynamic  response  of  the  actual  structure. 
Note  that  the  matrices  M,  C,  If,  F  of  the  second-order  Equation  (2.27)  can  he  related  to  the  first-order 
quadruple  (  A,  B,  C,  D)  as  given  in  Equations  (2.7)  and  (2.8).  The  transformation  from  second-order 
to  first-order  form  for  a  state  vector  [x,  xY  is  given  hy: 


A  = 


0 

-M-^K 


I 

-M-^C 


B  = 

C  = 


0 

M-^F 

[I  0] 


D  =  m 


(2.28a) 

(2.28h) 

(2.28c) 

(2.28d) 


Due  to  violation  of  modeling  assumptions  and  to  inherent  uncertainties  in  material  properties,  joint 
properties,  boundary  conditions,  etc.  in  developing  the  analytical  model,  there  will  exist  an  imperfect 
correlation  between  the  experimental  and  analytical  results.  If  the  correlation  results  are  unsatisfactory, 
an  adjustment  to  the  finite  element  model  is  necessary.  The  common  attribute  of  model  tuning 
techniques  is  that  they  attempt  to  minimize  the  required  modification  to  the  matrices,  assuming  the 
EEM  model  is  a  reasonable  approximation  to  the  physical  structure.  Additionally,  for  realistic  structures, 
the  EEM  model  will  include  a  large  number  of  degrees-of-freedom  (dimension  of  a;)  of  which  only  a 
small  subset  will  be  measured.  The  differences  in  the  techniques  stem  from  how  to  incorporate  the 
reduced  data  set  and  on  what  restrictions  to  make  on  the  set  of  allowable  perturbations  to  the  existing 
matrices  as  well  as  the  selection  of  what  error  metric  is  to  be  minimized.  Eurther  classifications  can 
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be  made  on  whether  damping  is  considered  in  both  the  model  and  the  measurements.  For  cases  where 
damping  is  neglected,  C  =  0  in  Equation  (2.27),  the  distinction  is  whether  adjustments  are  made  to 
the  mass  matrix,  stiffness  matrix,  or  both.  The  adjustments  are  the  result  of  the  solution  of  eigenvalue 
problems,  orthogonality  relationships,  and  eigenstructure  assignment  techniques.  Several  relevant 
techniques  will  now  be  reviewed. 

23.1  Sensitivity  Based.  A  straightforward  approach  to  adjusting  the  matrices  in  Equation 
(2.27)  is  to  establish  an  objective  function  which  is  a  measure  of  the  difference  between  the  experimental 
and  analytical  model  data. The  eigenstructure  (eigenvalues  and  eigenvectors)  has  been  widely 
accepted  as  an  acceptable  method  of  measuring  correspondence  between  models. This  is  due  to  the 
fact  that  the  eigenstructure  is  invariant  under  a  state  transformation  and  thus  allows  direct  comparisons 
of  different  realizations.  Sensitivity  is  the  quantihcation  of  changes  in  the  objective  function  due  to 
changes  in  the  design  variables.  The  sensitivity  of  the  objective  function  with  respect  to  the  design 
parameters  is  determined  for  each  parameter  in  the  matrices.  A  numerical  optimization  problem  is 
then  solved  by  searching  over  the  entire  parameter  space  using  the  sensitivity  information  to  determine 
search  directions.  An  advantage  of  this  method  is  that  the  updated  models  are  consistent  with  the  EEM 
formulation  and  thus  the  connectivity  is  preserved. 

While  this  method  can  yield  acceptable  results  for  small  problems,  the  major  drawback  is  the 
large  number  of  design  parameters  contained  in  the  search  space.  Eor  larger  problems,  the  selection 
of  a  suitable  set  of  parameters  to  search  over  may  not  be  intuitive.  Its  advantage  in  structural  damage 
identihcation,  however,  is  that  a  set  of  parameters  can  be  chosen  that  reflect  physical  failure  modes, 
such  as  searching  over  an  elastic  modulus  value  or  the  cross-sectional  area  of  each  member.  Preliminary 
knowledge  of  the  damage  location  may  be  required  for  large  structures  to  minimize  the  required  search 
space.  Hemez  and  Earhat  applied  this  procedure  to  damage  detection.  Success  of  the  detection 
algorithm  relied  on  the  extent  of  the  damage  and  the  sensitivity  to  the  chosen  design  variables. 

2.3.2  Orthogonality  Based.  If  the  underlying  physical  system  is  assumed  to  be  linear  and 
either  undamped  or  proportionally  damped,  then  the  mode  shapes  of  the  structure  will  be  orthogonal 
with  respect  to  the  mass  matrix.  Modal  orthogonality  can  then  be  represented  as: 


=  I 


(2.29) 
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where  $  is  the  matrix  composed  of  the  measured  eigenvectors  (mode  shapes).  Using  this  orthogonality 
requirement,  small  adjustments  can  then  he  made  to  the  mass  matrix  M  to  force  the  orthogonality 
requirement.  Note  that  confidence  in  the  analytical  mass  matrix  is  assumed  higher  than  the  confidence 
in  fhe  sfiffness  mafrix.  For  fhe  undamped  sfrucfure,  fhe  corresponding  eigenvalue  problem  can  fhen  he 
wriffen  as: 

=  A  (2.30) 


Adjusfmenls  fo  fhe  sfiffness  mafrix  are  fhen  made  such  fhaf  Equafion  (2.30)  is  safisfied.  In  Equafions 
(2.29)  and  (2.30),  hofh  <i>  and  A  are  required  fo  have  fhe  same  dimension  as  fhe  analyfical  model. 
As  previously  discussed,  fhe  sefs  of  measured  eigendafa  will  he  much  smaller  fhan  fhe  number  of 
degrees-of-freedom  in  fhe  analyfical  model.  Thus,  fo  apply  fhe  orfhogonalify  equafions,  eifher  a  modal 
expansion  on  fhe  measured  dafa  or  a  modal  reducfion  on  fhe  analyfical  model  musf  be  performed.  As 
presented  by  Berman,  a  modal  expansion  can  be  performed  on  fhe  fesf  dafa  as  follows. The  analyfical 
mafrices  are  firsl  ordered  such  fhaf  fhe  coordinates  associated  wifh  fhe  measured  degrees-of-freedom 
are  ordered  above  fhe  remaining  coordinates.  The  eigenvalue  problem  for  each  mode  can  fhen  be 
wriffen  as: 

[K  -  =  0  (2.31) 


This  equafion  can  fhen  be  parfifioned  as: 


ki  k2 

kl  k4 

The  unknown  lower  parfifion  of  fhe  eigenvector  can  fhen  be  solved  using  fhe  following: 


(2.32) 


(f>2i  =  -(ki  -  Xinii)  ^k^  - 


(2.33) 


This  resulf  is  repeated  for  each  measured  mode.  This  mefhod  represenfs  an  inferpolafion  of  fhe  measured 
modes  af  fhe  unmeasured  degrees-of-freedom.  The  accuracy  of  fhis  technique  is  clearly  dependenf  on 
how  accurately  fhe  analyfical  model  represented  fhe  physical  sfrucfure.  If  should  be  noted  fhaf  if  fhe 
frequency  dependence  in  fhe  above  equafion  were  ignored,  fhis  would  resulf  in  fhe  sfandard  Guyan 
reducfion  relafionship. 
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Having  obtained  full-length  eigenvectors,  the  adjusted  mass  matrix  is  then  formed  from  the 
solution  of  a  constrained  minimization  problem. An  objective  function  of  the  form: 

J  =  (2.34) 

is  used  which  penalizes  deviations  from  the  analytical  model.  In  the  original  work  the  Frobenius 
norm  was  chosen.  The  subscript  a  denotes  the  analytical  model.  The  orthogonality  constraint  is  then 
appended  with  Lagrange  multipliers  as  follows: 

n  n 

L  =  J  +  -  /),,  (2.35) 

i=i  j=i 

Differentiating  L  with  respect  to  M  and  setting  derivatives  equal  to  zero  yields  a  solution  to  the  adjusted 
mass  matrix  which  now  satishes  the  orthogonality  relationships  and  has  minimized  deviations  from 
the  analytical  mass  matrix  according  to  the  applied  norm  in  the  objective  function.  The  corrected  mass 
matrix  is  given  by: 

M  =  Ma^m~^(I  —  Ma  (2.36) 

with  m  dehned  as: 

m  =  (2.37) 

Although  not  presented  in  the  original  work,  an  iterative  scheme  could  be  used  becuase  the  interpolated 
value  of  <i>  will  change  for  the  adjusted  M.  Having  obtained  the  adjusted  mass  matrix,  the  stiffness 
matrix  is  then  obtained  in  a  similar  manner.  The  objective  function  is  formed  as  an  appropriate  weighted 
norm  on  the  difference  between  the  analytical  and  adjusted  stiffness  matrices.  In  Berman’s  work  which 
was  an  extension  of  the  stiffness  matrix  adjustment  method  of  Baruch,'^'  the  analytical  mass  matrix 
was  chosen  as  the  appropriate  weighting  function  for  the  individual  elements.  This  seems  appropriate 
only  when  elemental  stiffness  values  are  directly  related  to  the  mass  of  each  element.  Certainly, 
there  are  structures  which  contain  mass  elements  which  do  not  contribute  to  the  overall  stiffness  of  the 
structure.  With  the  norm  chosen,  the  eigenvalue  equation  is  then  appended  as  a  constraint  with  Lagrange 
multipliers.  An  additional  constraint  to  ensure  symmetry  of  the  stiffness  matrix  is  also  appended.  The 
adjusted  stiffness  matrix  is  then  found  from  a  differentiation  with  respect  to  the  elements  of  Ka  and 
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setting  them  equal  to  zero.  The  resulting  adjusted  stiffness  matrix  is  given  as: 


K  =  Ka  +  A  +  (2.38) 

where 

A  =  i  (2.39) 

Although  the  presentation  above  based  on  the  work  of  Berman  and  Baruch  did  not  consider  damping, 
it  could  be  expanded  to  include  the  proportional  damping  case.  The  advantage  of  this  method  is  the 
straightforward  mathematical  formulation  of  the  adjusted  matrices.  Furthermore,  this  method  will 
result  in  a  tuned  analytic  model  which  exactly  reproduces  the  experimental  data.  Its  shortcoming  is 
that  it  does  not  guarantee  the  closeness  to  unmeasured  modes  not  used  in  the  tuning  process.  This  is 
a  result  of  potentially  unrealistic  changes  in  the  stiffness  matrix  such  as  the  introduction  of  load  paths 
which  physically  do  not  exist.  The  introduction  of  fictitious  load  paths  is  undesirable  when  identifying 
damaged  elements  since  the  stiffness  changes  can  not  be  directly  related  to  physical  elements. 

The  orthogonality  based  approach  described  above  is  essentially  a  matrix  update  approach 
in  which  perturbation  mass,  damping,  or  stiffness  matrices  are  determined  that,  when  added  to  the 
analytical  matrices,  reproduces  the  measured  result.  Brock  originally  examined  this  problem  while 
enforcing  symmetry  and  positive-definiteness  of  the  solution.  Success  of  this  method  is  dependent  on 
the  ability  to  accurately  measure  or  reconstruct  the  modal  matrix  <i> .  When  using  limited  measurements, 
an  accurate  <i>  matrix  is  not  a  realistic  expectation. 

2.3.3  Connectivity  Based.  To  overcome  some  of  the  shortcomings  of  Baruch’s  method  of 
stiffness  matrix  adjustments,  Kabe  introduced  an  objective  function  which  ensures  stiffness  terms  are 
corrected  in  a  relative  manner  such  that  the  connectivity  of  the  analytical  model  is  preserved.  This 
method  was  also  expanded  to  include  a  weighting  function  which  assigns  a  confidence  level  fo  each 
ferm  in  fhe  sliffness  mafrix.  Kabe’s  algorifhm  uses  a  percenf  change  in  sliffness  value  cosf  funclion 
and  fhen  appends  a  consfrainf  on  fhe  sparsify  paffern  of  fhe  sliffness  mafrix  lo  preserve  fhe  original 
load  palhs.  The  addilion  of  fhe  slruclural  conneclivily  informalion  enables  Kabe’s  melhod  lo  idenlify 
sliffness  changes  exaclly  in  some  cases,  even  when  only  a  limiled  number  of  measuremenls  were  used. 
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A  mathematically  equivalent  but  more  intuitive  method  was  developed  by  Kammer  using  an  alternative 
matrix  minimization  formula.  Similar  work  has  been  investigated  by  Smith  and  Beattie.''’'" 

In  all  of  these  methods,  the  minimization  of  a  matrix  norm  (typically  the  Frobenius  norm)  of 
the  difference  between  the  measured  and  analytical  models  can  be  justified  when  the  goal  is  model 
refinement.  This  method  tends  to  make  small  changes  across  all  the  matrix  elements,  whereas  damage 
to  a  structural  element  will  create  a  large  localized  change. 

23.4  Residual  Force  Based.  Another  approach  which  has  been  applied  to  both  model 
tuning  and  damage  identification  is  the  residual  force  approach.  This  technique,  developed  by  Chen 
and  Garba,  is  based  on  computation  of  a  residual  force  vector  which  represents  the  mismatch  between 
the  analytical  model  and  the  modal  data. For  simplicity,  assume  an  undamped  structure.  The  solution 
to  the  second-order  homogeneous  equation: 

MxFKx  =  Q  (2.40) 

is  given  by  the  eigenvalues  and  eigenvectors  (A,,  <i>j)  in  Equation  (2.31).  If  the  mass  matrix  M  is 
assumed  correct,  and  the  stiffness  matrix  K  is  written  as: 

K  =  Ko  +  AK  (2.41) 

then  substitution  of  Equation  (2.41)  into  Equation  (2.31)  yields: 

AK^i  =  {\iM  -  (2.42) 

The  residual  force  vector  AK^i  is  essentially  equivalent  to  the  modal  force  error  proposed  by 
Ojalvo.  The  residual  force  vector  for  the  mode  can  then  be  written  in  terms  of  the  unknown 
elements  of  the  stiffness  matrix  Akij  and  the  connectivity  matrix  Ci  for  the  mode.  Construction 
of  Ci  is  based  on  nodal  geometry,  elemental  physical  properties,  and  the  measured  eigenvector  for  the 
mode.  An  example  construction  of  C,  is  given  in  Reference  10.  Equating  the  result  to  the  residual 
force  vector  yields: 

AK^i  =  CiAkij  (2.43) 
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The  length  of  the  residual  force  vector  corresponds  to  the  number  of  degrees-of-freedom  of  the  system, 
while  Akij  is  a  vector  whose  length  is  equal  to  the  number  of  independent  elements  in  the  stiffness 
matrix.  A  least  squares  solution  to  Akij  is  found  using: 

Akij  =  CfAK^i  (2.44) 


The  (•)’^  notation  denotes  the  pseudoinverse  operation.  If  multiple  modes  are  used,  the  equations  are 
stacked  so  that: 


Akij 


-cr 

-AK^i- 

^2 

AK^2 

(2.45) 


Advantages  of  this  method  include  its  ability  to  handle  a  subset  of  the  total  number  of  modes  of 
the  system.  However,  when  forming  the  least  squares  solution  without  a  full  set  of  eigendata,  there  is 
no  guarantee  that  a  realistic  Akij  will  result.  This  method  was  demonstrated  using  a  FEM  model  of  a 
78  element  triangular  truss.  The  author  used  3  iterations  and  constrained  Akij  values  to  lie  between 
0  and  100%  stiffness  reduction  from  the  unbroken  values.  The  use  of  a  reduced  length  eigenvector 
(when  the  number  of  sensors  is  less  then  the  degrees  of  freedom)  was  not  addressed. 

The  residual  force  vector  in  Equation  (2.42)  was  also  investigated  by  Kaouk  and  Zimmerman. 
Using  this  method,  the  problem  is  approached  in  a  decoupled  fashion,  using  the  residual  force  vector  to 
localize  the  damage  and  then  using  the  minimum  rank  update  to  determine  the  extent  of  damage.  This 
method  works  well  when  all  the  degrees-of-freedom  can  be  measured,  but  degrades  rapidly  when  the 
size  of  the  measurement  set  is  reduced.  The  determination  of  the  extent  of  damage  varied,  depending  on 
which  modes  were  used  in  the  algorithm.  Best  results  were  obtained  using  only  modes  with  significant 
changes  from  their  nominal  values. 

Additional  work  on  establishing  which  modes  to  use  in  damage  identification  algorithms  has 
been  presented  by  Kashangaki.  In  this  work,  calculation  of  the  strain  energy  associated  with  each 
mode  in  each  member  is  performed  using  the  analytical  model.  It  is  assumed  that  elements  which  are 
highly  strained  in  a  given  mode  will  be  more  readily  identified  using  fhe  eigendafa  for  fhaf  parficular 
mode.  This  informafion  could  easily  be  incorporated  info  Chen’s  work  once  an  assessmenf  of  inifial 
damage  locafion  was  performed. 
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2.3.5  Assigned  Eigenstructure  Based.  An  alternative  approach  to  model  tuning,  more 
frequently  formulated  in  a  controls  setting,  is  through  eigenstructure  assignment.  With  this  technique, 
a  set  of  desired  eigenvalues/eigenvectors  is  achieved  in  closed-loop  manner  hy  the  selection  of  a 
proper  set  of  feedback  gains.  The  following  paragraphs  discuss  some  general  results  for  eigenstructure 
assignment,  along  with  a  closed  form  solution  for  obtaining  the  desired  eigenstructure. 

Consider  a  linear  time-invariant  system  represented  as: 

X  =  Ax  -\-  Bu  (2.46a) 

y  =  Cx  (2.46b) 

u  =  Ky  (2.46c) 

The  dimensions  of  the  state,  control,  and  output  vectors  x,u  and  y  are  of  dimension  ra,  m  and  r 
respectively.  If  the  system  is  both  controllable  and  observable,  and  the  matrices  B  and  C  are  full  rank, 
then  the  following  results  can  be  proven : 

1.  The  position  of  max(m,  r)  closed-loop  eigenvalues  can  be  arbitrarily  assigned,  with  the 
stipulation  that  complex  eigenvalues  must  be  assigned  in  complex  conjugate  pairs. 

2.  A  total  of  min(m,  r)  elements  of  the  closed-loop  eigenvectors  can  be  assigned,  with  the 
same  complex  conjugate  stipulation. 

3.  The  assigned  eigenvector  associated  with  the  assigned  eigenvalue  A,,  must  lie  in  the 
subspace  spanned  by  (A,/  —  A)~^B. 

As  applied  to  structural  damage  identification,  the  matrix  A  is  obtained  from  the  analytical  mass 
and  stiffness  matrices  of  the  tuned  finite  element  model,  as  given  in  Equation  (2.28a).  The  matrix  B, 
which  represents  the  control  influence  matrix,  can  be  constructed  from  the  connectivity  matrix  of  the 
structure.  For  realistic  structures  this  matrix  can  always  be  chosen  to  be  full  rank  from  proper  selection 
of  the  structural  elements.  The  rank  of  the  measurement  matrix  C  will  primarily  depend  on  the  number 
of  sensors  chosen.  As  will  be  shown,  this  will  be  a  critical  factor  in  obtaining  the  desired  eigenstructure 
from  condition  2  above. 
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A  closed-form  solution  to  the  eigenstructure  assignment  problem  has  been  developed  by  Liebst 
and  Garrard.  The  solution  involves  determination  of  a  feedback  gain  matrix  G,  such  that  for  all 
desired  closed-loop  eigenvalue/eigenvector  pairs  (A,,  the  following  relation  holds. 

(A  +  BKC)^i  =  (2.47) 

This  is  equivalently  written  as  finding  the  vectors  Wi  such  that: 

(XJ  -  =  Bwi  (2.48) 

Determination  of  the  gain  matrix  G  once  the  w/s  have  all  been  calculated  is: 

G  =  W[C^-^  (2.49) 

where  W  =  [wiW2  ■  ■  ■]  and  assuming  the  matrix  product  is  non-singular. 

In  general,  the  desired  eigenvectors  may  not  be  achievable.  In  this  case,  it  is  desirable  to  select 
the  Wj  ’s  such  that  a  minimization  between  the  desired  and  achievable  eigenvectors  is  obtained.  For  the 
minimization,  a  weighted  cost  function  is  established  as: 

J,  =  ($,  -  (2.50) 

The  achievable  and  desired  closed-loop  eigenvectors  are  denoted  by  <i>j  and  <i>f  respectively.  The  (•)* 
notation  denotes  the  complex  conjugate  operation.  The  positive  definite  symmetric  weighting  matrix 
Pi  can  be  chosen  to  weight  certain  elements  of  the  vector  of  differences  in  the  desired  and  achievable 
eigenvectors  more  heavily  than  others.  Thus  for  a  given  desired  eigenvalue  A,,  it  is  desirable  to  minimize 
Ji  subject  to  Equation  (2.48).  The  solution  is  found  by  appending  Equation  (2.48)  to  Equation  (2.50) 
with  a  Eagrange  multiplier  //, . 

J,  =  +  //*  [(XJ  -  -  Bw,]  (2.51) 
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The  cost  Ji  is  minimized  in  equation  2.50  when: 


dwi 


-v*B  =  0 


(2.52) 


3  T 


Writing  equations  2.48,  2.52,  and  2.53  in  matrix  form  yields: 


or  equivalently. 


-(Xil-A)  -B  0 

-$r 

-  0  - 

0  0 

Wi 

= 

0 

1 

2a 

o 

AX 

1 

* 

1 _ 

-  U  - 

-$r 

-  0  - 

N. 

Wi 

= 

0 

-  U  - 

The  ohtainahle  A,  and  w,  are  then  given  hy: 


(2.53) 


(2.54) 


(2.55) 


-$r 

-  0  - 

Wi 

=  N-^ 

0 

-  U  - 

(2.56) 


It  can  he  shown  that  iV,  is  always  non-singular  with  positive  definite  Pi  and  B  full  rank,  even  if  A, 
is  not  moved  from  it’s  open-loop  value.  If  an  eigenvalue/vector  pair  (A,,  <i>j)  is  not  to  he  altered, 
setting  the  corresponding  w,  =  0  assures  that  the  open-loop  values  are  retained.  As  developed  above, 
this  algorithm  requires  the  selection  of  a  set  of  desired  eigenvalues/vectors.  Since  in  general  when 
performing  modal  tests  on  a  structure,  only  eigenvector  elements  corresponding  to  the  instrumented 
degrees-of-freedom  are  identified,  fhe  fechnique  as  presenfed  is  nof  suifahle  for  damage  defecfion.  The 
eigenvector  expansion  mefhod  presenfed  in  Equafion  (2.33)  is  nof  suifahle  to  use  wifh  damage  defecfion 
since  if  requires  knowledge  of  fhe  damaged  mass  and  sliffness  mafrices  which  are  unknown. 

Using  fhe  basic  concepfs  of  eigensfrucfure  assignmenf,  a  mefhod  for  correcfing  FEM  models 
using  eigensfrucfure  assignmenf  was  developed  by  Zimmerman  and  Widengren. Their  solufion 
mefhod  is  based  on  a  symmefric  eigensfrucfure  assignmenf  fechnique  in  which  fhe  symmefry  (nof 
necessarily  fhe  connecfivify)  of  fhe  change  in  fhe  sliffness  malrix  is  preserved.  A  slandard  oulpul 
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feedback  assignment  is  used,  but  a  set  of  pseudosensors  C  and  pseudoactuators  B  are  judiciously 
chosen  to  yield  a  symmetric  feedback  matrix.  The  pseudoactuator  matrix  is  chosen  such  that  the 
measured  eigenvectors  lie  in  the  achievable  eigenvector  subspace.  Thus  each  mode  has  it’s  own 
corresponding  pseudoactuator  matrix.  For  a  given  B,  the  corresponding  C  is  found  from  the  solution 
of  a  generalized  algebraic  Riccati  equation.  Since  there  will  exist  a  set  of  non-unique  real  solutions, 
this  set  must  be  evaluated  for  one  which  yields  the  ‘best’  adjusted  stiffness  matrix.  In  the  solution 
technique,  only  necessary  (but  not  sufficient)  conditions  are  used  and  hence  some  solutions  will  be 
asymmetric  and  can  be  immediately  discarded. 

Again,  the  assignment  algorithm  requires  full  length  eigenvectors.  The  eigenvector  elements 
corresponding  to  the  unmeasured  degrees-of-freedom  are  obtained  using  an  optimal  least-squares 
expansion  into  the  achievable  eigenvector  subspace.  An  alternative  to  the  least  squares  expansion 
has  also  been  developed  by  posing  it  as  an  orthogonal  Procrustes  problem  which  yields  similar,  but 
computationally  more  efficient  results. 

Using  the  symmetric  eigenstructure  assignment  technique  does  not  preserve  the  load  paths.  An 
iterative  scheme  was  introduced  to  zero  out  stiffness  matrix  elements  which  were  zero  in  the  analytical 
stiffness  matrix.  The  disadvantage  of  this  method  is  the  requirement  to  solve  the  generalized  algebraic 
Riccati  equation  and  then  sort  through  the  potential  solution  sets  in  an  iterative  fashion. 

An  alternative  to  the  symmetric  eigenstructure  assignment  is  presented  by  Lindner. In  this 
method,  eigenstructure  assignment  is  performed  for  each  element  of  the  truss,  using  the  elemental 
stiffness  value  as  the  pseudocontrol  variable.  Location  of  the  damaged  element  is  based  on  the 
assumption  that  for  each  eigenvalue  and  eigenvector  pair  (A,,  <i>j)  assigned,  only  the  assignment  using 
the  damaged  element  will  consistently  return  the  same  required  stiffness  change  to  assign  each  (A,,  <i>j) 
pair.  Hence  this  algorithm  requires  a  complete  search  over  all  the  elements,  but  in  a  non-iterative 
fashion.  A  full  set  of  sensor  information  (full  length  eigenvectors)  was  assumed  in  this  work.  An 
advantage  of  this  approach  is  that  it  does  not  require  the  inverse  connectivity  problem,  i.e.  finding  fhe 
elemenf(s)  corresponding  fo  fhe  change  in  fhe  sfiffness  mafrix.  This  is  done  inifially  in  fhe  problem 
formulalion.  A  disadvanfage  in  fhe  defecfion  scheme  is  fhaf  an  increasing  number  of  experimenfal 
modes  are  required  fo  heller  assess  which  elemenl  is  fhe  aclual  damaged  elemenl.  The  melhod  is 
also  unsuilable  for  mulliple  elemenl  failures,  since  each  eigenslruclure  assignmenl  uses  only  a  single 
pseudoconlrol  variable. 
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23.6  Realization  Theory  Based.  A  conceptually  straightforward  approach  to  damage 
detection  would  he  to  simply  compare  the  first-order  matrices  of  the  analytical  model  to  those  obtained 
using  an  identification  algorithm.  A  direct  comparison  however,  is  only  possible  when  the  two  systems 
are  represented  using  the  same  state-space  realization.  An  algorithm  developed  by  Alvin  and  Park 
called  the  common  basis-normalized  structural  identification  procedure  performs  this  transformation. 
This  procedure  requires  the  identified  model  fo  be  fhe  same  order  as  fhe  analyfical  model,  which  may 
be  difficull  fo  obfain.  When  smaller  realizafions  are  idenlified,  a  reduced  order  model  is  obfained.  A 
reducfion  of  fhe  analyfical  model  will  desfroy  fhe  connecfivify  informalion  and  make  damage  defecfion 
difficull.  In  subsequenl  work  by  fhe  same  aulhors,  Ihis  lechnique  was  applied  fo  an  eighl-bay  Imss 
slmclure,  inslrumenled  wilh  Ihree  degree-of-freedom  sensors  al  each  node. The  use  of  a  complete 
sel  of  sensors  al  each  node  allowed  a  posl  analysis  fo  determine  individual  elemenl  damage.  This 
melhod  was  also  shown  fo  be  very  sensilive  fo  experimenlal  imprecision. 

As  previously  menlioned,  when  performing  model  reducfion  fhe  sparsily  of  fhe  sliffness  malrix 
and  hence  connecfivify  informalion  is  nol  preserved.  The  advanlage  of  fhe  reduced  model  is  lhal  fhe 
number  of  sensors  can  correspond  direclly  fo  fhe  degrees-of-freedom  of  fhe  model.  A  compromise 
solulion  was  developed  by  Kim  in  which  an  intermediate  sel  of  coordinates  is  chosen,  greater  lhan 
fhe  number  of  measuremenls  bul  less  lhan  fhe  number  in  fhe  analyfical  model.  The  experimenlal 
modes  were  expanded  using  a  modal  expansion  lechnique,  and  fhe  analyfical  model  was  reduced  using 
a  slandard  model  reducfion  algorilhm.  Using  Ibis  intermediate  coordinate  sel,  an  oplimal  malrix  update 
was  performed  as  developed  by  Baruch. Using  Ihis  intermediate  sel,  fhe  localion  of  fhe  damage  can 
be  localized,  bul  can  no  longer  be  determined  fo  wilhin  an  individual  elemenl. 

2.4  Summary 

A  large  colleclion  of  techniques  for  system  idenlificalion,  model  Inning  and  damage  defecfion 
have  been  presenled.  For  system  idenlificalion,  fhe  ERA  melhod  was  singled  oul  as  fhe  melhod 
of  choice  for  obtaining  Ihe  modal  parameters  form  Ihe  measured  response  data.  Common  lo  mosl 
model  Inning  techniques  was  Ihe  minimizalion  of  a  malrix  norm  on  Ihe  perlurbalion  lo  Ihe  sliffness 
malrix  required  lo  malch  experimental  data.  For  damage  deleclion,  a  malrix  norm  minimizalion  is  nol 
necessarily  desirable.  Addilionally,  several  approaches  lo  Ihe  damage  deleclion  problem  have  been 
reviewed.  All  approaches  require  full-lenglh  eigenvectors  (or  conslruclion  Ihereof)  in  Iheir  formulalion. 
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Of  the  methods  presented,  it  is  the  author’s  opinion  that  the  assigned  eigenstructure  based  approach  has 
the  best  potential  to  be  computationally  efficient,  but  must  be  adapted  to  handle  limited  measurement 
data.  Using  the  desirable  aspects  of  each  method,  such  as  load  path  preservation  and  computational 
efficiency,  a  solufion  was  formulafed  in  fhe  subsequenf  chapfers  which  specifically  addresses  damage 
idenlificalion  when  only  a  minimal  number  of  sensors  are  used. 
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III.  Damage  Identification  Methodology 


The  terms  damage  detection  and  damage  identification  appear  frequently  in  the  literature.  In  this 
work,  the  term  detection  applies  to  the  process  of  monitoring,  typically  on  line,  the  measured  response 
of  the  system.  From  the  response,  a  decision  is  made  as  to  whether  or  not  ‘damage’  has  occurred. 
Damage  can  he  defined  as  a  failed  sensor  or  acfuafor,  or  an  acfual  change  in  fhe  dynamics  of  fhe 
sfrucfure.  In  confrasf  fo  defecfion,  sfrucfural  damage  idenlificalion,  as  used  herein,  refers  fo  fhe  off-line 
analysis  of  fhe  measured  response  fo  determine  damage  fo  individual  sfrucfural  elemenfs.  Clearly 
damage  defecfion  musf  occur  fo  begin  fhe  damage  idenfificafion  phase.  Currenf  faulf  defecfion  schemes 
include  moniforing  residuals  in  an  on-line  esfimafor.  The  residual  is  fhe  difference  hefween  fhe 
measured  oufpul  and  fhe  predicted  oufpul  of  an  on-line  esfimafor  such  as  a  Kalman  filter.  A  faulf 
is  indicated  hy  fhe  residual  level  rising  above  some  predetermined  fhreshold.  Differenl  fhreshold 
levels  are  invesfigafed  in  Faifakis.  Furfhermore,  since  fhe  residual  using  a  properly  funed  esfimafor 
should  appear  as  white  noise,  addifional  research  has  focused  on  monitoring  fhe  ‘whiteness’  of  fhe 
residuals.  In  fhis  research  efforl,  if  was  assumed  fhaf  a  suifable  damage  defecfion  algorifhm  is 
available  and  fhus  was  nol  furlher  invesfigafed.  Each  reporfed  fesf  case  began  wifh  fhe  assumpfion  fhaf 
damage  had  already  been  detected.  The  research  focused  on  fhe  idenlificalion  of  damaged  elemenl(s), 
using  fhe  measured  dala. 

Once  damage  has  been  delected,  fhe  off-line  damage  idenlificalion  process  is  inilialed.  The 
idenlificalion  process  can  be  divided  into  four  main  lasks:  (1)  idenlificalion  of  parlial  modal  properlies 
from  measured  dala  of  fhe  nominal  space  sfrucfure,  (2)  adjusling  fhe  FEM  model  to  malch  fhe  measured 
nominal  parlial  dala,  (3)  analyzing  fhe  exlenl  fo  which  sfrucfural  damage  can  be  localized  to  individual 
sfrucfural  elemenfs  using  fhe  measured  dala,  and  finally  (4)  fhe  idenlificalion  of  sfrucfural  damage 
using  measured  parlial  modal  dala  from  a  damaged  space  sfrucfure.  The  melhodology  for  each  lask  is 
overviewed  in  fhe  subsequenl  seclions  of  fhis  chapter,  while  fhe  malhemalical  developmenl  is  deferred 
fo  later  chapters.  The  four  lasks  are  based  on  a  EEM  model  of  fhe  space  sfrucfure  where  sfrucfural 
damping  is  assumed  negligible  and  is  omitted  from  fhe  analysis  for  simplicily.  The  free  vibralion  of 
fhe  undamaged  sfrucfure  is  modeled  as: 


M'x  -\-  Kx  =  D 


(3.1) 


3-1 


with  the  symmetric  mass  and  stiffness  matrices  M  and  K  G  a;  representing  the  state  vector 

and  X  denoting  a  double  time  differentiation  on  x.  The  eigenvalue  and  eigenvector  for  the  mode 
of  Equation  (3.1)  is  given  as  (A,,  <i>j)  whereas  the  measured  eigenvalue  and  partial  eigenvector  for  the 
same  mode  is  represented  as  The  relationship  between  the  n  dimensional  eigenvectors  <i>j  and 

the  s  dimensional  partial  eigenvectors  (fn  is  (fn  =  The  matrix  C  G  maps  the  full  length 

eigenvectors  into  the  partial  eigenvectors  corresponding  to  the  measured  degrees-of-freedom,  (j)=  C^. 
With  minimal  sensor  information  available,  a  natural  cost  function  representing  the  mismatch  between 
the  eigenstructure  of  the  finite  element  model  and  the  measured  eigendata  is: 

^  ( W  -  l)  +  ^  ^  h  -  ^ij)  (3.2) 

i=l  ^  ^  i=l j=l 

where  the  analytical  eigenvalue  for  the  mode  is  denoted  as  A,  and  denotes  the  element  of 
the  eigenvector  from  the  analytical  modal  matrix  <i>.  The  overbar  indicates  a  measured  quantity. 
The  positive  coefficients  a,  and  bij  allow  for  individual  weightings  in  the  objective  function.  The 
summation  upper  limits  r  and  s  represent  the  number  of  eigenvalues/eigenvectors,  and  elements  of  the 
eigenvectors,  respectively,  from  the  measured  data.  With  the  objective  function  defined,  each  fask  can 
now  be  discussed  by  examining  ifs  relafion  fo  Equafion  (3.2). 

3.1  Identification  of  Modal  Parameters  from  Measured  Data 

The  firsf  fask  is  fhe  exfracfion  of  fhe  sysfem  paramefers  (A,,  from  fhe  measured  response 
dafa.  This  is  accomplished  fhrough  fhe  use  of  fhe  ERA  mefhod.  ERA  is  a  lime  domain  approach  based 
on  fhe  system’s  impulse  response  funclions.  On  orbil,  Ihere  are  several  difficullies  associated  wilh 
direclly  measuring  impulse  response  funclions.  If  is  bolh  difficull  lo  apply  fhe  impulse,  and  difficull 
lo  oblain  adequate  signal/noise  levels  wilhoul  imparling  physical  damage.  An  allernalive  mefhod  is 
lo  measure  frequency  response  funclions,  and  Ihen  perform  an  inverse  Eourier  Iransform  lo  oblain 
Ihe  impulse  response  funclions.  The  frequency  response  funclions  are  averaged  before  performing 
Ihe  inverse  operalion  lo  minimize  Ihe  effecl  of  noise  corruplion.  Addilionally,  since  acceleralion 
measuremenls  are  used,  a  lime  inlegralion  (a  division  hy  juj  in  Ihe  frequency  domain)  is  used  lo  oblain 
velocity  measuremenls,  and  a  second  inlegralion  for  displacemenls.  The  ERA  melhod  is  Ihen  applied  lo 
Ihe  impulse  dala  lo  determine  a  discrele-lime  slate-space  realizalion.  This  realizalion  is  Ihen  converted 
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to  a  continuous  model  and  the  modal  properties  determined  via  an  eigenanalysis.  The  identified  modal 
properties  are  used  in  the  evaluation  of  Equation  (3.2)  either  initially  to  tune  the  baseline  FEM  model 
or  subsequently  to  identify  damage. 

3.2  Model  Tuning 

The  objective  of  model  tuning  is  the  determination  of  small  adjustments  to  the  matrices  in 
Equation  (3.1)  such  that  Equation  (3.2)  is  minimized.  An  additional  requirement  is  placed  on  the 
tuning  method  not  to  add  or  delete  load  paths  in  the  model.  This  is  imposed  for  two  reasons.  The  first  is 
so  the  model  will  correlate  well  with  measured  data  of  the  damaged  structure  when  the  model  is  used  to 
simulate  a  damaged  configuration.  The  second  is  to  ensure  that  changes  in  the  matrices  can  be  directly 
related  to  physical  elements  of  the  structure.  The  model  tuning  task  was  accomplished  using  ASTROS- 
ID,  an  automated  multidisciplinary  software  package.  The  method  employed  uses  a  mathematical 
optimization  strategy  to  minimize  deviations  between  measured  and  analytical  modal  frequencies  and 
partial  mode  shapes.  Search  directions  are  determined  based  on  eigenvalue  and  eigenvector  sensitivities 
to  design  variables.  A  mode  tracking  algorithm  is  also  incorporated  to  identify  and  account  for  mode 
switching  during  the  optimization  process.  It  will  be  demonstrated  in  Sections  6.4  and  8.2.4  that 
ASTROS-ID  can  also  be  utilized  for  damage  identification  by  restricting  the  allowable  changes  to  the 
matrices  in  Equation  (3.1),  only  allowing  changes  which  are  consistent  with  structural  damage. 

3.3  Damage  Identification  Using  Assigned  Partial  Eigenstructure 

An  alternative  to  the  sensitivity  based  approach  to  the  minimization  of  Equation  (3.2)  is  achieved 
through  the  use  of  an  assigned  partial  eigenstructure  (APE)  method.  Using  the  identified  eigensfrucfure 
of  fhe  damaged  sfmcfure  (A,,  and  fhe  funed  FEM  model  (M,  K),  fhe  APE  mefhod  is  applied  fo 
defermine  fhe  magnifude  of  fhe  combinafions  of  fhe  ficlilious  acfuafors  required  fo  mafch  fhe  measured 
dafa.  For  damage  idenlificalion,  if  is  assumed  fhaf  fhe  damage  is  confined  fo  changes  in  fhe  sfiffness 
mafrix  {K)  of  fhe  FEM  model  which  are  consisfenf  wifh  fhe  FEM  formulalion.  Addifionally,  only 
decreases  in  fhe  sfiffness  mafrix  are  permiffed.  The  resulfs  of  a  damage  localizafion  analysis  is  used 
fo  esfablish  fhe  inifial  search  space  for  fhe  given  measured  dafa.  The  minimizafion  of  Equafion  (3.2) 
is  accomplished  by  appending  fhe  eigensfrucfure  consfrainf  of  Equafion  (3.1)  along  wifh  fhe  sfrucfural 
consfrainf  fo  form  fhe  Eagrangian,  which  is  fhen  differenlialed  fo  determine  fhe  necessary  condifions. 
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When  only  a  subset  of  the  total  degrees-of-freedom  are  measured,  an  iterative  scheme  is  required  to 
solve  the  necessary  conditions. 

3.4  Sensor  Prioritization  and  Damage  Localization 

With  only  minimal  sensor  information  available,  two  questions  naturally  arise:  (1)  at  which 
locations  should  the  sensors  be  placed,  and  (2)  to  what  extent  can  damage  be  identified  with  the 
selected  sensor  locations?  An  eigenstructure  sensitivity  based  method  is  presented  to  answer  these 
questions.  The  method  presented  is  based  on  examining  the  first-order  partial  eigenstructure  sensitivity 
to  changes  in  the  structural  stiffness  of  each  element  of  a  finite  element  model.  No  a  priori  knowledge 
of  the  damage  location  is  assumed.  Two  aspects  of  the  partial  eigenstructure  sensitivity  are  explored. 
First  is  the  amount  by  which  variations  of  the  elemental  stiffness  values  change  the  measured  partial 
eigenstructure.  Independent  of  the  damage  detection  scheme  used,  elements  which  produce  little  or  no 
change  in  the  measured  data,  and  consequently  in  the  cost  function  J  of  Equation  (3.2),  will  be  difficult 
or  impossible  to  detect  when  damaged.  Second  is  the  direction  of  change  in  the  partial  eigenstructure. 
Elements  which  produce  similar  or  identical  changes  in  the  partial  eigenstructure  will  be  difficult  or 
impossible  to  distinguish  between  when  damaged.  Therefore,  sensor  locations  are  chosen  so  that  the 
change  in  the  measured  partial  eigenstructure  due  to  damage  is  maximized.  Eocalization  of  the  damage 
to  an  element(s)  is  based  on  both  the  amount  and  direction  of  change  to  the  partial  eigendata  for  the 
chosen  sensor  locations. 

3.5  Summary 

This  chapter  provided  an  overview  of  the  four  tasks  associated  with  damage  identification  and 
presented  the  methodology  that  will  be  used  to  accomplish  each  task.  In  the  following  chapters, 
the  mathematical  basis  for  each  task  will  be  presented  along  with  a  description  of  how  this  was 
implemented  in  software.  An  analytical  example  using  a  41 -element  8-bay  planar  free-free  truss  is 
presented  to  illustrate  each  task. 
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IV.  Identification  of  Modal  Parameters  from  Measured  Data 

4.1  Overview 

Independent  of  the  damage  identification  or  model  tuning  method  used  is  the  requirement  to 
measure  modal  data  System  identification  using  ERA  reconstructs  the  transfer  functions 

based  on  identification  of  the  system’s  impulse-response  functions.  A  disadvantage  with  this  method  is 
the  requirement  to  excite  the  system  with  an  impulse.  This  may  he  impractical  to  implement  on  certain 
systems,  and  it  may  he  difficult  to  obtain  response  levels  which  are  significantly  above  the  measurement 
noise  levels.  An  alternative  procedure  is  obtained  using  time-averaged  frequency  response  functions 
to  obtain  the  impulse  response  functions,  as  discussed  next. 

4.2  Obtaining  Measured  Data 

It  is  assumed  that  the  flexible  space  structure  is  equipped  with  at  least  one  disturbance  actuator 
(input)  and  is  instrumented  with  at  least  one  accelerometer  (output).  From  random  vibration  testing,  the 
averaged  frequency  response  functions  are  computed  for  each  input/output  combination.  This  is  done 
by  exciting  the  structure  with  a  band-limited  pseudo-random  noise  sequence  applied  to  each  actuator. 
The  method  below  details  the  computations  for  a  single  applied  input  series  u{t)  and  a  single  measured 
output  series  y{i).  The  extension  to  the  MIMO  case  is  simply  a  matter  of  subscripting  the  input/output 
pairs.  The  definition  of  the  discrete  Fourier  transform  of  each  measured  input  and  output  time  series 
u{t)  and  y{t),  consisting  of  N  points  sampled  with  sample  period  T  is  given  as: 

N-l 

U [k)  =  u(i)e~^  *  N  A;  =  0, 1, . . . ,  N  —  1  (4.1) 

8  =  0 

From  this,  the  auto  spectra  ^uu(k)  and  cross  spectra  ^uyik)  are  computed  using: 

'Suu(k)  =  ^U(kyU(k)  A;  =  0,l,...,N/2- 1  (4.2) 

and 

'^uy(k)  =  ^U(kyY(k)  A:  =  0,l,...,N/2-l  (4.3) 

where  (•)*  represents  the  complex  conjugate  operation  and  U(k)  and  Y (k)  are  the  discrete  Fourier 
transforms  of  the  sampled  input/output  time  series.  Averaged  spectral  estimates  ^uyik),  are  then 
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obtained  from  m  sample  sets  by: 


^  m 

=  A;  =  0,l,...,N/2- 1  (4.4) 

p=i 

The  averaged  frequency  response  functions  (FRF)  denoted  as  Huy{k),  are  then  obtained  from: 

Huy{k)  =  l^^  A:  =  0,l,...,N/2-l  (4.5) 

^nn{k) 

The  impulse  response  function  is  computed  for  each  input/output  combination  from  the  definition 
of  the  inverse  Fourier  transform,  given  as: 

N-l 

h{i)  =^Hy_y{k)e~^‘"~^  i  =  0, 1, . . . ,  N  -  1  (4.6) 

k  =  0 

This  “averaged”  impulse  response  function  can  then  be  used  as  an  input  to  a  time  domain  identification 
algorithm  such  as  the  ERA. 

4. 3  Eigensystem  Realization  Algorithm  ( ERA ) 

ERA  is  based  on  the  singular  value  decomposition  of  the  block  Hankel  matrix.  Consider  the 
MIMO  discrete-time  linear  system  described  by: 

x(k  -f  1)  =  Ax(k)  +  Bu(k)  (4.7a) 

y(^k)  =  Cx{k)  -f  Du{k)  (4.7b) 

where  (  A,  B,  C,  D)  are  of  dimension  (n  X  n),  (n  X  p),  (q  X  n),  and  {q  X  p)  respectively.  A  solution 
to  the  Markov  parameters  *  is  determined  from  a  unit  impulse  response  Y  as: 

Y{k)  =  CA’^-^B  (4.8) 

‘The  discrete  impulse  sequence  is  commonly  referred  to  in  the  literature  as  the  Markov  parameters. 
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The  ERA  method  determines  the  Markov  parameters  hy  forming  the  block  Hankel  matrix  H(k), 
composed  of  the  sampled  unit  impulse  response. 

-  Yik)  •••  Yik  +  n,)  - 

H(k-1)=  :  ■■■  :  (4.9) 

.Y^k  +  rir)  ■■■  Y^kYrir+ric). 

where  and  are  arbitrary  integers  satisfying  the  inequalities  qn^  >  n  and  prig  >  n.  The  singular 
value  decomposition  of  H  evaluated  at  A;  =  1  is  expressed  as: 

H{0)  =  PDQ'^  (4.10) 

The  matrices  P  and  Q  contain  the  left  and  right  singular  vectors  respectively,  and  D  is  a  diagonal 
matrix  containing  the  singular  values.  The  discrete-time  minimum-order  model  is  then  computed  from 
the  decomposition  as: 


(4.11a) 

B,  =  E]!‘^QlEp 

(4.11b) 

(4.11c) 

O 

II 

(4.1  Id) 

The  subscript  n  represents  the  first  n  columns  of  P  and  Q.  The  matrix  D„  is  a  diagonal  matrix 
composed  of  the  n  non-zero  singular  values.  is  [/p,0],  and  Ej  is  [lq,0]  where  Ip  and  Iq  are 
identity  matrices  of  order  p  and  q  respectively  and  0  is  the  zero  matrix.  A  transformation  from  the  z- 
plane  to  the  s-plane  based  on  the  sample  rate  of  the  impulse  data,  can  then  be  performed  if  a  continuous 
model  is  desired.  The  original  development  of  the  ERA  algorithm  is  given  in  Juang  and  Pappa.  The 
condensed  derivation  presented  above  was  based  on  work  by  Crassidas  et  al.  ™  Determining  D„  from 
D  in  Equation  (4.10)  when  using  noise  corrupted  data  may  be  difficult  since  there  are  no  longer  only 
n  non-zero  singular  values.  The  effects  of  noise  on  the  sampled  data  is  presented  in  subsequent  papers 
by  Juang  and  Pappa  ,  and  Akers  and  Bernstein. Alternative  strategies  have  been  investigated  to 
handle  noisy  measurements  and  are  found  in  the  literature, 
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The  ERA  identified  quadruple  {  D^)  is  used  as  the  basis  for  the  identified  eigenstruc- 

ture.  However,  when  comparing  eigenvectors  as  in  Equation  (3.2),  it  is  required  that  the  components 
of  the  eigenvectors  represent  the  same  coordinates  for  a  meaningful  comparison.  Clearly  the  eigen¬ 
vectors  of  the  matrix  will  vary  depending  upon  the  realization  used.  The  ERA  method  produces 
an  equivalent  (from  the  input/output  point  of  view)  hut  different  state-space  realization.  Although  the 
full-length  state  vector  can  never  he  recovered  if  only  limited  sensors  are  used,  the  measured  partial 
eigenvector  can  he  directly  related  to  the  instrumented  degrees-of-freedom  through  the  transformation: 

(4.12) 

The  matrix  is  the  output  matrix  determined  hy  ERA.  The  vector  <i>ci  is  the  eigenvector  of  the 
matrix  A^.  The  matrix  A^  is  the  state-space  A^  matrix  from  ERA  after  a  conversion  from  the  discrete 
to  the  continuous  domain  based  on  the  sample  rate.  The  relationship  between  continuous  and  discrete 
time  is  given  as: 


A,  = 

(4.13a) 

T 

=  / 

(4.13b) 

0 

Cd  =  C, 

(4.13c) 

II 

(4.13d) 

This  assumes  a  zero-order  hold  on  inputs  over  the  sample  period  T.  The  subscripts  c  and  d  denote 
continuous  and  discrete  respectively. 

4.4  Software  Implementation 

When  using  experimental  frequency  response  data  obtained  from  accelerometer  measurements, 
some  pre-processing  of  the  frequency  response  measurements  is  required.  The  first  is  a  double  time 
integration  to  obtain  displacement  measurements.  This  enables  the  identification  of  strictly  proper 
transfer  functions.  This  also  results  in  compatible  data  between  the  measured  and  EEM  data,  since  the 
EEM  eigenvector  data  represents  nodal  displacement  and  rotation  information.  The  time  integration  is 
accomplished  efficiently  by  a  division  of  the  frequency  response  function  by  the  complex  number  ju 
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where  lj  is  the  frequency  in  radians  corresponding  to  each  data  hin.  For  the  double  time  integration, 
the  divisor  is  simply  Next,  to  avoid  the  consequence  of  dividing  hy  zero,  the  first  few  data 

points  in  the  frequency  response  function  are  artificially  sef  fo  zero  and  no  division  is  performed.  This 
is  equivalenf  fo  removing  any  hias  from  fhe  measuremenfs.  The  firsl  few  hins  represenf  fhe  sfeady- 
sfafe  componenf.  Typically,  fhe  accuracy  of  fhe  accelerometers  for  very  low  frequencies  is  poor,  and 
Iherefore  no  valuable  dafa  is  losf  by  zeroing  oul  Ibis  porfion.  Lasfly,  for  fhe  parficular  analyzer  used 
in  fhis  experimenf  (see  Hardware  Secfion  8.3.1),  only  dafa  corresponding  fo  0  -  80%  of  fhe  Nyquisf''^' 
frequency  of  fhe  complex  frequency  response  is  oufpul  for  analysis.  Thus,  before  an  inverse  Fourier 
Iransform  can  be  performed,  fhe  dafa  is  padded  wifh  zeros  (acfing  as  low  pass  filler)  up  fo  fhe  Nyquisl 
frequency,  and  Ihen  reflected  aboul  fhe  real  and  imaginary  axes  fo  oblain  fhe  full  symmelric  speclrum. 
For  example,  dafa  consisling  of  4096  real  sampled  lime  poinls,  resulls  in  2048  unique  complex  poinls. 
Of  fhe  2048  complex  poinls,  only  1601  are  slalislically  reliable  due  fo  aliasing  effecls,  and  Ihus  fhis  is 
all  lhal  is  available  for  analysis.  To  oblain  a  discrele-lime  model  based  on  Ihe  same  sample  rale  as  Ihe 
recorded  dala  il  is  necessary  lo  use  Ihe  full  4096  poinl  symmelric  speclrum  for  use  in  Equalion  (4.6). 
To  achieve  Ihis,  Ihe  1601  complex  poinls  are  expanded  lo  2048  poinls  by  padding  wilh  zeros,  and  Ihen 
Ihe  symmelric  porlion  is  reconslrucled  lo  achieve  Ihe  full  4096  complex  poinls.  This  reconslrucled  sel 
of  4096  poinls  is  used  in  Ihe  compulalion  of  Ihe  impulse  response  funclion  as  given  in  Equalion  (4.6). 
This  is  repealed  for  each  measured  frequency  response  funclion. 

Erom  Ihe  impulse  response  funclions,  Ihe  block  Hankel  malrix  is  conslrucled  according  lo 
Equalion  (4.9),  using  Ihe  resulls  of  Equalion  (4.6).  The  size  of  Ihe  block  Hankel  malrix  is  chosen 
based  on  Ihe  desired  number  of  modes  lo  be  idenlifled  and  is  seleclable  as  a  parameter  in  Ihe  software. 
In  general,  Ihe  block  Hankel  malrix  is  nol  square,  unless  Ihere  are  an  equal  number  of  inpuls  and 
oulpuls.  Eor  compulalional  efficiency  only  Ihe  singular  values  and  vectors  corresponding  to  Ihe 
smaller  dimension  need  be  computed.  The  slate-space  represenlalion  is  found  from  Ihe  singular 
value  decomposilion  using  Equalions  (4.11).  After  a  conversion  from  discrete  to  conlinuous  lime 
using  MATEAB®’s  d2cin  algorilhm,  Ihe  measured  frequencies  and  shapes  are  oblained  from  an 
eigenanalysis  of  Ihe  malrix  and  Equalion  (4.12)  to  oblain  (A,,  cj)-).  The  sequence  of  evenls  required 
for  Ihe  idenliflcalion  of  Ihe  modal  parameters  is  depicted  in  Eigure  4. 1 . 

4.5  An  Analytical  Example 
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Continuous 


Figure  4.1  Sequence  of  events  for  identification  of  modal  parameters. 


Figure  4.2  41 -element  free-free  planar  truss  example  showing  degree-of-freedom  numbering. 
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Figure  4.3  41-element  free-free  planar  truss  example  showing  element  numbering. 
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To  illustrate  the  contribution  of  each  of  the  four  tasks  to  the  damage  identification  problem,  a 
41 -element  8-bay  planar  free-free  truss  was  chosen,  which  is  the  planar  analogue  of  the  Flexible  Truss 
Experiment  described  in  Section  8.3.  The  truss  was  modeled  using  forty-one  rod  elements,  with  two 
translational  degrees-of-freedom  per  node.  The  vertical  and  horizontal  members  were  constructed  of 
aluminum  with  a  cross-sectional  area  of  0.373  square  inches  and  thirty  inches  in  length.  The  diagonal 
members  were  modeled  as  Lexan  (polycarbonite),  with  a  cross-sectional  area  of  0.540  square  inches, 
an  elastic  modulus  20  times  less  stiff  than  aluminum  and  70%  percent  less  dense.  The  density  of  all 
members  was  equally  scaled  to  produce  natural  frequencies  for  the  first  few  flexible  modes  on  the  order 
of  10  Hertz.  This  configuration  was  also  chosen  because  it  exhibits  attributes  common  to  on-orbit 
flexible  structures  which  include:  a  low  fundamental  frequency,  rigid  body  motion,  minimal  structural 
damping,  and  structural  symmetry.  Redundant  diagonal  members  were  included  in  each  bay  to  help 
illustrate  the  concept  of  damage  localization.  Figure  4.2  shows  the  degree-of-freedom  numbering 
scheme  and  Figure  4.3  show  the  element  numbering  scheme.  An  excitation  actuator  was  simulated 
at  degree-of-freedom  #  2  for  vibration  testing.  Four  accelerometers  were  placed  at  degree-of-freedom 
numbers  2,  10,  18,  and  26. 

To  test  the  ERA  method,  a  simulation  of  the  measured  data  was  performed.  The  analytical  model 
was  used  to  generate  four  sampled-data  time  series  at  the  instrumented  degrees-of-freedom  to  a  random 
input  sequence.  A  second  random  sequence  was  added  to  each  output  time  series  to  simulate  the  effect 
of  noise  corruption.  A  metric  of  the  root-mean-square  (RMS)  value  of  the  noise  signal  divided  by  the 
RMS  value  of  the  original  signal  was  used  to  define  fhe  noise  corrupfion  level.  Fourier  Iransforms 
of  fhese  simulafed  signals  were  computed  and  averaged  fo  compute  fhe  averaged  frequency  response 
funclions  as  described  in  Secfion  4.2.  The  dafa  was  sampled  af  a  100  Hz  rate  using  1024  sample  poinfs 
and  25  averages.  One  percenf  of  crifical  damping  for  all  modes  was  included  in  fhe  analyfical  model  fo 
ensure  fhe  response  funclions  were  bounded.  Using  fhese  simulated  frequency  response  funclions,  fhe 
ERA  idenlifled  eigenvalues  and  eigen veclors  were  computed  and  compared  fo  fhe  FEM  resulls.  The 
resulls  are  labulaled  in  Table  4.1  for  fhe  flrsl  five  flexible  modes. 

A  sludy  was  conducled  fo  illuslrale  fhe  effecl  of  bolh  fhe  noise  level  on  fhe  measuremenls 
as  well  as  fhe  selected  size  of  fhe  block  Hankel  malrix.  The  number  of  sample  poinfs  used  in  fhe 
conslruclion  of  fhe  block  Hankel  malrix  for  a  fixed  number  of  inpuls/oulpuls  (1/4)  was  varied,  as  well 
as  fhe  measuremenl  noise  level.  Increasing  fhe  number  of  sample  poinfs  increases  fhe  compulalion 
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Table  4.1  41-element  free-free  planar  truss  eigenvalues. 


True 

ERA  (Hz)t 

Mode  # 

Frequency  (Hz) 

clean 

10%  noise 

I 

5.90 

5.90 

5.90 

2 

10.59 

10.59 

10.59 

3 

15.67 

15.67 

15.67 

4 

20.35 

20.35 

20.34 

5 

25.14 

25.14 

25.13 

fData  presented  for  75  sample  points  used  in  construction  of  the  Hankel  matrix. 

time  required  to  decompose  the  block  Hankel  matrix.  The  results  shown  in  Figure  4.4  illustrate  that, 
for  an  increase  in  noise,  a  larger  number  of  sample  points  is  required  to  achieve  the  same  relative 
error  in  the  measured  eigendata.  The  identification  error  was  the  value  of  J  as  defined  in  Equafion 
(3.2)  wifh  unify  weighfings.  For  fhe  acfual  fesf  configurafion,  fhe  number  of  sensors  available  and  fhe 
noise  level  are  predetermined  quanfifies.  Thus,  in  fhe  deferminafion  of  (A,,  if  is  desirable  fo  use 
as  many  sample  poinfs  as  compufafionally  feasible  in  fhe  consfrucfion  of  fhe  block  Hankel  mafrix.  All 
subsequenf  idenlificalion  using  ERA  will  adopf  fhis  approach. 

In  general,  a  similar  sfudy  can  be  performed  fo  illusfrale  fhe  effecl  of  adding  addifional  inpufs 
or  oufpufs  fo  exfracf  fhe  modal  quanfifies.  This  is  parficularly  imporfanf  when  fhere  are  modes  which 
are  eifher  unconfrollable  or  undefecfable  from  fhe  chosen  inpuf/oufpuf  sef.  For  fhis  analyfical  example 
however,  fhe  firsl  five  flexible  modes  are  bofh  confrollable  and  observable  from  fhe  inpuf/oulpuf  sef, 
and  fhus  no  furlher  ERA  analysis  on  fhis  model  was  performed.  In  a  separate  sfudy,  fhis  soffware  was 
applied  fo  an  experimenfal  large,  lighfly  damped  space  sfrucfure  fo  produce  a  slate-space  represenlalion. 
The  resulls  are  reported  in  Reference  69. 

4.6  Summary 

The  idenliflcalion  of  modal  parameters  from  measured  frequency  response  funclions  was  pre- 
senfed  using  fhe  ERA  mefhod.  Averaging  in  fhe  frequency  domain  was  performed  fo  mifigafe  noise 
effecls.  A  descripfion  of  how  fhis  algorifhm  is  implemented  in  soffware  was  presented  along  wifh  fhe 
resulls  of  fhe  mefhod  applied  fo  an  analyfical  example.  The  idenlifled  (A,  ,  cj)-)  pairs  can  Ihen  be  used 
eifher  fo  lune  fhe  FEM  model  as  described  in  Chapter  V,  or  idenlify  slruclural  damage  as  described  in 
Chapter  VI. 
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V.  Model  Tuning  Using  ASTROS-ID 

5.1  Overview 

Accurate  prediction  and  simulation  of  the  dynamic  behavior  of  large  flexible  space  structures 
require  analytical  models  which  agree  with  measured  data.  Unfortunately,  uncertainty  in  a  finite  element 
model  (FEM)  implies  less  than  perfect  correlation  between  analytical  and  measured  data.  When  the 
disagreement  is  deemed  unacceptable,  it  is  necessary  for  the  design  engineer  to  make  adjustments  to  the 
FEM.  For  large  problems  the  number  of  potential  parameters  to  adjust,  such  as  elemental  areas,  elastic 
moduli,  inertia  moments,  etc.,  quickly  becomes  overwhelming.  Thus  a  systematic  method  is  required 
to  ensure  that  the  adjustments  produce  the  desired  results.  To  this  end,  a  method  is  introduced  which 
poses  a  numerical  optimization  problem.  Namely,  given  a  set  of  measured  eigenvalues  and  partial 
eigenvectors,  determine  the  values  of  selected  physical  parameters  of  the  model  which  minimize  the 
weighted  deviations  from  the  analytical  eigenvalues  and  eigenvectors. 

ASTROS-ID  represents  an  automated  method  of  adjusting  analytical  finite  element  models  to 
measured  data.  The  algorithm  uses  a  mathematical  optimization  strategy  to  minimize  deviations 
between  measured  and  analytical  modal  frequencies  and  partial  mode  shapes.  A  mode  tracking 
algorithm  is  used  to  identify  and  account  for  mode  switching  during  the  optimization  process. 

A  critical  aspect  to  any  model  tuning  algorithm  is  its  practical  implementation.  Key  advantages 
to  the  method  presented  herein  include  its  ability  to  handle  a  small  subset  of  the  total  eigenstructure 
of  the  system  without  using  an  eigenvector  expansion  method,  and  the  ability  to  track  mode  switches 
during  the  tuning  process.  This  optimization  strategy  was  implemented  using  the  Automated  Structural 
Optimization  System  (ASTROS)  software  package,  developed  by  Wright  Eaboratory.  The  present 
research  work  is  an  extension  to  the  work  begun  by  Gibson  on  the  software  modules  ASTROS-ID,  an 
enhancement  to  ASTROS  to  incorporate  system  identification.  The  overall  intent  of  this  software 
enhancement  is  to  enable  the  user  to  input  a  set  of  desired  modal  frequencies  and  partial  modal 
vectors  and  then  iterate  on  a  set  of  structural  parameters  to  minimize  deviations  between  analytical  and 
experimental  measurements. 
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5.2  Theory 


Model  tuning  is  performed  by  minimizing  an  objective  function  based  on  a  weighted  sum  of 
deviations  from  measured  eigenvalues  and  partial  eigenvectors.  The  free  vibration  of  the  structure  is 
modeled  as: 


Mx  +  Kx  =  Q  (3.1) 

Assuming  real  eigenvalues  and  eigenvectors,  the  objective  function  as  previously  defined  is  given  as: 

^  ^  ^  ^  bij  {(^ij  -  ^ij)  (3.2) 

The  summation  upper  limits  r  and  s  represent  the  number  of  eigenvalues/eigenvectors,  and  elements 
of  the  eigenvectors,  respectively,  chosen  to  be  tuned.  The  minimization  of  the  objective  function  J  is 
carried  out  in  an  inner  and  outer  loop  fashion.  The  outer  loop  consists  of  performing  the  eigenanalysis, 
normalizing  and  matching  analytical  and  measured  modes,  calculating  design  parameter  sensitivities, 
updating  design  variables,  and  detecting  solution  convergence.  The  inner  loop  solves  an  approximate 
optimization  problem  after  each  outer-loop  eigenanalysis  using  the  new  design  variable  sensitivity 
information.  The  inner  loop  is  an  approximate  solution  because  the  sensitivity  information  is  used 
to  project  new  values  of  the  analytical  eigenvalues  and  eigenvectors  for  given  changes  in  the  design 
variables  without  recomputing  the  eigenanalysis.  Once  the  approximate  problem  is  solved  by  a  general 
constrained  optimization  method,  control  is  passed  to  the  outer  loop  where  the  variables  are  updated 
and  a  new  eigenanalysis  is  performed.  The  details  of  computing  the  sensitivity  information  as  well  as 
the  mode  normalization  and  tracking  is  discussed  next.  The  problem  considered  herein  assumes  a  real 
eigenanalysis  with  unique  eigenvalues,  and  for  simplicity  is  not  developed  for  the  case  where  damping 
is  included  in  the  finite  element  formulation. 

5.2.1  Eigenvalue  Sensitivity.  A  second-order  Taylor  series  expansion  is  used  to  project  new 
eigenvalue  values  to  given  design  variable  changes  contained  in  the  vector  Ag.  The  expansion  is  given 
as: 

Ai  =  Xgi  +  VXj  Ag  +  -Ag^(V^Xi)Ag  (5.1) 
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where  denotes  the  initial  eigenvalue  resulting  from  the  eigenanalysis.  The  element  of  the 
eigenvalue  gradient  vector  VA,  is  given  as: 


(5.2) 


where  a  comma  denotes  differentiation  with  respect  to  a  design  variable.  The  elements  of  the  Hessian 
matrix  V^A,  are  given  as: 


[Fi,j<^i,k  +  +  Ai 


with  Fij  defined  as: 


Fij  =  Kj  -  \iMj  -  \ijM 


(5.3) 


(5.4) 


The  notation  {•)ij  represents  differentiation  with  respect  to  the  design  variable  of  some  quantity 
(•),  associated  with  the  mode.  The  decision  to  include  a  second-order  approximation,  rather 
than  only  a  first  order,  was  due  to  the  fact  that  the  terms  appearing  in  the  second-order  eigenvalue 
gradient  are  already  computed  when  calculating  the  first-order  eigenvector  gradients.  Equations  (5.2) 
and  (5.3)  include  the  scalar  normalization  term  M^i.  Thus  the  eigenvector  normalization  will 
scale  the  eigenvalue  gradients.  Proper  choice  of  eigenvector  normalization  is  addressed  in  the  next 
section.  Design  sensitivities  of  the  mass  and  stiffness  matrices  M  and  K  are  known  explicitly 
from  the  finite  element  formulation,  and  can  thus  be  computed  using  either  analytical  derivatives  or 
finite-difference  methods.  Eigenvalue  sensitivity  for  each  mode  included  in  the  objective  function  is 
computed  according  to  Equations  (5.2)  and  (5.3).  The  eigenvector  derivatives  in  Equation  (5.3)  can  be 
computed  by  a  modal  expansion''*’'  or  more  efficiently  by  Nelson’s  method'^*"  when  only  a  subset  of 
eigenvectors  are  involved,  as  explained  next. 


5.2.2  Eigenvector  Sensitivity  and  Normalization.  A  first-order  Taylor  series  expansion  is 
used  to  project  the  new  value  of  the  eigenvector  based  on  the  current  eigenvector  value  <i>oj  and  the 
eigenvector  gradient  V<i>j.  The  expansion  is  computed  for  each  degree-of-freedom  in  the  eigenvector 
set.  The  expansion  of  the  term  of  the  eigenvector  for  the  mode  is  given  as: 


4>ik  =  4>oik  +  ^4>lk^{i 


(5.5) 
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The  eigenvector  gradient  is  found  by  first  differentiating  the  eigenvalue  equation  for  each  mode  with 
respect  to  each  design  variable: 

[K  -  \iM]  =  [\ijM  +  -  Kj]  (5.6) 

Since  the  bracketed  term  on  the  left  hand  side  is  necessarily  singular  by  definition,  the  solution  to 
is  found  by  employing  Nelson’s  method.  The  problem  is  solved  by  separating  into  the  sum  of  a 
particular  and  homogeneous  solution  given  as: 

^i,i  —  +  Vij  (5.7) 

Assuming  there  are  no  repeated  roots,  and  [K  —  AjM]  is  an  (n  X  n)  matrix,  then  its  rank  is  ra  —  1. 
Thus  the  homogeneous  solution  is  found  by  eliminating  a  row  and  column  and  then  performing  the 
inverse. 

%  =  [k  -  \m]  [a,- +  \iMj  -  1,-]  (5.8) 

and 

=  [^ij’  •  •  •  5  ^(;-i)i5  O5  •  •  •  7  ^(n-i)i]  (5-9) 

The  (•)  notation  denotes  matrices  reduced  by  one  row  and  column,  or  vectors  reduced  by  one  element. 
Nelson’s  method  removes  the  row  and  column  corresponding  to  the  maximum  entry  in  <i>j.  Equation 
(5.9)  shows  the  expansion  of  to  Vij  for  the  maximum  entry  occurring  in  the  entry.  Note  that 
the  matrix  inverse  in  Equation  (5.8)  need  not  be  explicitly  calculated.  Rather,  the  reduced  matrices 
can  be  used  in  Equation  (5.6)  and  solved  through  matrix  decomposition  followed  by  forward  and 
back  substitution.  To  solve  for  the  unknown  scalar  constant  in  Equation  (5.7),  a  normalization 
constraint  must  be  applied.  The  objective  of  the  optimization  is  to  minimize  differences  between 
measured  and  analytical  modes,  which  clearly  can  only  be  computed  when  the  eigenvectors  are 
normalized  in  the  same  manner.  To  achieve  this  objective,  a  point  normalization  scheme  is  used  in 
which  eigenvectors  are  normalized  such  that  the  degree-of-freedom  with  maximum  amplitude  in  the 
measurement  set  is  set  to  unity.  Presumably,  this  would  also  correspond  to  a  degree-of-freedom  in 
which  there  was  high  measurement  confidence  (i.e.,  high  signal/noise  rafio).  This  is  imporfanf  because 
fhe  analyfic  gradienf  of  fhis  degree-of-freedom  will  be  idenfically  zero.  If  is  nol  pracfical  nor  useful  for 
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optimization  to  normalize  eigenvectors  based  on  the  analytical  model,  such  as  by  mass  normalization, 
since  the  analytical  matrices  are  the  unknowns  in  the  optimization  routine.  However,  mass  normalized 
analytical  eigenvectors  are  useful  in  detecting  mode  switches  as  discussed  in  the  next  section.  After 
each  eigenanalysis  in  the  outer  loop,  the  eigenvectors  used  in  the  objective  function  must  be  point 
normalized  per  the  measurement  data.  Having  chosen  the  degree-of-freedom  normalization  point  for 
each  mode,  the  normalization  constraint  for  each  mode  can  be  expressed  as: 

=  l  (5.10) 

where  the  matrix  5,  contains  only  one  non-zero  entry,  a  one  in  the  row  and  column  corresponding  to 
the  normalization  point  for  the  mode.  Differentiating  the  constraint  Equation  (5.10)  with  respect  to 
each  design  variable  and  substituting  the  result  from  Equation  (5.7)  yields: 


^JSi[cij^i  +  Vij]  =  0 


(5.11) 


Due  to  the  special  form  of  the  matrix  5,  Equation  (5.11)  simplifies  to: 


=  -(V.:) 


ij  Jnorm.pt. 


(5.12) 


which  is  the  negative  of  the  element  in  vector  Vij  corresponding  to  the  normalization  point  for  the 
mode.  The  point  normalized  eigenvectors  must  be  used  in  calculating  the  term  M <i>j  in  Equations 
(5.2)  and  (5.3)  to  ensure  proper  scaling  of  the  eigenvalue  and  eigenvector  sensitivities.  The  eigenvector 
sensitivity  is  computed  for  each  mode  included  in  the  objective  function.  Eor  measured  modes  where 
only  frequency  information  (not  shape)  is  available,  a  first-order  Taylor  series  expansion  is  used  in  lieu 
of  Equation  (5.1). 

The  assumption  of  [K  —  AjM]  having  no  repeated  roots  for  use  in  Equation  (5.8)  can  potentially 
cause  difficulties  for  structures  with  closely  spaced  modes.  However,  no  difficulty  was  observed  for 
the  models  tuned  during  this  research  effort,  even  though  they  contained  closely  spaced  modes.  An 
examination  of  the  effects  of  repeated  roots,  and  the  required  extension  to  Nelson’s  method  to  handle 
this  case,  is  a  recommended  topic  of  future  research. 
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5.2.3  Mode  Switch  Detection.  For  convenience,  the  engineer  typically  will  assign  a  number¬ 
ing  scheme  to  the  set  of  measured  modes  from  an  experimental  analysis.  It  is  imperative  that  the  same 
numbering  scheme  be  employed  when  comparing  measured  modes  to  analytical  modes.  Prior  to  the 
start  of  the  optimization,  this  numbering  scheme  is  most  easily  facilitated  through  the  use  of  computer 
aided  software  capable  of  displaying  mode  shapes.  At  this  point,  the  designer  can  pair  up  measured 
and  analytical  modes.  Once  the  optimization  starts  however,  it  is  necessary  for  the  software  to  track 
the  analytical  mode  sequencing.  As  the  design  variables  change  it  is  highly  likely  that  mode  switching 
will  occur.  Without  mode  tracking,  an  optimization  scheme  would  become  hopelessly  lost  trying  to 
match  bending  modes  to  torsion  modes,  for  example,  and  vice-versa.  Mode  tracking  can  be  performed 
using  cross  orthogonality  checking.  If  the  eigenvectors  are  mass  normalized,  mode  tracking  can  be 
accomplished  by  computing  the  modal  correlation  coefficient  matrix  between  successive  eigenanalysis 
solutions.  The  cross  orthogonality  matrix  is  given  as: 

O  =  (5.13) 

The  modal  matrix  <i>  includes  only  the  eigenvectors  computed  in  the  eigenanalysis.  The  superscript 
n  denotes  the  iteration  number.  Assuming  mass  orthonormalization,  near  unity  values  in  O  indicate 
high  correlation  between  modes.  Mode  tracking  is  accomplished  by  searching  over  O  successively  for 
the  largest  absolute  value.  If  the  entry  is  in  a  diagonal  location,  this  mode  has  not  changed  between 
iterations.  An  absolute  maximum  in  a  non-diagonal  entry  indicates  the  mode  number  switch  between 
two  iterations  by  its  row  and  column  position.  After  each  mode  is  paired,  the  corresponding  row  and 
column  are  eliminated  from  O,  and  the  search  for  the  next  pair  is  accomplished  on  the  reduced  matrix 
until  all  modal  pairs  are  found.  As  discussed  previously,  mode  switch  detection  is  done  in  the  outer 
optimization  loop,  just  prior  to  point  normalization  and  eigensensitivity  calculations. 

5.3  Software  Implementation 

This  model  tuning  technique  was  implemented  using  the  structural  design  software  ASTROS. 
This  software  uses  an  executive  control  sequence  called  MAPOL  (Matrix  Analysis  Problem  Oriented 
Language)  to  develop  the  solution  sequence.  MAPOL  allows  the  user  to  incorporate  custom  software 
modules  into  the  solution  sequence  while  maintaining  full  access  to  the  ASTROS  solution  modules. 
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Figure  5.1  ASTROS-ID  decision  flow. 

For  the  problem  at  hand,  software  modules  written  in  FORTRAN  77  were  developed  (or  modified  from 
Gibson’s  original  work)  for  eigenvector  normalization,  mode  switch  testing,  eigenvalue  and  eigenvector 
sensitivities,  objective  function  evaluation,  and  processing  data  to  and  from  the  optimization  module. 
A  small  section  of  MAPOL  code  is  then  inserted  into  the  standard  ASTROS  solution  sequence  to 
control  the  calling  of  these  new  modules.  Tuning  information  (desired  frequencies,  mode  shapes  and 
the  corresponding  weighting  values)  are  appended  to  the  standard  bulk  data  input.  It  should  be  noted 
that  the  actual  software  implementation  was  less  restrictive  than  that  given  in  Equation  (3.2),  in  that  it 
allowed  for  a  different  number  of  measured  components  for  each  mode  if  desired.  This  is  applicable  in 
situations  where  numerous  frequencies  can  be  measured,  but  only  a  few  shapes.  Additional  information 
on  the  use  of  ASTROS  as  a  multidisciplinary  design  tool  can  be  found  in  References  55  and  5 1 .  The 
decision  chart  for  ASTROS-ID  is  shown  in  Figure  5.1.  A  brief  description  of  modifications  to  the 
original  software  modules  along  with  a  description  of  new  modules  is  contained  in  Appendix  A. 
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5.4  Selection  of  Tuning  Parameters 


Before  the  optimization  can  be  performed,  the  designer  must  choose  a  set  of  design  variables. 
In  ASTROS,  the  design  variable  for  a  beam  element  is  the  beam’s  cross-sectional  area,  because  weight 
minimization  is  typically  the  overall  objective.  In  model  tuning  however,  other  properties  such  as  the 
elastic  moduli,  mass,  or  both  may  be  the  desired  design  variables.  To  accommodate  the  ability  to  tune 
more  than  one  property  of  an  element,  superposition  can  be  used.  To  illustrate  this  method,  consider 
tuning  the  mass  and  elastic  modulus  of  a  beam  between  two  nodes.  Using  two  beam  elements,  one  with 
the  elastic  modulus  set  to  zero  and  a  second  with  the  density  property  set  to  zero,  will  have  the  combined 
effect  of  a  single  beam  element.  However,  now  both  properties  can  be  set  as  design  variables.  Note 
that  adjusting  the  cross-sectional  area  of  a  beam  with  a  zero  elastic  modulus  is  equivalent  to  adjusting 
the  density  property  of  that  beam.  A  similar  relation  is  true  for  the  beam  element  with  zero  mass 
density.  In  a  similar  fashion,  the  torsional  properties  of  the  beam  can  be  varied  independently.  Using 
this  technique  allows  the  user  a  wide  choice  of  design  variables. 

Selection  of  which  design  variables  to  vary  is  dependent  on  the  designer’s  objective  and  is 
problem  dependent.  If  in  the  development  of  the  finite  element  model  some  of  the  more  complicated 
geometries  were  simplified  using  equivalenf  buf  uncerfain  paramefers,  fhen  fhese  paramefers  are  fhe 
nafural  choice  for  fhe  design  variables.  If  however,  fhe  objecfive  is  damage  idenlificalion,  fhen  design 
variables  musf  be  chosen  which  relate  direcfly  fo  damagable  elemenfs.  Furfhermore,  a  design  variable 
may  be  associated  wifh  a  group  of  elemenfs.  For  funing,  if  may  be  desirable  fo  assign  a  single  design 
variable  fo  fhe  elasfic  modulus  of  all  fhe  longifudinal  elemenfs.  In  fhis  way,  any  symmefry  presenf  in 
fhe  sfrucfure  will  be  refained  fhroughoul  fhe  funing  process.  For  damage  idenfificafion,  if  is  required 
fhaf  a  unique  design  variable  be  assigned  fo  each  elemenf,  in  order  fo  isolate  fhe  damage.  In  general,  fhe 
search  space  (number  of  design  variables)  will  be  much  larger  for  fhe  damage  idenlificalion  problem 
lhan  for  model  funing.  A  melhod  of  reducing  fhe  inilial  search  space  is  presenled  in  Chapter  VII  on 
damage  localizalion.  A  discussion  of  fhe  choice  of  design  variables  is  presenled  for  fhe  analylical 
example  below,  and  for  fhe  experimenlal  models  in  fhe  Chapter  VIII. 

5.5  An  Analytical  Example 

To  illuslrale  fhe  use  of  model  funing  using  ASTROS-ID,  a  41-elemenl  free-free  planar  Iruss 
depicted  in  Figures  4.2  and  4.3  was  modeled  wifh  fhe  density  of  all  fhe  elemenfs  al  90%  of  fhe  Irue 
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Table  5.1  Frequency  values  for  the  tuned  41-element  free-free  planar  truss. 


Mode  # 

True  value 

Erequency  (H2 
Inilial  value 

0 

Tuned  value 

1 

5.90 

6.51 

5.90 

2 

10.59 

12.01 

10.59 

3 

15.67 

17.89 

15.67 

4 

20.35 

23.37 

20.35 

5 

25.14 

28.90 

25.14 

value  and  the  elastic  modulus  at  120%  of  the  true  value  for  the  diagonal  elements.  The  measured  data 
( Aj ,  (j)- )  was  taken  to  be  that  of  the  truss  using  the  true  density  and  stiffness  values  (see  Table  4.1).  The 
design  variables  were  chosen  as  follows.  Because  the  uncertainty  in  the  density  of  all  members  was 
assumed  equal,  the  density  of  all  structural  members  was  set  to  a  single  design  variable.  To  account 
for  the  possible  uncertainty  in  the  joint  connections,  the  elastic  modulus  of  the  vertical,  horizontal,  and 
diagonal  members  was  chosen  as  three  separate  design  variables.  For  this  set  of  design  variables,  the 
elastic  modulus  of  the  different  types  of  members  could  be  varied  independently,  but  all  members  of 
the  same  type  (horizontal,  vertical,  diagonal)  would  vary  simultaneously.  Note  that  the  selection  of 
these  four  design  variables  was  made  only  to  illustrate  the  tuning  method,  and  do  not  represent  any  true 
uncertainty  in  this  theoretical  model. 

The  measured  data  consisted  of  the  first  five  flexible  modes  af  fhe  same  four  insfrumenfed 
degrees-of-freedom  (2, 10, 18, 26),  as  used  in  Chapter  IV.  Equal  weighfing  was  placed  on  all  measured 
values  in  Equafion  (3.2).  The  resulfs  of  fhe  funing  are  lisfed  in  Table  5.1  for  fhe  changes  in  fhe  nafural 
frequencies.  The  corresponding  design  variable  changes  are  listed  in  Table  5.2.  Three  oufer-loop 
iferafions  were  required,  which  decreased  fhe  objecfive  funclion  seven  orders  of  magnifude.  To  fhe 
numerical  precision  of  fhe  inpuf  dafa,  fhe  objecfive  funcfion  was  zero.  As  can  be  seen  from  fhe  dafa  in 
Table  5.1,  fhe  frue  values  were  nol  completely  recovered,  allhough  fhe  objecfive  funclion  was  zero.  This 
is  a  resull  of  using  only  parlial  dafa  in  fhe  objecfive  funclion.  The  achieved  solulion  is  nol  necessarily 
unique.  This  is  parlicularly  frue  when  bolh  fhe  mass  and  sliffness  values  are  allowed  lo  vary,  as  was  fhe 
case  for  Ihis  example.  Eor  damage  idenliflcalion,  varialions  will  be  confined  lo  fhe  sliffness  malrix. 
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Table  5.2  Design  variable  values  for  the  tuned  41-element  free-free  planar  truss. 


Design  variable 

True 

Initial 

Tuned 

density  (all) 

1 

.9 

1.02 

horizontal  stiffness 

1 

1 

1.02 

vertical  stiffness 

1 

1 

1.00 

diagonal  stiffness 

1 

1.2 

1.02 

3  iterations  in  46.5  seconds  of  cpu  time. 

5.6  Summary 

An  algorithm  was  developed  to  tune  finite  element  models  to  measured  eigendata.  Careful 
attention  was  placed  on  eigenvector  normalization,  and  how  the  normalization  relates  to  the  computation 
of  the  eigenvalue  and  eigenvector  sensitivities.  Mode  switches  were  successfully  tracked  through  an 
analysis  of  the  correlation  coefficient  matrix.  Implementation  of  this  algorithm  in  ASTROS  makes 
this  tuning  procedure  readily  accessible  to  the  design  engineer  during  the  finite  element  design  phase. 
This  level  of  experimentally  validated  automated  tuning  represents  a  substantial  improvement  over 
current  capabilities.  Natural  applications  of  this  algorithm  include  the  identification  of  uncertain 
structural  parameters,  as  well  as  the  identification  of  damaged  structural  elements.  This  technique  was 
successfully  applied  to  an  analytical  model.  Experimental  results  using  ASTROS-ID  for  both  model 
tuning  and  damage  identification  are  contained  in  Chapter  VIII. 

ASTROS-ID  requires  an  eigenanalysis  and  gradient  evaluations  at  each  outer-loop  iteration, 
which  can  be  computationally  expensive.  An  alternative  method  to  minimizing  Equation  (3.2)  is 
pursued  using  an  assigned  eigenstructure  method  developed  in  the  next  chapter.  Reduction  of  the 
initial  search  space,  which  will  ease  the  computational  burden,  is  presented  in  Chapter  VII. 
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VI.  Damage  Identification  Using  Assigned  Partial  Eigenstructure 

6.1  Overview 

Presented  is  an  algorithm  to  identify  individual  damaged  structural  elements  of  large  flexible 
space  structures  using  on-orbit  data.  Using  measured  partial  eigendata,  a  control  theoretic  approach  is 
applied  in  which  fictitious  actuators  corresponding  to  each  structural  element  are  assumed.  Using  the 
identified  parfial  eigensfrucfure  ( A, ,  cj)- )  and  fhese  flcfifious  acfuafors,  an  assigned  parfial  eigensfrucfure 
(APE)  fechnique  is  employed  on  fhe  analyfical  model  of  fhe  undamaged  sfrucfure.  An  idenfiflcafion 
of  sfrucfural  damage  is  obfained  direcfly  from  fhe  required  confrol  of  fhe  flcfifious  acfuafors  (sfiffness 
adjusfmenf)  fo  achieve  fhe  desired  eigensfrucfure.  This  mefhod  represenfs  an  alfernafive  fo  fhe  gradienf 
search  fechnique  using  ASTROS-ID  presenfed  in  Chapfer  V. 

6.2  Theory 

Wifh  minimal  sensor  informalion  available,  a  nafural  cosf  funclion  represenfing  fhe  mismafch 
befween  fhe  eigensfrucfure  of  fhe  flnife  elemenf  model  and  fhe  measured  eigendafa  is: 

^  ^  a,  -  l)  +  ^  ^  hj  Uij  -  ^ij)  (3.2) 

i=i  AAj-  /  i=i  j=i 

wifh  all  quanfifies  as  previously  defined.  A  minimizafion  of  Equafion  (3.2)  was  developed  in  Chapter  V 
fo  perform  model  updafing  using  ASTROS-ID,  in  which  fhe  minimizafion  is  solved  using  fhe  eigenvalue 
and  eigenvector  sensifivifies  af  each  iferafion  step.  An  alfernafive  fo  compufing  fhese  sensifivifies  (which 
sfill  requires  an  eigenanalysis  af  each  iferafion  step)  is  accomplished  using  APE,  developed  below. 

The  APE  mefhod  is  based  on  minimizing  fhe  same  cosf  funcfion  given  in  Equafion  (3.2).  Two 
inifial  assumpfions  are  made.  Eirsf,  sfrucfural  damage  is  confined  only  to  changes  in  fhe  sfiffness 
properfies  of  fhe  sfrucfure.  Second,  sfrucfural  damping  is  negligible.  These  fwo  assumpfions  are 
consisfenf  wifh  mosf  on-orbif  damage  scenarios  of  large  flexible  space  sfrucfures.  The  free  vibrafion 
of  fhe  undamaged  sfrucfure  is  modeled  as: 


M'x  -\-  Kx  =  6 


(3.1) 
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with  the  symmetric  mass  and  stiffness  matrices  M,  K  G  a;  representing  the  state  vector  and 

(•)  denoting  a  double  time  differentiation.  With  damage  confined  to  the  stiffness  matrix,  the  damaged 
structure  is  modeled  as: 

Mx  +  (K  —  AK)x  =  0  (6.1) 

where  AK  represents  an  unknown  perturbation  to  the  stiffness  as  the  result  of  structural  damage. 
The  eigenvalue  and  eigenvector  for  the  mode  of  Equation  (6.1)  is  given  as  (A,,  <i>j)  whereas 
the  measured  eigenvalue  and  partial  eigenvector  for  the  same  mode  is  represented  as  (A,,  (f)-}.  The 
relationship  between  the  n  dimensional  eigenvector  <i>j  and  the  partial  eigenvector  (fn  is  (fn  = 

The  matrix  C  G  maps  the  full  length  eigenvectors  into  the  partial  eigenvectors  corresponding  to 
the  measured  degrees-of-freedom.  For  the  APE  method,  the  cost  function  in  Equation  (3.2)  based  on 
the  errors  between  the  finite  element  model  and  the  measured  modes  is  rewritten  in  vector  notation  and 
is  given  as: 

J  =  ^(A  -  Af  A(A  -  A)  +  i  -  Ji)  (6.2) 

i  =  l 

where  A  =  [Ai,  A2, . . . ,  A^]  for  the  r  measured  modes,  with  s  degrees-of-freedom  measured  for  each 
mode.  The  positive-definite  constant  matrices  A  and  W  can  be  used  to  weight  the  contribution  of  each 
term  in  the  overall  cost  function.  This  cost  function  is  then  minimized  subject  to  the  eigenstructure 
constraint: 

(-AiM  +  K  -  AK)^i  =  0  (6.3) 

where  only  r  of  the  (A,,  <i>j)  coupled  with  only  s  components  of  <i>j  are  able  to  be  measured.  Addi¬ 
tionally,  to  ensure  AK  is  consistent  with  the  finite  element  formulation,  the  structural  constraint  is 
represented  as: 

AK  =  BGB'^  (6.4) 

where  B  is  constructed  from  the  nodal  connectivity  information  and  the  elemental  parameters,  and  G 
is  a  diagonal  matrix  composed  of  the  fraction  of  damage  for  each  element.  The  parameterization  of 
AK  is  best  illustrated  using  a  simple  example.  For  a  two  degree-of-freedom  spring  mass  system  as 
shown  in  Figure  6.1,  the  equations-of-motion  obtained  from  an  application  of  Newton’s  second  law 
are: 

niiXi  =  —kiXi  —  k^ixi  —  X2)  (6.5) 
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m2X2  =  -k2{x2  -  Xi) 

These  equations  can  then  he  written  in  matrix  form  as: 


'  nil 

0  ■ 

X  -f 

ki  +  k2 

-k2' 

0 

m2 . 

.  -k2 

k2  . 

a;  =  0 


where  the  stiffness  matrix  K  can  then  he  written  as: 


^it^2  ”^2]  [1  r^i  01  [1 


—  k2  k2 


0  10  A;,  0  1 


The  matrix  pre-multiplying  the  diagonal  matrix  in  Equation  (6.8)  contains  the  structural  connectivity 
information.  For  example,  the  first  column  corresponds  to  spring  ki,  which  is  connected  only  to 
degree-of-freedom  Xi,  and  hence  has  only  a  single  non-zero  entry  in  row  one.  The  second  column 
corresponds  to  spring  k2,  which  is  connected  to  hoth  degree-of-freedom  Xi  and  X2,  and  hence  hoth 
rows  have  non-zero  entries.  Direction  cosines  are  used  to  determine  the  values  of  the  non-zero  entries, 
based  on  the  relative  position  of  the  nodes.  For  this  linear  example,  the  direction  cosines  are  plus  and 
minus  one.  In  the  same  fashion,  the  stiffness  perturbation  matrix  Alf  can  be  expressed  as: 

' Vh  \gi  01  \^/Vl 

Alf  =  =  BGB^  (6.9) 

.  0  J  L  0  £(2  J  L  0  _ 
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Note  that  with  0  <  <  1,  any  comhination  of  decreases  in  the  spring  elements  can  he  modeled  in 

AK.  For  a  truss  structure  constructed  from  p  rod  elements,  B  is  written  as: 


B  =  [B,---Bp]  (6.10) 

A- E- 

Bi  =  J  — ^  •  [0, . .  .,0,Ci,C2,C3,0, . .  .,0,  -Cl,  -C2,  -C3,0, . .  .,0]^  (6.11) 

G  =  diag(pi  •  • -pp)  0  <  Pi  <  1  (6.12) 

with  Cl ,  C2 ,  C3  representing  the  direction  cosines  for  the  element,  inserted  at  the  degrees-of-freedom 
associated  with  the  element.  The  variables  A,,  Ei  and  Li  are  the  cross-sectional  area,  elastic 
modulus,  and  length  of  the  element  respectively.  A  value  of  pi=  0  corresponds  to  an  undamaged 
element  whereas,  gi  =  1  corresponds  to  a  complete  loss  of  stiffness  to  the  element.  For  heam 
elements  with  six  degrees-of-freedom  per  node,  the  expression  for  Bi  becomes: 


Bi  =  Rt 


_2_ 

L, 


-1 


2_ 

L, 


-1 


-2 

L, 


-2 

L, 


-1 


-1  -1 


-1 


•diag^^(Ai  3/ii  3/2*-  Pij)  In  hi 


1/2 


(6.13) 


where  Ri  is  the  rotation  matrix  between  the  element’s  local  coordinate  frame  in  which  the  inertia 
properties  h  and  h  are  defined,  and  the  global  coordinate  system.  The  variables  p,  and  j)  are  Poisson’s 
ratio  and  the  torsional  stiffness  of  the  element  respectively.  Only  the  non-zero  portion  of  Bi  is 
shown,  which  occurs  in  the  rows  corresponding  to  the  global  degrees-of-freedom  associated  with  the 
element.  For  beam  elements,  G  is  now  a  block  diagonal  matrix,  with  each  diagonal  block  linked 
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to  a  single  design  variable  gi  (gi  ■  Iq,  where  Iq  is  the  6x6  identity  matrix).  Although  developed 
explicitly  herein  only  for  rod  and  beam  elements,  any  element’s  symmetric  matrix  Alf,  can  be  written 
as  BGiB^  using  the  non-zero  singular  values  and  singular  vectors  of  Alf,.  The  minimization  of 
the  cost  function  in  Equation  (6.2)  is  solved  by  forming  the  Lagrangian  and  establishing  necessary 
conditions.  The  appended  cost  function  J,  for  the  mode  is: 

Ji  =  ~  ~  ~  4'i)  +  li  —  BGB^  (6.14) 

where  //,  is  a  vector  of  Lagrange  multipliers  for  the  mode.  With  (fn  =  this  becomes: 

Y  =  -  Tif  +  -  ¥7)  +  +  K  -  BGB^)^,  (6.15) 

It  is  then  assumed:  3  6?  9  A,  =  A,  ,  V  )  =  1, . . . ,  r  implying  there  are  sufficient  design  variables 
(structural  elements)  to  achieve  the  measured  eigenvalues.  This  is  satisfied  if  fhe  measured  dafa 
is  consisfenf  wifh  acfual  damage.  For  fhe  case  of  noise  corrupfed  measuremenfs.  A,  ~  A,  ,  V  f  = 
1, .  . . ,  r,  and  is  assumed  fo  confribufe  negligibly  fo  fhe  cosf  funclion.  Wifh  fhese  assumpfions,  Equafion 
(6.15)  reduces  fo: 

J,  =  -  ¥7)  +  vj {-\M  +  K  -  BGB'^)^,  (6.16) 


The  necessary  condifions  for  fhe  minimizafion  become: 


^  =  (-A,M  +  K  -  BGB'^)^,  =  0 

avi 


dY 

dG 


d 


{vjBGB^^i)  =  0 


C^WiC{^i  -  -f  lufi-XiM  +  K  -  BGB'^)]^  =  0 


which  can  be  rewriffen  as: 


-\M  +  K)^i  -  BGB'^^i  =  0 


(6.17) 

(6.18) 

(6.19) 

(6.20) 
(6.21) 
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(6.22) 


+  (-A,M  +  K)v,  -  BGB'^v,  = 

To  pose  this  non-linear  optimization  problem  as  an  approximate  linear  problem,  it  is  necessary  to 
introduce  the  matrix  operator  P(a,  (3).  In  the  case  where  there  is  only  one  column  of  the  matrix  B 
associated  with  each  design  variable,  such  as  for  spring  and  rod  elements  as  given  in  Equations  (6.10 
through  6.12),  P(a,  (3)  is  defined  as: 

n 

P(a,(3)  with  P  and  (3  G  ,  a  G  P^j  =  (6.23) 

k  =  l 

where  is  the  row  and  the  column  of  the  matrix  P.  In  terms  of  the  operator  P,  terms  of  the 
form  (3T (3'^ a  can  be  written  as: 

f3Tf3’^a  =P(a,f3)^  where  7  =  diag(r)  ,  F  G  3?^^^(diagonal)  (6.24) 

When  there  are  multiple  columns  of  B  associated  with  a  single  design  variable,  such  as  given  in  Equation 
(6.13)  for  beam  elements,  an  additional  summation  of  the  columns  of  P  is  required  for  each  design 
variable.  In  the  case  where  there  are  g  columns  of  P  for  each  design  variable,  P  G  G  G 

and  hence  P  G  according  to  Equation  (6.23).  The  summation  is  then  defined  as: 

3  -  2g  _  pq  ^ 

P  =  [J:Pj  Z  Pj  ■■■  E  Pj]  ,  (6.25) 

i  =  l  i=3  +  l  j=(p-l)q  +  l 

where  Pj  is  fhe  j column  of  P.  Nofe  fhaf  P  =  P  for  fhe  case  where  g  =  1 .  Eurfhermore,  consfrucfing 
P  from  P  gives  fhe  design  engineer  fhe  abilify  fo  link  mulfiple  elemenfs  fo  a  single  design  variable 
gi  as  desired.  In  terms  of  fhe  operator  P,  fhe  following  subslifufions  can  be  made  in  fhe  necessary 
condifions: 

BGB'^^i  =  P{^i,B)g  (6.26) 

^  (//f  PGP^^,)  =  P(<i>,,  Pf //,.  (6.27) 

BGB^Vi  =  P{vi^B)g  (6.28) 

The  necessary  condifions,  wifh  fhese  subsfifufions  become: 

{-\M  +  K)^i  -  P($i,  B)g  =  0  (6.29) 
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=  0 


(6.30) 


+  {-\M  +  K)v,  -  P{v„B)g  =  (6.31) 

which  can  now  he  written  in  matrix  form  for  the  mode  as: 


-(-\M  +  K) 

0  -P($„P)- 

0 

0 

0 

Vi 

= 

0 

.  C^WiC 

(_A,M  +  P')  -P(//„P). 

-  9  - 

Then  since  g,  the  vector  of  fractional  structural  damage  for  each  element  {gi)  is  the  same  for  each 
measured  mode,  the  necessary  conditions  can  he  assembled  as: 


N 


r$ii 

Vr 

L  g  J 


0 

0 


0 
0 


(6.33) 


where 


■(-AiM  +  lf) 

P(<i>i,  pf 

-P(^uB)- 

C'^WiC 

(-AiM  +  K) 

-P(vuB) 

{-\M  P  K) 

-P($,,P) 

P(<i>.,  Pf 

C^WrC  (-KM  +  K) 

- 1 

1 

(6.34) 

and  Ne  §^^(2>.+p)  rx(2«r+p)  ^  representing  an  over-determined  set  of  linear  equations  whenever  r  >  1. 
Only  the  non-zero  entries  are  shown.  The  desired  solution  vector  g  is  then  found  from  a  least  squares 
solution  to  Equation  (6.33)  using  a  decomposition  and  back  substitution.  Since  iV  =  , //,  ), 

an  iterative  scheme  is  introduced  updating  (<i>j ,  //,  )  with  the  results  of  the  previous  iteration.  The  initial 
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guess  at  Vi)  is  to  use  the  nominal  vector  <i>j  from  the  undamaged  model,  with  the  measured  cjyi 
elements  inserted  at  the  measured  degrees-of-freedom.  It  is  assumed  that  structural  damage  did  not 
result  in  catastrophic  damage,  and  thus  the  nominal  eigenvectors  are  a  reasonable  initial  guess.  Some 
typical  deviations  in  the  eigenvectors  resulting  from  damage  to  a  small  number  of  structural  elements 
are  shown  in  Section  6.4.  Vector  normalization  and  sign  convention  is  accounted  for  by  setting 
—  |^z  |  —  ^  and  *  (j)i  ^  0.  The  Lagrange  multiplier  vector  Vi  is  initially  assumed 

to  be  zero.  During  the  iteration  process,  values  of  the  damage  fraction  gi  outside  the  allowable  range 
are  removed  from  subsequent  iterations,  further  reducing  the  parameter  search  space.  The  weighting 
matrix  W,  is  nominally  set  to  the  identity  matrix  (scaled  such  that  1 1 W,  1 1  ~  1 1  A'  1 1 )  corresponding  to  the 
case  where  all  measurements  are  assumed  at  the  same  level  of  uncertainty. 


Having  presented  both  the  ASTROS-ID  and  the  APE  methodology  for  the  minimization  of 
Equation  (3.2),  a  comparison  can  now  be  made,  illustrating  the  different  nature  of  the  two  solution 
techniques.  Because  the  APE  method  assumes  the  eigenvalues  can  be  achieved  exactly,  it  is  equivalent 
to  placing  a  very  high  weighting  on  the  eigenvalues.  Additional  design  degrees-of-freedom  are  then 
used  to  achieve  the  desired  eigenvector  components.  However,  if  the  weighting  on  the  eigenvalues  in 
ASTROS-ID  is  too  large  as  compared  to  the  eigenvector  weightings,  then  the  eigenvector  information 
is  negligable  in  the  cost  function,  and  it  is  equivalent  to  only  using  the  eigenvalue  information.  This 
tends  to  spread  the  assigned  damage  over  several  elements.  Eurthermore,  since  their  are  no  constraints 
placed  on  the  solution  vector  in  the  APE  method  until  after  the  decomposition  and  back  substitution  is 
formed,  unrealistic  stiffness  changes  often  result,  which  are  then  easily  identified  and  removed  from  the 
search  space.  This  establishes  a  natural  and  simple  method  to  reduce  the  search  space  for  subsequent 
iterations.  No  such  search  space  reduction  is  currently  available  for  ASTROS-ID,  and  is  listed  as  a 
topic  for  future  research.  Note  that  in  either  case,  convergence  to  the  global  minimum  is  not  guaranteed. 

Solution  techniques  used  to  minimize  Equation  (3.2)  can  be  classified  as  eifher  a  mafhemafical 
programming  approach  or  an  opfimalify  criferion  approach.  Mafhemafical  programming,  also  referred 
fo  as  a  direcf  search  mefhod,  is  based  on  esfablishing  search  direcfions  fo  arrive  af  fhe  solufion.  This 
technique  is  used  in  ASTROS-ID.  In  confrasf,  opfimalify  criteria  mefhods  seek  fo  direcfly  safisfy  fhe 
opfimalify  criteria,  iferafively  using  a  recursion  formula.  Salisfaclion  of  Equafion  (6.33)  is  fhe  opfimalify 
criferion  used  in  APE.  An  overview  of  fhe  fwo  approaches  as  applied  fo  sfrucfural  opfimizafion  can  be 
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Figure  6.2  Assigned  partial  eigenstructure  decision  flow. 

found  in  the  text  by  Haftka  and  Giirdal.  Additional  work  using  optimality  criteria  is  presented  in  the 
work  by  Venkayya. 

6.3  Software  Implementation 

Structural  damage  detection  using  APE  proceeds  as  follows.  First,  the  matrix  B  is  constructed 
according  to  Equation  (6. 1 1)  for  rod  elements,  or  Equation  (6. 13)  for  beam  elements.  A  Bi  is  constructed 
for  each  element  in  the  search  space.  Next,  using  the  mass  and  stiffness  matrices  M  and  K  from  a  tuned 
finite  element  model  (one  in  which  the  measured  data  of  the  undamaged  structure  correlates  well  with 
the  model)  and  the  measured  eigendata  of  the  damaged  structure.  Equation  (6.33)  is  constructed.  A  QR 
decomposition  and  back-substitution  is  used  to  solve  for  the  achievable  eigenvectors  <i>j,  the  Eagrange 
multipliers  //,  ,  and  the  damage  fractions  pi.  Elements  in  g  outside  an  allowable  range  (i.e.,  an  increase 
in  stiffness  or  negative  stiffness)  are  removed  from  the  search  space  by  removing  the  corresponding 
columns  of  B.  At  each  iteration  step,  the  matrix  N  is  updates  using  either  the  new  ,  //,  )  solution 
pair  or  the  new  reduced  B  matrix.  The  solution  sequence  is  repeated  until  convergence.  The  decision 
chart  for  this  algorithm  is  shown  in  Figure  6.2. 
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For  the  initial  construction  of  N  using  the  measured  (A,,  cj)-)  data,  two  preprocessing  steps  are 
required.  The  first  is  a  pairing  of  the  measured  mode  shapes  with  the  nominal  mode  shapes.  This  is 
accomplished  hy  normalizing  the  oo-norm  of  the  partial  mode  shapes  to  unity  and  then  checking  the 
cross-orthogonality  relation  as  given  in  the  following: 

O  =  ,  (^and(^  G  (6.35) 


Modes  are  paired  based  on  the  maximum  values  in  the  row  and  column  positions  of  the  matrix  O .  A 
second  check  is  performed  to  verify  the  measured  modal  frequencies  of  the  damaged  structure  are  at 
or  helow  the  corresponding  frequencies  of  the  undamaged  analytical  model,  i.e..  A,  <  A,  V  i.  For 
each  mode,  measured  frequencies  above  the  analytical  frequencies  are  set  to  the  analytical  values. 
Structural  damage,  when  confined  fo  decreases  in  fhe  sfiffness  mafrix,  can  only  decrease  fhe  nafural 
frequencies.  This  requiremenf  can  easily  be  shown  by  esfablishing  fhe  negafive  definileness  of  changes 
in  fhe  eigenvalues  fo  changes  in  fhe  sfiffness  mafrix.  As  will  be  developed  in  Chapter  VII,  fhese 
changes  fo  firsl  order  can  be  expressed  as: 

=  (6.36) 

For  fhis  firsl-order  model,  if  is  sufficienl  fo  show  fhe  negafive  definifeness  of  fhe  eigenvalue  wifh 
respecf  fo  fhe  sfiffness  value  change.  This  relafionship  is  given  as: 


dgj 


(6.37) 


which,  using  fhe  definilion  of  P  from  Equafion  (6.23),  valid  for  rod  elemenfs,  and  carrying  oul  fhe 
vecfor  mulfiplicafion  yields: 

dX  ^  ^ 

^  ^ (6.38) 

%  titi 

Rearranging  fhe  summafion  yields: 

dX-  "  " 

^  ^  ^kiBkj  (6.39) 

;=i  k  =  l 
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which  simplifies  to: 

dXi  / "  \ ' 

^  =  (6.40, 

which  is  a  negative  definite  quantity.  If  the  more  general  definition  of  P  is  used  as  given  in  Equation 
(6.25),  an  additional  summation  is  required  over  the  corresponding  columns  of  B  associated  with  the 
design  variable.  Again,  this  is  a  negative  definite  quantity  since  it  is  the  sum  of  negative  definite 
terms  from  Equation  (6.40). 

The  APE  iterative  solution  technique  involves  the  least  squares  solution  of  the  matrix  N  in 
Equation  (6.34).  This  matrix  is  classified  as  a  large  sparse  rectangular  matrix  which  is  possibly 
singular,  and  therefore  the  solution  technique  is  tailored  to  this  matrix  classification.  A  sparse  QR 
algorithm  is  used  to  decompose  N  and  then  a  back-substitution  to  compute  the  solution  vector.  Due 
to  the  size  of  iV,  a  sparse  solver  is  required  for  all  but  small  pedagogical  problems.  Column  pivoting 
in  the  QR  decomposition  is  not  incorporated  since  pivoting  does  not  preserve  the  matrix  sparsity  and 
increases  the  number  of  fill-ins.  Eor  the  case  when  N  is  singular,  the  nature  of  the  solution  is  different 
than  that  of  the  Moore-Penrose  pseudo-inverse.  The  solution  vector  contains  as  many  zero  entries  as  the 
rank  deficiency  of  N .  This  is  a  desirable  attribute  when  determining  damage  detection,  since  typically 
damage  is  localized  in  the  structure,  and  hence  only  a  small  number  of  elements  have  non-zero  damage 
fractions.  The  APE  algorithm  was  coded  using  MATEAB®,  with  portions  written  in  Eortran  to  handle 
the  sparse  matrix  manipulations  to  speed  processing  time.  MATEAB’s  spqrmex  algorithm  was  used 
to  perform  the  decomposition.  In  practice,  the  range  of  allowable  damage  fractions  (0  <  <  1)  was 

widened  to  ensured  ^j  ’s  were  not  discarded  prematurely  before  convergence. 

An  additional  consideration  in  damage  detection  algorithms  is  that  of  uniqueness.  With  only 
partial  modal  data  available,  the  problem  is  generally  ill-posed  and  hence  there  may  exist  multiple 
damage  fractions  which  result  in  the  same  partial  eigendata.  This  problem  is  exacerbated  when 
noise  corrupted  measurements  are  used.  The  use  of  the  structural  constraint  as  given  in  Equation 
(6.4),  combined  with  the  restriction  that  the  damage  fraction  g  lie  within  an  allowable  range,  helps  to 
minimize  the  problem  of  non-uniqueness.  Eor  a  given  problem,  solution  uniqueness  is  dependent  on 
the  number  of  modes,  the  number  of  eigenvector  components  measured,  the  quality  (signal/noise)  of 
the  measurements,  the  finite  element  model,  and  the  numerical  accuracies  of  the  detection  algorithm. 
To  account  for  non-unique  solutions  in  the  APE  algorithm,  an  off-line  sensitivity  analysis  is  performed 
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Figure  6.3  41-element  free-free  planar  truss  showing  element  and  degree-of-freedom  numbering. 


once  the  size  and  quality  of  the  measurement  set  has  been  determined.  Using  the  results  of  this 
analysis,  the  parameter  search  space  is  confined  to  elements  which  produce  unique  and  identifiable 
eigensfrucfures.  Elemenfs  whose  damage  resulfs  in  an  idenlifiable  change  in  fhe  measured  parfial 
eigensfrucfure  are  defined  as  defecfable.  Nexf,  APE  symmefric  elemenfs  are  defined  as  damaged 
sfrucfural  elemenfs  which  produce  fhe  same  measured  parfial  eigensfrucfure.  The  resulfs  of  fhe  APE 
algorifhm  indicafe  damage  fracfions  which  are  confined  fo  an  APE  symmefric  elemenf  sef.  Eurfher 
refinemenf  fo  an  individual  elemenf  wifhin  fhe  sef  is  nof  possible  wifhouf  addifional  or  higher  qualify 
measuremenf  dafa. 

Numerous  fesfs  of  fhe  APE  algorifhm  were  conducfed  using  finife  elemenf  models  of  sizes 
ranging  befween  16  and  192  degrees-of-freedom  consfrucfed  of  8  fo  104  sfrucfural  elemenfs  using 
eifher  spring,  rod,  or  beam  elemenfs.  In  all  cases,  when  using  analyfical  eigendafa  simulafing  bofh 
perfecl  measuremenfs  and  model  correlafion,  fhe  APE  algorifhm  converged  fo  fhe  correcf  damaged 
elemenf(s),  and  indicafed  fhe  correcf  percenf  of  damage.  In  each  case,  less  fhan  fen  percenf  of  fhe 
fofal  degrees-of-freedom  were  included  in  fhe  measuremenf  sef,  and  only  a  small  number  of  fhe  nafural 
frequencies.  An  addifional  affribule  of  fhe  algorifhm  is  fhaf,  in  fhe  case  where  fhe  enfire  eigenvector 
can  be  measured  exacfly,  fhe  exacf  solufion  is  obfained  wifhouf  iferafing. 


6.4  An  Analytical  Example 

The  APE  mefhod  was  applied  to  fhe  41-elemenl  free-free  planar  fruss  shown  in  Eigure  6.3. 
Sfrucfural  damage  was  arbifrarily  defined  as  a  50%  reducfion  in  sfiffness  fo  elemenf  #7  and  a  30% 
reducfion  fo  sfiffness  in  elemenf  #18.  An  eigenanalysis  of  fhe  EEM  model  using  fhe  reduced  sfiffness 
values  produced  fhe  (A,,  measured  dafa.  The  frequencies  and  parfial  shapes  corresponding  to 


6-12 


Table  6.1  Damage  identification  results  on  the  41-element  free-free  planar  truss. 


True 

APE 

ASTROS-ID 

element  # 

%  damage 

element  # 

%  damage 

element  # 

%  damage 

7 

50 

7 

50.0 

7 

34.3 

18 

30 

18 

30.0 

10 

34.4 

18 

19.1 

19 

12.6 

cpu  time 

14.6  (sec.) 

540.2  (sec.) 

the  same  instrumented  degrees-of-freedom  (2,  10,  18,  26)  as  used  in  Chapter  IV  for  the  first  five 
flexible  modes  were  used.  Figures  6.4  and  6.5  show  the  resulting  change  in  modal  frequency  and 
shape  from  structural  damage.  The  damaged  configuration  is  plotted  using  dashed  lines,  while  the 
nominal  configuration  is  plotted  using  solid  lines.  As  can  be  seen  from  the  figures,  the  change  in  the 
shapes  due  to  the  damage  is  minimal,  justifying  the  use  of  the  nominal  eigenvectors  for  the  unmeasured 
degrees-of-freedom  to  initiate  the  APE  algorithm.  The  41  fictitious  actuators  corresponding  to  the  41 
structural  elements  were  used  to  construct  B  according  to  Equation  (6.10). 

The  results  of  the  APE  method  are  listed  in  Table  6.1.  Nineteen  iterations  were  used  in  the 
solution  process,  starting  from  an  initial  search  space  of  41  elements.  An  analysis  using  ASTROS-ID 
with  41  design  variables  corresponding  to  the  stiffness  of  the  41  elements  was  also  performed.  The 
results  are  listed  in  Table  6.1,  along  with  the  required  CPU  times  for  each  method.  Pour  outer-loop 
iterations  were  used  in  obtaining  the  ASTROS-ID  results.  Equal  weighting  on  the  eigenvalue  and 
eigenvector  components  were  used  in  the  objective  function,  since  all  quantities  are  known  exactly. 
Note  that  although  both  methods  were  run  on  the  same  CPU  platform  (Sparc- 10  workstation),  the 
methods  were  implemented  in  two  different  software  environments.  Because  of  this,  there  may  be 
some  small  percentage  of  the  CPU  times  reported  that  is  associated  with  processing  unique  to  the 
software  environment  and  not  the  solution  computation.  Table  6. 1  does  illustrate  the  difference  in  the 
nature  of  the  solutions.  Using  the  APE  method,  the  solution  space  is  reduced  to  the  fewest  number  of 
elements  which  can  achieve  the  measured  eigenstructure.  Using  ASTROS-ID,  the  damage  is  equally 
assigned  between  elements  which  have  the  same  eigenvalue  and  eigenvector  sensitivities.  Elements 
number  7  and  number  10  have  exactly  the  same  eigenvalue  sensitivity  and  nearly  the  same  eigenvector 
sensitivity  to  the  measured  data.  The  same  is  true  for  elements  18  and  19.  In  the  case  where  uncertainty 
exists  in  the  measured  data,  it  becomes  increasingly  difficult,  if  not  impossible  to  distinguish  between 
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which  element  to  correctly  assign  the  damaged  value.  A  slight  perturbation  in  the  measured  data  can 
result  in  the  APE  method  assigning  a  50%  damage  to  element  10  rather  than  element  7.  The  problem 
of  distinguishing  between  elements  with  similar  sensitivities  is  termed  ‘damage  localization’  and  is 
discussed  in  Chapter  VII. 

6.5  Summary 

A  method  was  presented  to  achieve  a  set  of  measured  eigenvalues  and  eigenvectors  through 
changes  in  the  stiffness  matrix  which  are  consistent  with  the  finite  element  formulation.  The  method 
does  not  require  the  use  of  full-length  eigenvectors  in  its  formulation,  and  is  suitable  for  use  with 
minimally  instrumented  structures.  The  method  was  demonstrated  on  an  analytical  model,  accurately 
identifying  the  damaged  elements.  A  comparison  between  ASTROS-ID  and  APE  was  also  performed. 
Common  to  both  methods  is  the  question  of  which  degrees-of-freedom  to  instrument  to  facilitate 
damage  identification,  and  to  what  extent  can  the  damage  be  localized  from  the  measured  data.  These 
two  questions  are  addressed  in  Chapter  VII. 
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VII.  Sensor  Prioritization  and  Damage  Localization 

7.1  Overview 

For  large  space  structures,  practicality  dictates  that  only  partial  modal  data  can  be  measured.  With 
minimal  sensor  information,  two  questions  naturally  arise:  at  which  locations  should  the  sensors  be 
placed,  and  to  what  extent  can  damage  be  identified  with  the  selected  sensor  locations?  An  eigenstruc- 
ture  sensitivity  based  method  is  presented  to  answer  these  questions.  While  sensor  placement  methods 
have  been  presented  by  Lim  ,  Kammer  ,  and  Liu  and  Tasker  ,  they  focused  on  maximizing 
either  controllability  or  observability  and  not  damage  detection.  Kashangaki  introduced  a  modal 
sensitivity  parameter  as  a  quantitative  measure  of  the  eigenvalue  and  eigenvector  sensitivity,  and  used 
this  to  determine  which  modes  should  be  used  in  a  damage  detection  scheme.  In  practice  however,  for 
a  given  complex  structure,  only  a  few  of  the  lower  frequency  global  modes  can  be  accurately  identified. 
Af  higher  frequencies,  fhe  separafion  of  local  and  global  modes  becomes  increasingly  difficull,  if  nol 
impossible.  Furfhermore,  only  a  few  degrees-of-freedom  can  be  insfrumenfed.  Therefore,  emphasis 
herein  is  placed  on  priorifizing  sensor  locafions  and  on  fhe  ability  fo  localize  damage  from  parfial 
eigendafa  for  a  given  number  of  modes,  and  nol  on  which  modes  lo  measure. 

The  melhod  presenled  is  based  on  examining  fhe  firsl-order  parfial  eigenslruclure  sensilivily  lo 
changes  in  fhe  slruclural  sliffness  of  each  elemenl  of  a  finite  elemenl  model.  No  a  priori  knowledge 
of  fhe  damage  localion  is  assumed.  Two  aspecls  of  fhe  parfial  eigenslruclure  sensilivily  are  explored. 
Firsl,  is  Ihe  amounl  by  which  varialions  of  Ihe  elemenlal  sliffness  values  change  Ihe  measured  parfial 
eigenslruclure.  Independenl  of  Ihe  damage  deleclion  scheme  used,  elemenls  which  produce  little  or 
no  change  in  Ihe  measured  dala  will  be  difficull  or  impossible  lo  deled  when  damaged.  Second,  is  Ihe 
direclion  of  change  in  Ihe  parfial  eigenslruclure.  Elemenls  which  produce  similar  or  idenlical  changes 
in  Ihe  parfial  eigenslruclure,  will  be  difficull  or  impossible  lo  dislinguish  belween  when  damaged. 
Therefore,  sensor  localions  are  chosen  so  lhal  Ihe  change  in  Ihe  measured  parfial  eigenslruclure  due  lo 
damage  is  maximized.  Localizalion  of  Ihe  damage  lo  an  elemenl(s)  is  based  on  bolh  Ihe  amounl  and 
direclion  of  change  lo  Ihe  parfial  eigendala  for  Ihe  chosen  sensor  localions. 
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7.2  Sensitivity  Theory 


Eigenvalue  and  eigenvector  sensitivity  to  changes  in  structural  elements  will  be  based  upon  the 
finite  element  model  of  the  structure.  For  on-orbit  damage  scenarios  of  large  flexible  space  structures, 
two  assumptions  are  made.  First,  structural  damage  is  confined  to  changes  in  the  stiffness  properties 
of  the  structure.  Second,  structural  damping  is  negligible.  With  these  assumptions,  the  free  vibration 
of  the  structure  is  modeled  as: 


Mx  -\-  (K  —  AK)x  =  0  (6.1) 

with  all  variables  as  previously  defined.  The  eigenvalue  and  eigenvector  for  the  mode  of  Equation 
(6.1)  is  given  as  (A,,  <i>j)  whereas  the  measured  eigenvalue  and  partial  eigenvector  for  the  same  mode 
is  represented  as  The  relationship  between  the  n  dimensional  eigenvectors  <i>j  and  the  s 

dimensional  partial  eigenvectors  (fn  is  (fn  =  The  matrix  C  G  maps  the  full  length 

eigenvector  into  the  partial  eigenvector  corresponding  to  the  measured  degrees-of-freedom.  With 
minimal  sensor  information  available,  a  natural  cost  function  representing  the  mismatch  between  the 
eigenstructure  of  the  finite  element  model  and  the  measured  eigendata  is: 

^  ( W  -  l)  +  ^  ^  h  (<^6  -  ^ij)  (3-2) 

Of  interest  for  sensor  location  determination  is  how  to  choose  the  matrix  C  such  that  structural  damage 
results  in  observable  changes  in  A,  and  cj)-  and  hence  in  J.  Once  C  is  determined,  damage  localization  is 
concerned  with  the  uniqueness  of  changes  in  A,  and  cj)-  for  variations  in  the  matrix  AK.  As  developed 
in  Chapter  VI,  the  structural  constraint  can  be  imposed  on  AK  by  expressing  it  as: 

AK  =  BGB'^  (6.4) 

where  B  is  constructed  from  the  nodal  connectivity  information  and  the  elemental  parameters.  G  is 
a  diagonal  matrix  (6?  =  diag(^i  . .  .gp))  composed  of  the  fraction  of  damage  for  each  element  (gi). 
A  value  of  gi  =  0  corresponds  to  an  undamaged  element  whereas  gi  =  1  corresponds  to  a  complete 
loss  of  stiffness  to  the  element.  The  eigenvalue  and  eigenvector  sensitivity  to  structural  damage 
is  computed  based  on  the  method  presented  by  Fox  and  Kapoor. The  sensitivity  calculations  are 
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consistent  with  the  method  presented  in  Chapter  V,  however  with  damage  being  confined  to  the  stiffness 
matrix  the  calculations  are  further  simplified.  Furfhermore,  fhe  mafrix  operafor  P(a,  (3)  as  defined  in 
Chapfer  IV  is  used  fo  develop  fhe  sensifivify  equafions.  The  eigenvalue  equafion  is  wriffen  as: 


(-XiM  +  K  -  BGB'^)^i  =  0 


(7.1) 


Wifh  fhe  assumpfion  fhaf  changes  from  sfrucfural  damage  are  confined  fo  fhe  perfurhafion  mafrix  Alf , 
Equafion  (7.1)  is  differenlialed,  which  after  simplifying  resulfs  in: 

M^,^  +  ^{BGB'^^i)  =  {)  (7.2) 

dg  dg 

As  developed  in  Chapfer  VI,  fhe  mafrix  operafor  P(q;,  /3)  is  defined  as: 

n 

P(a,(3)  with  P,  (3  E  G  ,  P^j  =  ^2'^kl3ijl3kj  (6.23) 

k  =  l 

where  is  fhe  row  and  fhe  column  of  fhe  mafrix  P.  In  terms  of  fhe  operator  P,  fhe  mafrix 
producf  in  parenfhesis  in  Equafion  (7.2)  can  he  wriffen  as: 

BGB^^i  =  P{^i,B)g  where  ^  =  diag(G')  ,  6?  G  3?^^^(diagonal)  (7.3) 


Each  eigenvector  is  normalized  so  fhaf  <i>f  M^i  =  1.  Premulfiplying  Equafion  (7.2)  hy  <i>f  and  using 
fhe  operafor  P,  fhe  eigenvalue  sensifivify  from  Equafion  (7.2)  can  he  wriffen  as: 

f)\  f)\ 

^  =  -ci>fP(<i>,,P)  ,  ^g5?ix^  (7.4) 

dg  dg 

In  a  similar  fashion,  eigenvector  sensifivify  is  computed  hy  differenfiafing  Equafion  (7.1)  and  using  fhe 
resulfs  of  Equafion  (7.4).  The  eigenvector  sensifivify  for  fhe  mode  is: 

[K  -  \M  -  BGB^]  ^  +  /]  P($„  P)  ,  ^  G  3?"^"  (7.5) 

dg  dg 

The  mafrix  /  denofes  fhe  n  X  n  idenfify  mafrix.  The  mefhod  infroduced  hy  Nelson  is  used  to  solve 
Equafion  (7.5).  This  is  necessary  due  fo  fhe  singularifyoffhe  mafrix  [If  —  AjM  —  BGB^].  Assuming 
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no  repeated  roots,  this  method  involves  separating  the  solution  into  a  particular  and  homogeneous 
solution,  where: 


dg 


+  Vi 


(7.6) 


with: 

V;-  =  [k  -  \M  -  BGB^]  [Mi, if  +  /]  P($„  B)  ill) 


and 

(7.8) 

The  (•)  notation  represents  matrices  reduced  hy  one  row  and  one  column,  or  vectors  reduced  hy  one 
element.  Nelson’s  method  removes  the  row  and  column  corresponding  to  the  maximum  entry  in  <!>,  . 
Equation  (7.8)  corresponds  to  the  maximum  entry  occurring  at  the  element.  For  computational 
efficiency,  a  decomposition  and  substitution  is  preferable  to  explicitly  computing  the  inverse.  The 
constant  row  vector  c,  is  given  as: 

c,  =  -^jMV,  (7.9) 


With  the  first-order  eigenvalue  and  eigenvector  sensitivities  defined,  lef  V  A  be  fhe  mafrix  whose 
row  and  column  enfry  is  defined  as: 


_  OA,  1_ 


(7.10) 


The  term  A,  is  infroduced  fo  correcf  fhe  scaling  of  fhe  differenl  modes.  Wifh  fhis  definilion,  each 
column  of  VA  corresponds  fo  differenl  slruclural  elemenls  and  each  row  lo  a  differenl  mode.  The 
eigenvalue  change  from  changes  in  fhe  slruclural  elemenls  Ag  G  lo  lirsl  order,  is  given  as: 


AA  =  VXAg 


(7.11) 


where 

AA  =  [AAi, . . . ,  AAr]^  and  AXi  =  V— — A  ,Vl  =  l,...,r  (7.12) 

^oi 

The  vector  A  A  G  5?’’^  ^  consisls  of  Ihe  fraclional  changes  to  Ihe  r  measured  eigenvalues  of  Ihe  slruclure 
due  to  damage.  Similarly,  Ihe  vector  AA^  G  conlains  Ihe  eigenvalues  of  Ihe  slruclure  evaluated 
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at  BGgB^ .  For  the  undamaged  structure,  BGgB^  =  0.  Since  Equation  (7.11)  is  valid  only  for 
small  changes  of  Ag,  it  is  not  possible  to  use  it  directly  to  determine  damaged  elements.  However, 
for  the  purpose  used  herein.  Equation  (7.11)  is  adequate  to  examine  the  relationship  between  A  A  and 
Ag.  Information  on  the  amount  and  direction  of  changes  in  A  A  are  contained  in  the  matrix  VA.  With 
Agi  G  [0,  l],rows  of  VA  with  negligible  norms  will  contribute  negligibly  to  changes  in  AA.  Rows  in 
VA  which  are  similar  or  identical  to  one  another  will  have  values  of  Agi  that  affect  A  A  similarly  or 
identically,  and  hence  will  be  indistinguishable  from  one  another.  Eor  this  analysis,  it  is  assumed  that 
11^112  is  small  and  that,  although  individual  elements  in  A^  may  be  close  to  unity,  the  overall  effect  on 
the  global  nature  of  the  structure  is  small,  i.e.  no  catastrophic  failures. 

Using  the  results  from  Equation  (7.5),  the  partial  eigenvector  sensitivity  for  the  mode  is 
defined  as: 

V4>,  =  ,  V4>,  G  (7.13) 

dg 

where  C  is  determined  from  the  measured  degrees-of-freedom  as  previously  defined.  Wifh  fhis 
definifion,  changes  in  fhe  parfial  eigenvector  for  fhe  mode,  to  firsf  order  in  Ag,  is  given  as: 

A(i)i  =  V(f)iAg  where  A(f)i  =  (/){  -  (f)i^  (7.14) 

The  vector  (fn  G  is  fhe  parfial  eigenvector  for  fhe  insfrumenfed  degrees-of-freedom,  and  G 
confains  fhe  parfial  eigenvector  of  fhe  finite  elemenf  model  evaluated  af  BGgB^ .  Similar  to 
fhe  eigenvalue  case,  informalion  on  fhe  amounf  and  direcfion  of  change  of  A(j)i  is  confained  in  fhe 
mafrix  V (fn .  Nofe  fhaf  fhere  is  one  mafrix  V (fn  for  each  measured  mode.  Informalion  on  which  to 
base  bolh  sensor  localion  and  damage  localizalion  is  confained  in  fhe  malrices  VA  and  V(^j.  Two 
properlies  of  Ihese  malrices  are  invesligaled,  which  are  referred  to  as  fhe  deleclabilily  and  colinearily. 
Deleclabilily  is  a  measure  of  fhe  amounf  of  change  which  occurs  from  changes  in  a  design  variable, 
whereas  colinearily  is  a  measure  of  fhe  direcfion  of  change. 

7.2.1  Sensor  Location  Prioritization.  Inilially  wifh  C  =  I,  y4’i  from  Equalion  (7.13) 
confains  informalion  indicaling  which  degrees-of-freedom  to  inslrumenl.  As  previously  discussed,  on 
orbil  only  a  few  of  fhe  low  frequency  modes  can  be  measured.  Given  Ibis  fad  and  fhe  problem  of 
solulion  non-uniqueness  associated  wifh  using  parfial  measuremenl  dala,  if  is  assumed  lhal  any  mode 
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which  can  he  measured  should  he  used  in  a  damage  detection  scheme.  Given  then  that  there  are  r 
measured  modes,  the  detectability  in  the  measured  eigenvectors  at  the  degree-of-freedom  from 
changes  in  the  p  elements  of  the  structure  is  defined  as: 

p  r 

=  (7.15) 

k=l i=l 

The  vector  D^=  is  then  sorted  in  descending  order,  initially  prioritizing  the  sensor 

locations  based  on  detectability.  A  threshold  is  set  based  on  the  measurement  uncertainty  and  the  finite 
element  modeling  errors.  Values  of  below  this  threshold  indicate  degrees-of-freedom  which  are 
unaffected  by  structural  damage  for  the  measured  modes.  Next,  a  colinearity  check  is  made  to  determine 
degrees-of-freedom  which  yield  similar  information  on  the  damaged  elements.  The  colinearity,  denoted 
S,  between  any  two  vectors  a  and  /3  is  defined  as: 

=  0^/3  ||q;||2  =  II/3II2  =  1  (7.16) 


which  is  simply  fhe  cosine  of  fhe  angle  befween  fhe  fwo  vectors.  A  value  of  =  1  indicate  perfecl 
colinearify  whereas  =  0  indicates  orfhogonal  vectors.  Wifh  fhis  definilion,  fhe  colinearify  of  fhe 
eigenvector  sensifivify  befween  measured  degrees-of-freedom  /  and  m  is  defined  as: 


5.,^  = 


-  y:  [V®, .  vof] 


i  =  l 


Im 


(7.17) 


Again  a  fhreshold  from  unify  is  sef  based  on  fhe  measuremenf  uncerfainfy  and  fhe  finite  elemenf 
modeling  errors.  Values  in  5^G  wifhin  fhis  fhreshold  are  declared  colinear,  indicafing  fhaf  af 

fhese  degrees-of-freedom  for  fhe  r  measured  modes,  fhe  changes  in  fhe  eigenvector  are  indisfinguishable 
from  one  anofher  to  changes  in  fhe  sfrucfural  elemenfs.  Using  fhis  informalion,  mulfiple  colinear  enfries 
in  vector  are  removed,  leaving  only  one  enfry  from  each  colinear  grouping.  The  remaining  firsf 
s  elemenfs  of  vector  represenf  fhe  priorifized  s  degrees-of-freedom  to  place  sensors.  Wifh 
defined  as  fhe  firsf  s  elemenfs  of  fhe  reduced  and  sorted  vector  D^,  fhe  mafrix  C  is  chosen  such  fhaf 
=  CD^  is  safisfied.  An  analyfical  example  is  presented  in  a  subsequenf  secfion,  following  fhe 
discussion  of  damage  localizafion. 
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7.2.2  Damage  Localization.  Given  the  r  modes  measured  at  the  s  degrees-of-freedom 
as  determined  above,  damage  localization  determines  the  extent  to  which  damage  can  he  isolated  to 
individual  elements.  Similar  detectahility  and  colinearity  metrics  are  used,  which  are  now  restricted  to 
the  instrumented  degrees-of-freedom.  The  detectahility  in  the  measured  eigenvalues  from  changes  in 
the  structural  element  is  defined  as: 


r 

=  E|VA,,|  (7.18) 

i  =  l 

where  Dx=  [Dxi  ■  ■  ■  ■  Eigenvalue  colinearity  information  is  contained  in  Sx-  The  row  and 

column  of  Sx  indicates  the  colinearity  of  the  eigenvalue  changes  between  the  and  structural 
elements  and  is  defined  as: 

=  [VA^-VA]^.,  (7.19) 

The  eigenvalue  colinearify  is  independenf  of  fhe  degree-of-freedom  af  which  if  is  measured.  Similarly, 
defecfabilily  in  fhe  measured  eigenvectors  from  fhe  sfrucfural  elemenf  is  defined  as: 

5  r 

=  (7-20) 

1=1  8  =  1 

where  Eigenvector  colinearify  informafion  is  confained  in  S^.  The  row 

and  A;*^  column  of  indicates  fhe  colinearify  of  fhe  eigenvector  changes  befween  fhe  and  A;*^ 
sfrucfural  elemenfs  and  is  defined  as: 


= 


-  y:  i^of  ■  V®,] 


i  =  l 


■i  jk 


(7.21) 


Note  fhe  similarity  befween  Equafions  (7.15)  and  (7.17)  and  Equafions  (7.20)  and  (7.21).  Eor  fhis 
reason,  sensor  priorifizafion  and  damage  localizafion  are  considered  dual  problems,  eifher  of  which  can 
be  defermined  wifh  only  a  slighf  modificalion  to  fhe  same  algorifhm.  Nofafionally,  defeclabilily  D  and 
colinearify  S  were  mulfiply  defined,  once  for  sensor  priorifizafion  and  again  for  damage  defecfion.  If 
should  be  clear  from  fhe  confexf  of  fhe  problem  which  definilion  applies. 

Wifh  fhese  definilions,  damage  localizafion  proceeds  as  follows.  Elemenfs  in  Dx  and 
which  are  below  fhe  modeling  and  measuremenf  uncerfainly  fhreshold  level  are  declared  undefecfable. 
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Damage  in  an  undetectable  element  cannot  be  identified  from  the  measured  data.  Of  the  remaining 
elements,  colinear  elements,  as  indicated  by  elements  in  Sx  and  above  the  uncertainty  level  are 
indistinguishable  from  one  another.  From  the  measured  data,  structural  damage  can  only  be  localized 
to  a  colinear  group,  and  not  to  an  individual  element  within  the  group.  Elements  contained  in  a  colinear 
group  are  referred  to  as  symmetric  elements. 

7.3  Software  Implementation 

The  sensor  prioritization  and  damage  localization  method  was  implemented  using  MATLAB® 
software.  For  a  given  number  of  modes,  the  eigenvalue  and  eigenvector  sensitivities  are  computed 
using  Equations  (7.4)  and  (7.6)  respectively.  From  these,  the  detectability  metrics  Dx  and  are 
computed  using  Equations  (7.18)  and  (7.20).  The  colinearity  metrics  Sx  and  are  computed  using 
either  Equations  (7.15)  and  (7.17)  for  sensor  prioritization  or  using  the  transposes  as  given  in  Equations 
(7.19)  and  (7.21)  for  damage  localization.  These  metrics  are  then  compared  against  threshold  values 
based  on  model  and  measurement  uncertainty.  The  uncertainty  is  determined  by  how  well  the  finite 
element  model  correlates  to  the  measured  data,  for  nominal  as  well  as  damage  configurations.  The 
detectability  threshold  was  established  as  a  percentage  of  the  maximum  element  in  the  vector  D.  For 
colinearity  the  threshold  was  a  percentage  decrease  from  unity  value.  With  the  thresholds  established, 
the  elements  of  the  structure  are  then  classihed  as  either  undetectable  (U)  using  Dx  and  D^,  symmetric 
(S),  or  identihable  (I).  For  computational  efficiency,  detectability  is  checked  hrst.  Any  elements 
with  values  below  the  detectability  threshold  are  classihed  as  U  and  are  removed  from  the  sensitivity 
gradient  matrices  before  forming  the  colinearity  metrics.  Colinearity  groupings  are  then  determined 
from  the  colinearity  matrix  by  replacing  the  entries  of  the  matrix  with  either  ones  or  zeros  based  on 
being  above  or  below  the  colinearity  threshold.  In  this  way,  a  nonzero  entry  in  the  row  and 
column  indicates  symmetry  between  the  and  elements.  Note  that  only  the  entries  above  the  main 
diagonal  need  be  computed.  Based  on  these  entries,  the  elements  are  classihed  as  either  S  or  I.  For 
a  colinear  grouping  of  elements  as  determined  by  Sx  and  S^,  one  element  is  classihed  as  identihable 
and  the  remaining  as  symmetric.  The  selection  of  the  weighting  between  emphasis  on  Sx  and  on 
S^  is  dependent  on  the  damage  identihcation  scheme  used.  For  damage  identihcation  based  on  the 
cost  function  minimization,  such  as  given  in  Equation  (3.2),  the  metric  results  should  be  combined 
consistently  with  the  cost  function  weighting  coefficients.  For  example,  a  high  relative  weighting  on 
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the  eigenvalues  corresponds  to  an  increased  emphasis  on  the  Sx  metric.  A  decision  flow  chart  for  the 
damage  localization  process  as  implemented  for  use  with  the  APE  method  is  shown  in  Figure  7.1.  The 
decision  chart  reflects  the  emphasis  the  APE  method  places  on  eigenvalues  over  the  eigenvectors. 

Independent  of  the  algorithm  used  for  damage  identification,  the  advantage  of  the  U/S/I  classifi¬ 
cation  is  apparent.  It  quickly  indicates  which  damaged  elements  cannot  he  detected  from  the  measured 
data.  Furthermore,  only  elements  in  the  I  classification  need  he  included  in  the  search  space.  It  has 
been  observed  that  decreasing  the  search  space  precludes  singularities  of  the  N  matrix  in  Equation 
(6.34)  and  significantly  reduces  the  required  solution  time. 

7.4  An  Analytical  Example 

To  demonstrate  the  use  of  the  detectability  and  colinearity  metrics,  an  analytical  example  of  both 
the  sensor  prioritization  and  damage  localization  was  performed  using  the  same  41 -element  free-free 
planar  truss  assembly  as  shown  in  Figure  7.2  and  described  previously  in  Chapter  IV. 

The  first  analysis  was  an  examination  of  the  relationship  between  increasing  the  number  of 
measured  modes  and  increasing  the  number  of  sensors.  Table  7.1  contains  the  tabulated  results.  The 
data  clearly  shows  that,  if  possible,  increasing  the  number  of  modes  measured  is  preferable  to  increasing 
the  number  of  sensors  to  enhance  damage  localization.  For  all  cases,  the  threshold  values  were  fixed  af 
10%  for  defecfabilily  and  5%  colinearify.  If  should  be  noted  fhaf  fhe  fabulaled  resulfs  in  some  insfances 
show  fhaf  adding  addifional  sensors  had  an  adverse  effecl  on  idenfiflabilify.  This  frend  is  an  arlifacl  of 
using  differenl  lengfh  vectors,  due  to  a  differenl  number  of  sensors,  compared  againsf  a  fixed  fhreshold. 
The  dafa  as  presenfed  is  intended  only  to  show  fhe  overall  frends. 

To  demonsfrafe  sensor  priorifizafion,  fhe  locafions  of  fhree  sensors  were  selected  using  a  fixed 
number  of  modes  and  fhe  same  fhresholds  as  sfafed  above.  As  shown  in  Table  7.2,  fhe  priorifized 
locafions  increase  fhe  number  of  idenfiflable  elemenfs  over  fwo  randomly  chosen  sensor  locafions. 

A  fhird  analysis  was  performed  to  demonsfrafe  fhe  validify  of  using  only  flrsf-order  sensifivifies, 
evaluated  af  fhe  nominal  conflgurafion,  over  fhe  range  of  A^.  Table  7.3  lisfs  fhe  symmefric  elemenfs  as 
defermined  by  Sx  for  fhe  flrsf  five  flexible  modes,  for  fhe  same  fhreshold  values  as  used  above  and  using 
a  single  sensor.  As  an  example,  fhe  resulfs  indicate  elemenfs  8, 9, 33,  and  34  are  symmefric  to  fhe  parfial 
measured  dafa,  and  fhus  only  one  elemenf  should  be  included  in  fhe  damage  defecfion  search  space. 
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Table  7.1  Damage  localization  results  for  a  41-element  free-free  planar  truss. 


Measured  Flexible  Modes 

#  Sensors 

1 

1,2 

1,2,3 

1, 2,3,4 

1,2,3,4,5 

1 

29/9/3 t 

21/15/5 

21/15/5 

9/25/7 

9/24/8 

2 

29/9/3 

21/8/12 

21/3/17 

9/6/26 

9/1/31 

3 

29/2/10 

21/2/18 

21/0/20 

9/3/29 

9/1/31 

4 

29/0/12 

21/1/19 

21/0/20 

9/4/28 

9/0/32 

5 

29/0/12 

21/3/17 

21/0/20 

9/3/29 

9/0/32 

6 

29/0/12 

21/3/17 

21/0/20 

9/2/30 

9/0/32 

Locations + 

34,2,1,3,33,35 

1,35,3,33,9,27 

1,35,33,3,36,2 

1,35,33,3,36,2 

1,35,3,33,16,22 

fData  presented  in  U/S/I  format  where  U  denotes  the  number  of  undetectable  elements,  S  the  number  of  symmetric  elements, 
and  I  the  number  of  identifiable  elements.  The  sum  of  U,  S,  and  I  equals  the  total  number  of  elements  of  the  structure, 
f  The  locations  of  the  sensors  were  chosen  using  the  prioritization  method  presented  and  are  reported  by  degree-of-freedom 
number  in  prioritized  order. 


Table  7.2  Damage  localization  for  different  sensor  locations  for  a  41-element  free-free  planar  truss. 


Method 

Sensor  Location^ 

Damage  Localization  (U/S/I)+ 

Prioritized 

1,35,3 

9/1/31 

Random 

3,11,19 

9/5/27 

Random 

2,4,34 

9/3/29 

fReported  by  degree-of-freedom  number. 
fBased  on  measuring  the  first  five  flexible  modes. 


7-11 


Table  7.3  Damage  localization  based  on  eigenvalue  sensitivity  using  the  first  5  flexible  modes,  for 
the  41-element  free-free  planar  truss. 


Element  # 

Equivalent  Symmetric  Elements  (5*^)^ 

2 

5,37,40 

3 

4,38,39 

7 

10,32,35 

8 

9,33,34 

12 

15,27,30 

13 

14,28,29 

17 

20,22,25 

18 

19,23,24 

Undetectable  Elements 

1,6,11,16,21,26,31,36,41 

fResults  are  independent  of  selected  sensor  location. 

The  measured  eigenvalues  with  damage  to  element  8,  gs  G  [0, 1]  corresponds  to  the  same  eigenvalues 
for  an  equivalent  break  in  either  elements  9,  33,  or  34.  For  an  example  with  multiple  breaks,  elements  8 
and  30  were  reduced  by  70%  and  50%  respectively.  Using  any  combination  of  two  elements,  selecting 
one  element  from  the  symmetric  set  (8,9,33,34),  and  one  element  from  the  symmetric  set  (12,15,27,30), 
reduced  by  70%  and  50%  respectively,  the  changes  in  the  measured  eigenvalues  were  within  2%  of 
each  other.  This  is  well  within  the  uncertainty  level  established  in  the  threshold  for  Sx  used  to  select 
the  symmetric  sets. 

7.5  Summary 

A  method  was  presented  which  prioritizes  the  degrees-of-freedom  to  instrument  when  used  to 
collect  modal  data  for  damage  detection.  It  was  shown  that  this  method  can  also  be  used  to  determine 
the  extent  to  which  damage  can  be  localized  from  these  sensor  locations.  The  method  represents  a 
computationally  attractive  alternative  to  an  exhaustive  search  over  the  parameter  space  and  is  a  valuable 
tool  during  the  design  phase  to  determine  measurement  and/or  modeling  accuracy  requirements.  An 
analytical  example  was  presented  which  showed  that  the  extent  to  which  damage  can  be  localized  is 
limited  by  the  amount  and  quality  of  the  measured  data.  It  was  also  shown  that  increasing  the  number 
of  measured  modes  is  of  greater  benefit  than  increasing  the  number  of  sensors. 

Having  developed  the  four  tasks  associated  with  damage  identification,  each  task  was  combined 
into  an  integrated  software  package  and  programmed  in  MATLAB®.  A  description  of  this  software 
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tool  is  provided  in  Appendix  B.  Using  this  software,  the  four  tasks  of  damage  identification  were 
applied  to  experimental  structures,  the  details  of  which  are  discussed  in  the  Chapter  VIII. 
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VIII.  Experimental  Validation  of  Theory 


8.1  Test  Objective 

The  test  objective  is  the  experimental  demonstration  of  the  ability  to  identify  structural  damage 
from  a  simulated  damaged  space  truss  using  limited  measurement  data.  In  choosing  an  experimental 
apparatus  to  validate  a  new  technique,  it  is  important  for  the  experiment  to  exhibit  the  common  dynamic 
characteristics  of  the  intended  application.  For  large  flexible  space  structures,  these  characteristics 
include  a  low  fundamental  frequency  (on  the  order  of  1  Hz)  and  high  modal  density  at  low  frequencies 
with  low  modal  damping  ratios,  and  a  truss-like  structure.  Furthermore,  it  is  important  that  both 
the  excitation  actuator  and  measurement  sensors  be  non-grounded  since  only  these  type  devices  are 
applicable  to  space-based  applications.  Lastly,  to  demonstrate  the  detection  of  actual  damage,  it 
is  desirable  to  have  a  structure  which  contains  structural  elements  which  can  easily  be  modified  or 
removed  to  simulate  a  failure.  Two  different  test  beds  where  chosen  to  validate  the  theory.  The  first 
validation  was  through  the  use  of  experimental  test  data  obtained  from  NASA’s  8-bay  truss  test  bed, 
and  the  second  using  AFIT’s  six-meter  Flexible  Truss  Experiment. 

8.2  NASA  Test  Data  Analysis 

8.2.1  Hardware.  NASA’s  8-bay  truss  test  bed  consisted  of  eight  cubic  bays  of  a  hybrid 
space  truss  cantilevered  from  a  rigid  backstop  plate.  This  configuration  represents  a  scaled  section  of 
the  proposed  International  Space  Station.  Each  bay  is  a  half  meter  in  length,  constructed  of  aluminum 
members.  A  typical  joint  configuration  is  shown  in  Eigure  8.1.  The  truss  was  fully  instrumented  with 
one  triaxial  accelerometer  at  each  of  its  32  unconstrained  nodes.  Disturbance  excitation  was  achieved 
using  two  ground-based  dynamic  shakers  attached  at  two  different  node  points.  A  complete  description 
of  the  hardware  and  the  testing  procedure  is  contained  in  the  work  by  Kashangaki.  ™ 

8.2.2  Model  Tuning.  The  NASA  truss  was  modeled  using  104  rod  elements,  with  three 
translational  degrees-of-freedom  per  node.  The  material  properties  of  all  elements  were  identical. 
Eumped  masses  were  incorporated  into  the  model  to  account  for  the  mass  of  the  node  balls,  standoffs, 
sleeves,  collars,  and  instrumentation.  The  data  received  from  NASA  (T.  Kashangaki)  consisted  of  the 
identified  nafural  frequencies  and  shapes  for  fhe  flrsf  five  flexible  modes  of  fhe  fruss,  for  fhe  nominal 
as  well  as  several  damage  conflgurafions.  The  flrsf  five  modes  in  numerical  order  consisted  of  fhe 
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Figure  8.1  Joint  construction  for  the  NASA  8-bay  truss  experiment. 


first  bending  modes  about  the  ‘X’  and  ‘Z’  axes,  the  first  torsional  mode  about  the  ‘Y’  axis,  and  the 
second  bending  modes  about  the  ‘X’  and  ‘Z’  axes.  A  typical  frequency  response  function  is  depicted 
in  Figure  8.2,  measured  at  the  free  end  of  the  truss.  Both  the  nominal  and  damaged  measurements  are 
shown.  As  can  be  seen  in  the  figure,  identification  of  global  modes  past  80  Hertz  is  difficult  due  to  the 
presence  of  local  element  modes  in  this  regime.  Hence  only  the  first  five  modes  could  be  accurately 
measured.  Also  depicted  in  the  frequency  response  functions  is  the  effect  of  structural  damage.  For 
this  particular  damage  case  which  corresponds  to  the  removal  of  a  longeron  element  at  the  cantilevered 
end,  damage  is  manifest  in  the  frequency  response  by  a  separation  of  the  first  two  flexible  modes  only. 

Nine  damage  cases  were  tested.  The  damage  configurations  were  the  full  removal  of  one  or  two 
elements  of  the  truss.  The  different  damage  configurations  are  shown  in  Figures  8.3,  8.4,  and  8.5.  For 
the  nominal  configuration,  i.e.  no  damage,  the  results  of  the  finite  element  analysis  and  the  measured 
data  are  compared  in  Table  8.1.  Although  the  FEM  model  is  in  fair  agreement  with  the  measured  data, 
any  disagreement  will  result  in  the  damage  identification  method  assigning  a  percentage  of  damage  to 
an  element(s)  to  account  for  the  disagreement.  Therefore,  tuning  was  performed  using  ASTROS-ID 
to  ensure  the  initial  disagreement  was  as  minimal  as  possible.  A  total  of  110  design  variables  were 
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Figure  8.2  Measured  frequency  response  function  for  the  NASA  truss. 


CASE  A 


CASE  B 


CASE  C 


Figure  8.3  Damage  configurations  for  the  NASA  truss,  cases  A,  B,  and  C 
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Table  8.1  Measured  natural  frequencies  of  NASA’s  8-bay  truss. 


Mode  # 

Freqi 

Measured 

lency  (H: 
Inilial 

z) 

Tuned 

%  diff 

Inilial 

erence 

Tuned 

1 

13.88 

13.79 

13.88 

0.65 

0.00 

2 

14.48 

14.31 

14.47 

1.19 

0.07 

3 

48.41 

50.53 

48.40 

4.37 

0.02 

4 

64.03 

65.98 

64.03 

3.05 

0.00 

5 

67.46 

71.20 

67.52 

5.55 

0.08 

used  to  tune  the  model,  which  included  the  elastic  modulus  of  the  104  rod  elements.  Two  design 
variables  were  used  for  the  mass  of  the  rods,  one  for  the  half  meter  length  rods  and  a  second  for  the 
diagonal  members  (1  /v/2  in  length).  The  remaining  four  design  variables  were  used  to  adjust  the  mass 
properties  of  the  lumped  nodal  masses.  There  were  four  different  nodal  configurations,  depending  on 
how  many  members  were  connected  at  a  node  point  and  whether  or  not  an  external  shaker  attachment 
was  included.  To  tune  the  model,  all  five  measured  modes  were  used  wifh  equal  weighfing  on  all 
frequencies.  Eigenvector  informalion  was  nol  included  in  fhe  Inning  process  for  several  reasons. 
Firsf,  fhe  priorifized  sensor  locafions  were  nol  chosen  unlil  after  fhe  model  was  luned,  and  Ihus  which 
elemenls  of  fhe  eigenvector  lo  include  in  fhe  Inning  process  was  yel  unknown.  (Eigenvalue  informalion 
is  independenl  of  fhe  sensor  locafions,  assuming  a  sensor  is  nol  located  al  a  node  poinl  for  a  parlicular 
mode.)  Nole  lhal  allhough  Ihe  full  lenglh  eigenvectors  are  available  from  Ihe  lesl  dala,  Ihis  does  nol 
represenl  a  realislic  on-orbil  capability  and  Ihus  complete  knowledge  of  Ihe  eigenvectors  was  nol  used. 
A  second  reason  for  nol  including  eigenvector  informalion  was  Ihe  desire  to  mainlain  Ihe  symmelry 
of  Ihe  slruclure.  Wilh  110  design  variables,  Ihere  was  adequate  design  freedom  to  nearly  achieve 
any  parlial  measured  eigenvectors.  Thus  any  measuremenl  error  in  Ihe  parlial  vectors  could  alter  Ihe 
symmelry  of  Ihe  model.  The  resull  of  Ihe  luning  is  listed  in  Table  8. 1 .  The  Inning  process  required  23.4 
minutes  of  CPU  lime  for  four  outer-loop  ileralions,  and  resulted  in  a  decrease  in  Ihe  objeclive  funclion 
by  four  orders  of  magnilude.  All  design  variables  remained  wilhin  10%  of  Iheir  nominal  values. 

Using  Ihe  luned  model,  an  eigenanalysis  for  each  of  Ihe  nine  damage  cases  was  performed. 
The  resulls,  along  wilh  Ihe  experimenlally  measured  dala  from  Ihe  damaged  slruclure,  is  presented 
in  Table  8.2  for  comparison.  In  order  for  Ihe  objeclive  funclion  minimizalion  to  be  successful,  il  is 
imporlanl  lhal  Ihe  analylical  model  wilh  simulated  damage  correlate  well  wilh  Ihe  measured  dala  of  Ihe 
damaged  slruclure.  Unforlunalely,  luning  to  Ihe  nominal  as  well  as  Ihe  damaged  dala  is  only  possible 
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Table  8.2  Changes  in  natural  frequencies  from  damage  on  the  NASA  truss. 


Damage  case 

Mode  #1 

1 

Mode  #2 

%equency  (Hz) 
Mode  #3 

Mode  #4 

Mode  #5 

Nominal 

13.88/  13.88'i' 

14.48  / 14.47 

48.41/48.40 

64.03  /  64.03 

67.46  /  67.52 

A 

13.94/13.88 

9.50  /  9.43 

48.52  /  48.40 

64.16/64.03 

65.91/64.98 

B 

13.47/13.16 

14.12/14.25 

35.65  /  34.25 

60.18/59.06 

65.86/65.90 

C 

13.97/13.88 

11.39/11.29 

48.53/48.40 

64.50  /  64.03 

59.90/58.51 

D 

13.21/13.16 

14.44  / 14.24 

36.68  /  35.92 

61.35/61.20 

66.95  /  65.97 

E 

13.96/13.88 

11.42/11.30 

48.57  /  48.40 

64.61  /  64.03 

59.91/58.60 

F 

13.94/13.88 

9.50  /  9.42 

48.50  /  48.40 

64.08  /  64.03 

65.80/64.91 

G 

12.29/12.26 

14.50/14.47 

48.67  /  48.40 

50.65  /  48.60 

67.76  /  67.52 

H 

13.73/13.70 

14.55  / 14.47 

48.68  /  48.40 

54.76  /  54.30 

67.71/67.52 

I 

13.74/13.58 

9.86  /  9.74 

36.66/35.89 

63.35  /  62.56 

58.86/57.46 

fData  presented  in  (measured  /  FEM  simulated  )  format,  where  the  first  number  represents  the  measured  frequency  and  the 
second  is  the  result  of  an  eigenanalysis  on  the  FEM  model  with  the  damaged  element(s)  removed. 


with  a  priori  knowledge  of  the  damage.  An  attempt  to  ensure  the  simulated  damaged  analytical  model 
correlated  well  with  the  measured  data  was  the  rationale  behind  maintaining  symmetry  during  the 
tuning  process. 

8.2.3  Sensor  Prioritization  and  Damage  Localization  Analysis.  After  tuning  the  analyt¬ 
ical  model,  a  prioritization  of  the  degrees-of-freedom  to  instrument  was  performed  as  developed  in 
Chapter  VII.  To  demonstrate  the  capabilities  of  the  APE  software,  a  small  number  of  sensors  (8)  were 
chosen.  Threshold  values  of  10%  for  detectability  and  7%  percent  for  colinearity  were  used  for  both  the 
sensor  prioritization  and  the  damage  localization  method.  These  threshold  values  represent  the  assumed 
combined  uncertainty  in  both  the  measurement  error  and  the  modeling  error.  These  values  were  chosen 
by  performing  three  analyses  for  threshold  values  of  5,  7,  and  10%,  and  then  comparing  overall  results 
against  values  in  Table  8.2.  The  eight  prioritized  sensor  locations  are  shown  in  Figure  8.6.  These  eight 
degree-of-freedom  locations  were  used  to  construct  the  eight  elements  of  the  partial  eigenvectors  for 
the  damage  identification  process. 

Having  identified  fhe  sensor  locafions,  a  damage  localizafion  analysis  was  performed.  The 
resulfs  of  fhis  analysis  are  confained  in  Table  8.3,  lisfing  fhe  undefecfable,  symmefric,  and  idenlifiable 
elemenfs.  Table  8.4  presenfs  a  descripfion  of  fhe  elemenf  numbering  used.  The  resulfs  show  fhaf  using 
only  fhe  firsl  five  modes  and  fhe  8  componenf  eigenvectors,  64  of  fhe  elemenfs  are  undefecfable  from  fhe 
measured  dafa.  This  indicates  fhaf  changes  in  fhe  measured  dafa  are  insignificanl  from  damage  in  fhese 
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Figure  8.6  Prioritized  sensor  locations  for  the  NASA  truss. 

elements.  These  results  are  consistent  with  a  similar  analysis  on  this  truss  presented  hy  Kashangaki, 
Smith  and  Lim''*’'  showing  that  95%  of  the  total  strain  energy  associated  with  the  first  six  modes 
was  contained  in  only  40  elements.  The  unidentifiahle  elements  are  categorized  as  either  battens,  or 
elements  located  near  the  free  end  of  the  truss.  The  remaining  40  elements  of  the  localization  analysis 
are  divided  among  23  symmetric  groupings  containing  one,  two  or  four  elements.  One  element  from 
each  of  the  23  symmetric  groups  is  used  to  define  fhe  inifial  search  space  for  fhe  idenfificafion  process. 


8.2.4  Damage  Identification  Results.  Using  fhe  resulfs  of  fhe  sensor  priorifizafion  fo  define 
fhe  measured  dafa,  fhe  funed  analytical  model,  and  fhe  damage  localization  analysis  fo  define  fhe 
inifial  search  space,  damage  identification  using  APE  was  performed.  The  resulfs  are  confained  in 
Table  8.5.  On  average,  fhe  resulfs  were  achieved  in  one  minufe  requiring  20  iferafions.  In  six  of  fhe 
nine  cases,  fhe  damage  was  localized  fo  a  single  elemenf  or  a  single  symmefric  grouping.  For  Case  D 
and  Case  E,  a  repealed  use  of  fhe  APE  melhod  using  fhe  resulfs  of  fhe  firsl  identification  application 
as  the  initial  search  space,  was  able  to  localize  the  damage  down  to  the  single  correct  element  or 
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Table  8.3  Damage  localization  results  for  the  NASA  8-bay  truss. 


Element  # 

Equivalent  Symmetric  Elements 

32 

36, 45,  49 

34 

38 

47 

51,60,64 

58 

62,71,75 

59 

61 

63 

65 

72 

73 

77,  86,  90 

74 

76 

78 

84 

88,  97,  101 

85 

87 

89 

91 

98 

99 

103 

100 

102 

104 

Undetectable  Elements 

1  -3 1 ,33,35,37,39-44,46,48,50,52-57,66-70,79-83,92-96 

Table  8.4  Element  numbering  and  descriptions  for  the  NASA  8-bay  truss. 


Bay^  # 

Eongeron 

Description 

Diagonal 

Batten + 

1 

6,  8,  10,  12 

7,  9,  11,  13 

1,2,  3,4,5 

2 

19,21,23,25 

20,  22,  24,  26 

14,  15,  16,  17, 

18 

3 

32,  34,  36,  38 

33,  35,  37,  39 

27,  28,  29,  30, 

31 

4 

45,47, 49,51 

46, 48,  50,  52 

40,41,42, 43, 

44 

5 

58,  60,  62,  64 

59,61,63,65 

53,  54,  55,  56, 

57 

6 

71,73,75,77 

72,  74,  76,  78 

66,  67,  68,  69, 

70 

7 

84,  86,  88,  90 

85,87,  89,91 

79,  80,  81,  82, 

83 

8 

97,  99,  101,  103 

98,  100,  102,  104 

92,  93,  94,  95, 

96 

fBays  are  numbered  consecutively  starting  from  the  free  end. 
^Includes  diagonal  members  in  the  batten  plane. 


single  symmetric  group.  A  capability  could  be  incorporated  into  the  APE  algorithm  to  adaptively  re¬ 
initialize  the  algorithm.  Also  noted  in  the  these  two  cases  is  assigned  damage  values  which  are  greater 
than  100%.  Unlike  ASTROS-ID,  the  APE  method  cannot  constrain  the  values  of  gi  to  lie  strictly 
within  zero  and  one.  These  constraints  can  only  be  enforced  by  removing  from  the  search  space  the 
corresponding  elements  with  ^/s  outside  the  allowable  range.  As  previously  discussed  in  Chapter  VI, 
to  avoid  prematurely  discarding  elements  and  allow  for  modeling  error,  the  allowable  range  of  g  should 
be  widened.  Eor  the  results  reported,  the  allowable  range  was  set  to  (0  <  gi  <  2).  An  additional 
capability  to  slowly  reduce  the  allowable  range  as  the  iteration  progresses  could  be  incorporated.  The 
exact  method  to  accomplish  both  the  re-initialization  and  allowable  range  reduction  are  referred  to  as 
algorithm  percolation  methods,  and  is  listed  as  a  topic  of  future  recommended  research.  Eor  damage 
Case  I,  the  compound  break,  damage  to  element  #71  was  not  identihed.  Damage  to  element  #71  (a 
longeron  in  the  sixth  bay)  was  assigned  to  a  longeron  in  either  bay  three  or  four.  This  difficulty  is 
in  part  due  to  the  fact  that  the  measured  data  for  this  damage  case  does  not  correlate  well  with  the 
simulated  damaged  analytical  model,  as  indicated  by  the  values  given  in  Table  8.2.  The  true  culprit, 
modeling  error  or  measurement  error,  cannot  be  determined  from  the  known  information.  Anytime  the 
simulated  damage  to  the  analytical  model  does  not  agree  with  the  measured  data  for  the  same  damage 
conhguration,  any  method  based  on  matching  the  partial  measured  data  will  have  difficulty  in  obtaining 
the  true  solution. 

Eor  comparison  purposes,  ASTROS-ID  was  performed  on  Cases  D,  H,  and  I  using  the  same  23 
element  initial  search  space  as  used  for  the  APE  method.  The  results  are  reported  in  Table  8.6.  In  the 
objective  function.  Equation  (3.2),  the  eigenvalues  were  ah  assigned  equal  weighting  of  100  for  ah 
hve  modes  (a,  =  100),  and  ah  eigenvector  components  were  assigned  unity  weighting  (bij  =  1).  Eor 
Case  D,  ASTROS-ID  identihed  damage  to  diagonal  elements  in  the  last  four  bays  of  the  truss,  but  did 
not  isolate  it  down  to  a  single  element.  Eor  Case  H,  damage  was  identihed  to  longeron  elements  in  the 
third  and  fourth  bay  of  the  truss.  As  for  APE,  ASTROS-ID  could  not  correctly  idenitify  damage  Case 
I,  due  to  the  measuement  and/or  modeling  error  as  previosuly  discussed.  Eor  Case  I,  both  methods  did 
identify  damage  in  both  longerons  and  diagonal  elements. 

Comparing  the  overall  results,  the  ASTROS-ID  method  tended  to  spread  the  assigned  damage 
over  several  elements,  whereas  the  APE  method  typically  assigned  damage  to  a  single  element  or 
symmetric  group.  This  comparison  illustrates  the  different  nature  of  the  two  solution  techniques,  as 
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Table  8.5  APE  identification  results  on  the  NASA  truss. 


True  damage 

APE  Identihed 

Damage  case 

Element  # 

%  damage 

Element  # 

%  damage 

cpu  time  (sec.) 

A 

84 

100 

84(88,97,101)t 

85 

33 

B 

85 

100 

85 

96 

80 

C 

71 

100 

58(62,71,75) 

94 

65 

D 

78 

100 

78 

98 

76 

104 

no 

E 

62 

100 

58(62,71,75) 

no 

56 

32(36,45,49) 

103 

P 

97 

100 

84(88,97,101) 

89 

32 

G 

51 

100 

47(51,60,64) 

88 

85 

H 

34 

100 

34(38) 

97 

43 

I 

71 

100 

32(36,45,49) 

80 

63 

78 

100 

78 

99 

fData  presented  in  I(S)  format  where  I  is  the  number  of  the  identified  element  and  S  is  the  symmetric  element  numbers. 

previously  discussed  in  Section  6.2.  Differences  in  the  solutions  are  also  in  part  due  to  the  method 
used  to  establish  the  initial  search  space  from  the  sensitivity  analysis.  Both  methods  used  the  results  of 
the  same  damage  localization  analysis,  however  the  analysis  as  described  in  Section  7.3  was  tailored 
for  the  APE  method.  In  general,  an  order  of  magnitude  increase  in  CPU  time  was  required  for  the 
ASTROS-ID  solution. 


8.3  Flexible  Truss  Experiment  (FTE) 

8. 3. 1  Hardware.  The  six-meter  PTE  was  assembled  at  the  Air  Porce  Institute  of  Technology 
from  excess  hardware  received  from  the  Structural  Dynamics  Branch  of  Wright  Eaboratory  after 
termination  of  the  12-Meter  Truss  Active  Control  Experiment. The  hardware  consists  of  the  truss 
assembly,  actuators  and  their  power  drivers,  accelerometer  sensors,  and  real-time  digital  control  and 
signal  processing  equipment.  As  presented  in  the  documentation  on  the  12-meter  truss,  there  was 
considerable  difficulty  in  obtaining  a  model  of  the  12-meter  truss  which  correlated  well  with  measured 
data.  Although  only  half  the  structure  is  currently  used  due  to  physical  space  limitations,  the  difficulty 
in  modeling  the  structure  makes  it  ideally  suited  to  use  in  validating  a  model  tuning  algorithm.  It 
should  also  be  noted  that  the  experiment’s  name  is  somewhat  of  a  misnomer.  Although  the  word  truss 
is  used,  the  structure  is  actually  a  frame  structure  with  rigid  connections  between  members.  The  term 
truss  typically  is  used  for  structures  with  pinned  connections. 
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Table  8.6  ASTROS-ID  identification  results  on  the  NASA  trass. 


True  damage 

ASTROS-ID  Identified 

Damage  case 

Element  # 

%  damage 

Element  # 

%  damage 

cpu  time  (sec.) 

D 

78 

100 

61 

65 

793 

63 

36 

65 

93 

76 

15 

78 

88 

91 

64 

104 

24 

H 

34 

100 

34(38) 

96 

602 

47(51,60,64) 

73 

I 

71 

100 

32(36,45,49) 

40 

541 

78 

100 

58(62,71,75) 

83 

61 

84 

63 

50 

65 

94 

78 

90 

84(88,97,101) 

94 

8. 3. 1.1  Truss  Description.  The  basic  structure  of  the  experiment  is  a  lightly  damped 
six-meter  trass,  cantilevered  vertically  from  a  rigid  support  base.  The  FTE  is  depicted  in  Figure  8.7. 
The  trass  is  composed  of  two  equal  length  frames  of  welded  tubular  aluminum  alloy  longerons  and 
battens  with  bolt-in  tubular  Lexan  diagonals  in  a  back-to-back  “K”  pattern.  A  section  of  the  FTE  is 
depicted  in  Figure  8.8.  The  assembled  trass  has  a  square  cross  section  of  20  inches  on  a  side.  The 
longerons  are  made  from  6061-T6  aluminum  alloy  tubes  with  a  1.5-inch-square  cross  section  and 
0.065  inch  wall  thickness.  The  battens  are  6061-T6  tubes  with  0.5-inch-square  cross  section  and  0.063 
inch  wall  thickness.  The  diagonal  members  are  Fexan  tubing  (270,000  psi  elastic  modulus)  with  a 
1.5-inch-diameter  circular  cross  section  with  0.125-inch  wall  thickness.  The  diagonals  have  aluminum 
end  fittings  which  are  fastened  to  the  trass  with  two  bolts  and  a  half-clevis  joint  at  both  ends.  The  truss 
has  4  bays  in  each  of  the  two  sections  for  a  total  of  eight  bays.  The  two  sections  are  bolted  together 
with  two  bolts  at  each  longeron  end.  Four  bolts  at  the  base  of  each  longeron  secure  the  trass  to  a  l-inch 
thick  aluminum  plate  which  is  securely  bolted  to  the  laboratory  floor.  The  bolt-in  diagonals  allow  quick 
structural  modifications  to  simulate  a  damaged  structure. 
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Figure  8.7  Six-Meter  Flexible  Truss  Experiment,  showing  initial  sensor  locations. 


Figure  8.8  Joint  construction  for  the  FTE. 
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Linear  Bearing  Linear  Bearing 

Figure  8.9  Linear  actuators  for  the  FTE. 

83.1.2  Actuators.  Input  excitation  to  the  truss  for  frequency  response  testing  is 
provided  through  the  use  of  two  linear  momentum  exchange  (reaction  mass)  actuators.  Each  actuator  is 
capable  of  approximately  1 -pound-force  peak  output.  Power  is  provided  to  the  actuators  hy  individual 
current  drive  circuits.  The  actuators  use  a  linear  dc  motor  with  the  armature  fixed  to  a  hase  and 
the  permanent  magnetic  field  suspended  on  shafls  and  linear  hearings.  Two  linear  springs  provide 
fhe  centering  force  for  fhe  mass.  The  resonanf  frequency  for  each  acfuafor  is  approximately  0.9  Hz, 
wifh  an  effeclive  viscous  damping  rafio  of  approximately  ten  percenf  of  crifical.  Eor  comparison,  fhe 
fundamenlal  frequency  of  fhe  fruss  is  approximately  7  Hz.  Mounfing  plafes  were  fahricaled  fo  attach 
fhe  acfuafors  atop  fwo  longerons  on  fhe  free-end  of  fhe  fruss.  Driving  fhe  acfuafors  in  phase  wifh  one 
anofher  excites  fhe  bending  modes  in  one  axis.  The  torsional  modes  are  excited  by  driving  fhe  acfuafors 
180  degrees  oul-of-phase.  A  skefch  of  fhe  acfuafor  is  shown  in  Eigure  8.9. 

8. 3. 1.3  Sensors  and  Supporting  Equipment.  Accelerafion  measuremenfs  are  made 
using  eighf  Sunsfrand  QA-1400  single-axis  inerfial  accelerometers.  These  accelerometers  were  chosen 
for  fheir  high  sensifivify  and  low  noise  characferisfics.  Inifially,  fhe  eighf  accelerometers  were  placed 
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Figure  8.10  Low  frequency  mode  shapes  for  the  FTE. 

along  one  longeron  at  the  node  location  at  each  bay.  For  signal  processing,  a  Tektronix  2642A  Fourier 
Analyzer  is  used  to  measure  and  average  the  frequency  response  functions. 

8.3.2  Model  Tuning. 

8.3.2. 1  Finite  Element  Model.  The  FTE  was  modeled  in  ASTROS  using  five  different 
types  of  beam  elements  for  a  total  of  96  elements.  The  five  types  were:  aluminum  vertical  longeron 
elements,  aluminum  horizontal  batten  elements,  Eexan  diagonal  elements,  aluminum  horizontal  mid¬ 
batten  elements,  and  aluminum  horizontal  top-batten  elements.  The  different  batten  configurations 
are  due  to  the  fact  that  the  ETE  is  constructed  of  two  4-bay  sections  which  can  be  bolted  together. 
Eumped  masses  were  included  to  account  for  the  actuators,  top-plates  (actuator  attach  points),  “K” 
brackets  (used  to  secure  the  diagonal  members)  and  mid-plates  (bolt  assemblies  used  to  secure  the 
two  4-bay  sections  together).  Care  was  taken  to  accurately  obtain  the  mass  properties,  area  properties 
and  moments  of  the  structural  elements.  The  elastic  properties  were  determined  from  laboratory  tests 
on  individual  elements.  These  parameters  were  all  used  to  construct  a  baseline  data  deck  for  input 
into  ASTROS.  Using  the  results  of  the  baseline  finite  element  analysis  along  with  software  written  in 
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MATLAB®to  display  and  animate  mode  shapes,  a  characterization  of  the  low  frequency  (helow  70  Hz) 
behavior  of  the  FTE  was  performed.  Mode  shapes  were  classified  into  four  categories  as:  X-hending, 
Y-hending,  torsion,  and  breathing.  An  illustration  of  these  four  different  mode  shapes  are  shown  in 
Figure  8.10.  Unlike  the  NASA  truss,  the  FTE  did  not  have  diagonal  members  in  the  batten  plane  and 
thus  exhibited  a  breathing  mode.  The  breathing  mode  corresponds  to  an  expansion  and  contraction  of 
the  truss  in  the  batten  plane,  similar  to  the  motion  of  the  human  rib  cage.  The  diagonal  battens  in  the 
NASA  truss  sufficiently  stiffened  the  batten  plane  to  preclude  this  type  of  motion  at  low  frequencies. 
The  presence  of  this  fourth  type  of  mode  increased  the  modal  density  and  consequently  increased  the 
difficulty  in  the  modal  parameter  identification  process.  Due  to  the  near  symmetry,  the  Y-bending  mode 
shapes  are  similar  to  the  X-bending  shapes  and  are  not  shown.  Note  that  the  frequencies  appearing 
in  the  figure  are  for  fhe  analyfical  model  prior  fo  funing.  The  nexf  step  was  fo  compare  fhe  analyfical 
resulfs  fo  measured  dafa. 

8.3.22  Testing  Procedure.  Experimenfal  measuremenfs  were  performed  on  fhe  FTE 
using  random  vibrafion  fesfing.  Measured  eigendafa  of  fhe  FTE  were  obfained  from  fhe  16  FRFs 
befween  fhe  fwo  linear  acfuafors  and  eighf  single-axis  accelerometers  placed  as  shown  in  Figure  8.7. 
The  inpuf  excifafion  fo  each  acfuafor  was  a  band  limited  (0  -  50  Herfz)  pseudo-random  signal.  There 
were  4096  discrete  sample  poinfs  recorded  by  fhe  specfrum  analyzer,  providing  a  frequency  resolufion 
of  0.031  Herfz.  Frequency  averaged  fransfer  funcfions  befween  fhe  inpuf  excifafion  and  fhe  eighf  ac¬ 
celerometers  were  measured.  The  inverse  discrete  Fourier  Iransforms  of  fhese  fransfer  funcfions  yielded 
fhe  impulse  response  funcfions  which  were  inpuf  info  ERA  fo  obfain  measured  modal  frequencies  and 
shapes,  as  described  in  Chapter  IV.  Disfinguishing  befween  bending  modes  and  torsion  modes  was 
easily  facililaled  by  excifing  fhe  sfrucfure  in  bofh  bending  (acfuafors  in  phase),  and  torsion  (acfuafors 
180  degrees  oul-of-phase).  Figure  8.11  shows  fhe  resulfing  averaged  FRF  oufpuf  af  accelerometer  #8 
for  fhe  fwo  differenl  types  of  excifafion,  measured  using  a  100  Hz  bandwidfh. 

8. 3.2. 3  Sensor  Prioritization.  The  sensor  priorifizafion  mefhod  presented  in  Chap¬ 
ter  VII  was  applied  to  fhe  FTE  to  determine  fhe  priorifized  locafion  of  fhe  eighf  single-axis  accelerom¬ 
eters  for  damage  idenfificafion.  During  inifial  fesfing,  if  was  determined  fhaf  fhe  inifial  locafions  of 
fhe  eighf  sensors,  along  one  longeron  as  was  shown  in  Figure  8.7,  was  nol  a  good  choice.  For  fhese 
sensor  locafions,  if  was  nof  possible  to  measure  fhe  Y-bending  mode.  Addifionally,  several  of  fhese 
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Figure  8.11  Comparison  of  bending  and  torsional  excitation  on  the  FTE. 

locations  yielded  small  gradient  values  as  determined  by  the  sensitivity  analysis.  For  the  analysis,  the 
first  eight  modes  were  used,  which  represented  all  the  global  modes  within  the  frequency  band  (0-50 
Hz)  measured.  These  modes  were  used  to  prioritize  the  eight  sensor  locations.  The  prioritized  locations 
of  the  eight  sensors  are  shown  in  Figure  8.12.  The  output  of  these  eight  sensors  was  used  to  obtain  the 
partial  mode  shapes  used  in  the  tuning  algorithm. 

8.3. 2. 4  Measured  Data.  Numerous  problems  were  encountered  in  obtaining  the 
measured  modal  properties.  The  first  was  in  the  choice  of  the  excitation.  The  original  choice  was  to 
use  the  two  actuators  in-phase  and  out-of-phase  as  previously  described.  This  presented  a  difficulty  in 
measuring  the  breathing  modes.  Examination  of  the  transfer  functions  depicted  in  Eigure  8.1 1,  illustrate 
the  absence  of  the  breathing  modes  (~22  Hz,  ~34  Hz)  in  the  measured  transfer  functions.  Although 
ASTROS-ID  does  not  require  that  all  modes  be  included  in  the  objective  function,  without  inclusion 
of  the  breathing  modes,  numerous  mode  swaps  ocurred  as  the  breathing  modes  “wandered”  during  the 
tuning  process.  Eor  damage  identification  to  be  successful,  it  is  important  that  the  simulated  analytical 
damage  closely  match  the  measured  data  of  the  damaged  structure  as  previously  discussed.  This  is  best 
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Figure  8.12  Prioritized  sensor  locations  on  the  FTE. 

accomplished  hy  ensuring  all  modes  including  the  breathing  mode  match  the  physical  structure  as  hest 
as  possible.  An  additional  discussion  on  this  topic  is  contained  in  Reference  11. 

To  measure  the  breathing  mode,  the  actuators  atop  the  longerons  were  rotated  ninety  degrees. 
This  allowed  direct  excitation  of  the  breathing  mode  by  exciting  out-of-phase,  and  X-bending  excitation 
by  exciting  in-phase.  In  this  configuration  however,  both  Y-bending  and  torsion  modes  were  not  excited 
and  hence  could  not  be  accurately  measured.  The  simple  solution  seemed  to  be  to  take  two  measurement 
sets  in  the  first  excitation  configuration,  and  two  additional  measurement  sets  with  the  actuators  rotated 
ninety  degrees.  While  all  modes  were  now  clearly  identifiable  in  fhe  fransfer  funclions,  a  new  problem 
was  observed.  The  inerfia  properfies  of  fhe  recfangular  linear  acfuafor’s  affecfed  fhe  measured  response 
depending  on  fheir  orienfafion,  and  fhus  a  shiffing  of  fhe  modes  resulfed  whenever  fhe  acfuafor  posifion 
was  changed.  A  compromise  was  achieved  by  using  fhe  fwo  acfuafors  afop  fwo  longerons  diagonally 
from  one  anofher  and  mounfed  ninefy  degrees  fo  one  anofher.  In  fhis  configurafion  each  acufafor  was 
excifed  independenfly,  one  excifing  X-bending  and  breafhing,  and  fhe  second  excifing  Y-bending  and 
forsion.  This  sfill  did  nof  represenf  an  ideal  sifuafion  for  fhe  following  reason.  As  one  acfuafor  was 
being  excifed,  fhe  ofher  acfuafor’s  proof  mass  was  fixed.  This  effecfively  changed  fhe  apparenf  inerfia 
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properties  atop  either  longeron,  depending  upon  which  actuator  was  being  excited.  Again,  a  small 
shift  was  observed  in  the  measured  frequencies,  however,  this  was  determined  to  be  an  acceptable 
compromise,  not  requiring  actuator  re-design.  Note  that  the  difficulties  with  the  actuator  excitation 
all  arose  as  the  result  of  using  non-ground  based  excitation,  which  represented  a  significant  mass  at 
the  tip  of  the  cantilevered  frame.  For  future  space  application  studies,  actuator  design  is  an  important 
consideration  in  obtaining  accurate  modal  data. 

Having  finally  determined  fhe  sensor  and  acfuafor  posifions,  fhe  modal  exfracfion  mefhod  using 
ERA  was  performed.  Due  fo  fhe  closely  spaced  modes,  an  addifional  problem  was  encounfered. 
Allhough  fhe  modal  frequencies  were  easily  idenlified,  fhe  exlracled  mode  shapes  were  observed  as 
linear  combinalions  of  one  anolher  for  fhe  closely  spaced  modes.  This  was  especially  Iroublesome 
during  fhe  idenlificalion  of  parlial  mode  shapes  for  fhe  damaged  slruclure.  Damage  Case  2,  as  described 
in  Seclion  8.3.4  was  chosen  lo  illuslrale  fhe  problem.  The  measured  modes  along  wilh  fhe  analylical 
modes  for  fhe  eighl  sensor  localions  are  depicted  in  Figure  8.13.  As  shown  in  fhe  figure,  Ihere 
is  considerable  discrepancy  in  fhe  dala,  parlicularly  for  fhe  closely  spaced  fourlh,  fiflh,  and  sixlh 
modes.  Addifional  lesling  using  four  independenl  excilalions,  where  each  mode  was  individually 
excited,  yielded  close  resulls  lo  fhe  analylical  shapes.  However,  when  only  Iwo  aclualors  were  used, 
for  fhe  reasons  expounded  upon  above,  fhe  modes  appeared  as  linear  combinalions  of  one  anolher. 
To  delermine  fhe  linear  combinalion,  a  leasl  squares  solulion  was  used  lo  determine  Ihe  coeflicienls 
required  lo  besl  malch  Ihe  measured  modes  lo  Ihe  modes  of  Ihe  simulated  damaged  analylical  model. 
Two  coeflicienls  were  computed  lo  separate  modes  1  and  2,  and  Ihree  differenl  coefficienls  were  used 
lo  separate  modes  3,  4,  and  5.  The  de-coupled  mode  shapes  are  shown  in  Figure  8.14.  The  abscissa 
represenls  Ihe  eighl  sensor  measuremenls  numbered  one  Ihrough  eighl.  A  line  is  shown  connecling 
Ihe  eighl  discrete  poinls  lo  aid  in  visualizalion.  Note  lhal  Ibis  de-coupling  melhod  was  only  possible 
because  of  Ihe  known  localion  of  Ihe  damage  for  Ihe  simulalion.  In  aclual  praclice,  a  re-design  of 
Ihe  inpul  aclualors  would  be  required  lo  correclly  oblain  Ihe  measured  modes  wilhoul  resorling  lo 
Ihe  analylical  model.  This  problem  was  nol  observed  in  Ihe  dala  on  Ihe  NASA  Imss,  which  did  nol 
use  non-ground  based  aclualors,  nor  exhibil  a  brealhing  mode.  In  all  cases,  Ihere  was  no  Irouble  in 
oblaining  Ihe  torsion  shapes,  which  were  clearly  separated  from  Ihe  bending  and  brealhing  shapes. 
Addifional  research  on  determining  parlial  mode  shapes  for  closely  spaced  modes  is  recommended  for 
fulure  work. 
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Figure  8.13  Analytical  and  measured  mode  shapes  before  de-coupling,  Case  2. 
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Figure  8.14  Analytical  and  measured  mode  shapes  after  de-coupling,  Case  2. 
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Figure  8.15  Comparison  of  initial  analytical  model  and  measured  frequency  response. 

83.2.5  Model  Tuning  Results.  Despite  the  effort  expended  in  constructing  the  baseline 
finite  element  model  of  the  FTE,  it  did  not  agree  well  with  the  measured  data.  The  poor  correlation 
between  the  baseline  analytical  and  measured  data  is  depicted  in  Figure  8.15,  which  shows  the  transfer 
function  between  a  colocated  sensor  and  actuator  at  the  top  of  the  FTE.  Prior  to  computing  the  analytical 
transfer  functions,  damping  was  included  in  the  model  by  assuming  a  value  of  one  half  of  one  percent  of 
critical  damping  for  all  modes  (a  typical  measured  value).  This  was  done  only  to  avoid  the  unbounded 
resonant  spikes  due  to  an  undamped  model,  and  thus  the  height  of  the  resonant  peaks  between  analytical 
and  measured  data  is  insignificant.  It  is  clear  from  Eigure  8.15  that  the  baseline  finite  element  model 
does  not  adequately  represent  the  measured  dynamic  characteristics  of  the  ETE  and  would  hence  benefit 
greatly  from  model  tuning. 

Using  ASTROS-ID,  the  baseline  finite  element  model  was  tuned  to  the  measured  data.  Eor  the 
tuning,  the  frequencies  and  shapes  of  the  first  eight  modes  were  included  in  the  objective  function. 
Twenty  three  design  variables  were  chosen  for  the  tuning  process.  A  description  of  the  design  variables 
as  well  as  the  tuned  values,  are  given  in  Table  8.7.  Values  in  the  table  are  normalized,  such  that  the 
initially  assumed  value  for  each  design  variable  is  unity.  To  account  for  the  added  stiffness  in  the 


8-20 


Table  8.7  Design  variable  values  in  the  FTE  finite  element  model. 


DV# 

Descripfion 

Tuned  value 

1 

baffens  (elasfic  modulus) 

1.00 

2 

mid-baltens(elaslic  modulus) 

1.00 

3 

longerons  (elasfic  modulus) 

0.80+ 

4 

diagonals  (elasfic  modulus) 

1.20+ 

5 

fop-ballens  (elasfic  modulus) 

0.99 

6 

baffens  (/i) 

0.79 

7 

mid-baffens  (/i) 

0.67 

8 

longerons  (/i) 

0.58 

9 

diagonals  (/i) 

1.18 

10 

fop-ballens  (/i) 

0.84 

11 

baffens  {I 2) 

1.11 

12 

mid-ballens  {I 2) 

1.14 

13 

longerons  {I 2) 

1.99 

14 

diagonals  {I 2) 

1.00 

15 

fop-ballens  {I 2) 

1.00 

16 

baffens  (mass) 

1.51 

17 

mid-ballens  (mass) 

1.26 

18 

longerons  (mass) 

0.42 

19 

diagonals  (mass) 

2.81 

20 

“K”  brackels  (mass) 

0.50+ 

21 

aclualors  (mass) 

1.00 

22 

mid-plales  (mass) 

0.84 

23 

lop-plales  (mass) 

1.28 

fValues  were  normalized  such  that  initial  values  are  all  unity. 
fValue  on  boundary  of  allowable  excursion  limit  for  this  variable. 


joints  resulting  from  the  welded  assemblies,  the  bending  stiffness  of  the  elements  was  chosen  as  a 
design  variable  by  allowing  the  Ii  and  I2  properties  to  vary  in  addition  to  the  elastic  modulus  for  each 
element.  To  account  for  the  symmetry  of  the  structure,  all  common  elements  were  linked  to  a  single 
design  variable.  The  tuning  process  converged  in  six  outer-loop  iterations  requiring  18  minutes  of 
CPU  time.  Convergence  was  defined  as  a  less  fhan  one  half  of  one  percenf  change  in  fhe  objecfive 
funclion  value  befween  fwo  consecufive  ouler-loop  iferafions.  The  resulfs  of  fhis  Inning  are  shown  in 
Figure  8.16,  depicted  for  fhe  same  acfuafor/sensor  combinafion  as  used  in  Figure  8.15. 

The  resulfs  show  excellenf  agreemenf  over  fhe  40  Hz  frequency  band  spanned  by  fhe  modes 
included  in  fhe  objecfive  funcfion.  The  firsf  four  columns  of  Table  8.8  lisfs  fhe  frequencies  for  fhe 
firsl  eighf  modes,  comparing  fhe  measured  dafa  fo  fhe  funed  finile  elemenf  model.  The  objecfive 
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Figure  8.16  Comparison  of  the  tuned  analytical  model  and  measured  frequency  response. 

function.  Equation  (3.2),  included  the  terms  a,  and  bij  to  weight  individual  contrihutions  to  the  overall 
value.  For  the  tuning,  all  a,  coefficients  were  set  to  100  and  all  bij  were  set  to  unity.  This  represents 
equal  confidence  in  fhe  measuremenf  dafa  for  all  frequencies  included  in  fhe  ohjecfive  funcfion,  wifh  a 
sfronger  emphasis  placed  on  funing  fhe  modal  frequencies  fhan  fhe  shapes.  In  general,  selecfion  of  fhe 
weighfing  coefficienfs  is  dependenf  upon  fhe  confidence  wifh  which  each  mode  is  measured,  as  well 
as  fhe  designer’s  desire  fo  minimize  a  selecfed  porfion  of  fhe  correlafion  error.  For  fhis  case,  emphasis 
was  placed  on  minimizing  fhe  frequency  correlafion  errors. 

8.3.3  Damage  Localization  Analysis.  Using  fhe  funed  analyfical  model,  fhe  damage  lo- 
calizafion  mefhod  presenfed  in  Chapfer  VII  was  applied  fo  fhe  FTE.  Based  on  repeafahilify  of  fhe 
measuremenf  dafa,  and  fhe  exisfing  correlafion  errors  hefween  fhe  finife  elemenf  model  and  fhe  mea¬ 
sured  dafa,  fhe  defecfahilify  fhreshold  values  were  sef  af  10%.  The  fhreshold  for  Sx  was  sef  af  7% 
and  fhe  5^  fhreshold  af  25%.  This  reflecfed  fhe  increased  confidence  in  measuring  fhe  eigenvalues 
over  fhe  eigenvectors,  reflecfing  fhe  difficulfy  in  measuring  fhe  closely  spaced  modes.  Using  fhe  eighf 
sensor  locafions  chosen  in  Secfion  8. 3. 2. 3  and  fhe  firsf  eighf  modes,  fhe  exfenf  fo  which  damage  can  he 
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Table  8.8  Analytical  and  measured  natural  frequencies  of  the  FTE. 


Mode  # 

Descripfion 

Nominal 

Erequency  (Hz) 

Dam 

Case  1 

aged 

Case  2 

EE  funed 

measured 

EE 

measured 

EE 

measured 

1 

Isf  Y  Bending 

6.32 

6.32 

6.32 

6.34 

6.23 

6.20 

2 

Isf  X  Bending 

6.33 

6.33 

5.85 

5.86 

6.32 

6.30 

3 

Isf  Torsion 

12.51 

12.64 

11.38 

11.27 

12.27 

12.37 

4 

Isf  Breafhing 

22.26 

22.29 

22.26 

23.31 

22.26 

22.10 

5 

2nd  Y  Bending 

23.97 

24.02 

22.01 

21.48 

23.17 

23.48 

6 

2nd  X  Bending 

24.09 

24.18 

24.03 

24.15 

24.04 

24.20 

7 

2nd  Breafhing 

34.01 

34.03 

33.36 

33.63 

33.71 

33.79 

8 

2nd  Torsion 

36.46 

36.31 

36.07 

35.89 

36.15 

36.02 

localized  was  determined  by  classifying  the  96  elements  using  the  U/S/I  format.  The  results  indicate 
that  26  elements  cannot  be  identified,  and  the  remaining  70  elements  are  arranged  in  15  symmetric 
groups.  The  26  undetectable  elements  were  all  batten  elements.  The  remaining  15  symmetric  groups 
consisted  of  either  longerons,  battens,  or  diagonal  elements  localized  to  either  a  single  bay  or  adjacent 
bays.  This  analysis  was  used  to  define  fhe  inifial  search  space  for  a  damage  idenfificafion  algorifhm 
using  APE.  The  resulfs  of  fhe  damage  localizafion  analysis  are  given  in  Table  8.9.  A  descripfion  of  fhe 
elemenf  numbering  is  confained  in  Table  8.10. 

83.4  Damage  Identification  Results.  Eor  fhe  APE  mefhod,  fhe  funed  elemenfal  parameters 
of  fhe  finite  elemenf  beam  model  were  used  fo  consfrucf  fhe  mafrix  B  according  fo  Equafion  (6.13) 
along  wifh  fhe  funed  mass  and  sfiffness  mafrices  M  and  K .  The  measured  parfial  eigendafa  of  modes 
1  fhrough  8  were  fhen  used  fo  determine  sfrucfural  damage.  Two  damage  configurafions  were  fesfed. 
The  firsf  (Case  1)  was  fhe  full  removal  of  a  diagonal  sfruf  in  fhe  sixfh  bay  (as  measured  from  fhe  free 
end).  The  second  (Case  2)  damage  configuralion  was  fhe  replacemenf  of  a  diagonal  sfruf  in  fhe  sevenfh 
bay  wifh  one  in  which  ~  50  percenf  of  fhe  sfruf’s  cross  secfional  area  was  removed.  The  analyfical 
and  measured  nafural  frequencies  for  fhe  firsf  eighf  modes  are  confained  in  Table  8.8.  Mismafches 
befween  fhe  analyfical  and  measured  frequencies  are  fhe  resulf  of  bofh  measuremenf  uncerfainly  as 
well  as  modeling  errors. 

Using  fhe  resulfs  of  fhe  damage  localizafion  analysis,  fhe  search  space  consisted  of  fifteen 
elemenfs,  one  elemenf  from  each  APE  symmefric  group.  Since  no  a  priori  knowledge  of  fhe  damage 
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Table  8.9  Damage  localization  results  for  the  FTE. 


Element  # 

Equivalent  Symmetric  Elements 

17 

18,  19,  20 

29 

30,31  32 

33 

34,  35,  36 

37 

39 

38 

40, 42,  44 

41 

43, 46,  48 

45 

47, 49,  52 

50 

51,53,55 

54 

56,  57,  59 

58 

60,  62,  63 

69 

70,71,72,  73,74,  75,76 

77 

78,  79,  80 

81 

90,  95,  96 

82 

89,  94,  97 

83 

84,  85,  86,  87,  88,  91,  92,  93,  98,  99,  100 

Undetectable  Elements 

5  -  16,21  -28,61,64-68 

Table  8.10  Element  numbering  and  descriptions  for  the  ETE. 


Bay^  # 

Eongeron 

Description 

Diagonal 

Batten 

1 

65,  66,  67,  68 

85,  86,91,  100 

33,  34,  35,  36 

2 

61,62,63,64 

84,  87,  92,  99 

29,  30,31,32 

3 

57,  58,  59,  60 

83,  88,  93,  98 

25,  26,  27,  28 

4 

53,  54,  55,  56 

82,  89,  94,  97 

21,22,  23,24 

5 

49,  50,51,52 

81,90,  95,96 

17,  18,  19,  20 

6 

45, 46,  47, 48 

77,  78,  79,  80 

13,  14,  15,  16 

7 

41,42,  43,44 

73,  74,  75,  76 

9,  10,  11,  12 

8 

37,  38,  39, 40 

69,  70,71,72 

5,  6,  7,  8 

fBays  are  numbered  consecutively  starting  from  the  free  end. 
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Figure  8.17  APE  identified  damage  regions  on  the  EXE. 
Table  8.1 1  APE  damage  identification  results  on  the  ETE. 


True 

APE 

APE 

cpu 

Test 

failed 

identified 

symmetric 

time 

Case 

element  # 

element  #s 

element  #s 

(sec.) 

1 

77  -  100% 

77  -  104% 

78  79  80 

530 

2 

75  -  50% 

69  -  40% 

70  71  72  73  74  75  76 

543 

location  is  assumed,  and  the  measurement  and  model  uncertainty  for  both  test  cases  are  identical,  the 
search  space  for  each  test  case  is  identical. 

The  results  of  the  two  test  cases  are  presented  in  Table  8.11  and  depicted  in  Eigure  8.17.  In 
each  case,  the  damage  was  correctly  localized  to  a  small  area  of  the  truss  containing  the  true  damaged 
element.  Eurther  damage  identification  refinement  to  the  exact  element  would  require  either  additional 
measured  data,  or  a  closer  correlation  between  measured  data  and  the  analytical  model.  Improving 
model  correlation  requires  a  higher  fidelity  model  and/or  less  measurement  uncertainty.  In  each  test 
case,  the  correct  element  and  exact  amount  of  damage  was  determined  when  noise  free  analytical 
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Figure  8.18  Effect  of  structural  damage  on  the  measured  response. 

simulated  data  was  used  in  lieu  of  the  measured  data.  A  representative  frequency  response,  showing 
the  effect  of  damage  (Case  1)  is  depicted  in  Figure  8.18. 

8.4  Summary 

The  damage  identification  process  was  illustrated  on  two  experimental  structures  and  the  results 
reported.  The  first  was  a  cantilevered  truss  modeled  with  104  rod  elements  with  a  total  of  96  degrees- 
of-freedom.  The  measured  data  consisted  of  the  first  five  flexible  modal  frequencies,  and  only  eighf 
componenfs  of  fhe  five  corresponding  eigenvectors.  The  second  experimenf  was  a  canfilevered  frame 
assembly  modeled  wifh  96  beam  elemenfs  wifh  a  fofal  of  192  degrees-of-freedom.  The  measured 
dafa  consisted  of  fhe  flrsf  eighf  flexible  modal  frequencies,  and  only  eighf  componenfs  of  fhe  eighf 
corresponding  eigenvectors.  In  each  fesf  case,  fhe  sfrucfural  damage  could  be  localized  to  a  small 
porfion  of  fhe  sfrucfure.  The  exfenf  to  which  damage  can  be  localized  was  limited  by  bofh  model 
fldelify  and  accuracy  of  fhe  measured  modes. 
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IX.  Conclusions  and  Recommendations 


9. 1  Research  Conclusions 

A  method  was  presented  to  identify  damaged  structural  elements  from  limited  measurement 
data.  The  problem  was  broken  into  four  distinct  tasks:  identification  of  modal  parameters,  model 
tuning,  damage  localization,  and  damage  identification.  The  research  showed  that  structural  damage 
can  be  identified  using  only  a  small  measured  subsef  of  fhe  eigensfrucfure.  Experimenfal  fesfs  were 
conducted  on  fwo  separate  sfrucfures.  The  firsl  was  a  canfilevered  fruss  modeled  wifh  104  rod  elemenfs 
wifh  a  fofal  of  96  degrees-of-freedom.  The  measured  dafa  consisted  of  fhe  firsl  five  flexible  modal 
frequencies,  and  only  eighl  componenls  of  fhe  flve  corresponding  eigenvectors.  The  second  experimenl 
was  a  canfilevered  frame  assembly  modeled  wifh  96  beam  elemenfs  wifh  a  lolal  of  192  degrees-of- 
freedom.  The  measured  dafa  consisted  of  fhe  firsl  eighl  flexible  modal  frequencies,  and  only  eighl 
componenls  of  fhe  eighl  corresponding  eigenvectors.  Two  key  factors  in  fhe  ability  lo  idenlify  fhe 
damage  are:  fhe  accuracy  of  fhe  measured  dafa,  and  fhe  fidelity  of  fhe  analylical  model.  In  fhe  case  of 
perfecl  measuremenls  and  perfecl  model  correlalion,  fhe  exacl  damage  can  be  idenlifled.  To  accounl 
for  imperfecl  measuremenls  and  model  correlalion  errors,  fhe  concepl  of  damage  localizalion  was 
inlroduced.  In  such  cases,  damage  can  be  localized  to  a  small  sub-seclion  of  fhe  slruclure. 

Also  invesligaled  was  fhe  relalion  belween  increasing  fhe  number  of  sensors  and  increasing 
fhe  number  of  measured  modes.  The  resulls  showed  lhal  more  informalion  on  slruclural  damage  is 
gained  from  fhe  ability  lo  measure  an  addilional  modal  frequency  lhan  from  fhe  abilily  to  measure 
one  more  degree-of-freedom.  Previous  sludies  invesligaled  which  modes  should  be  used  in  a  damage 
idenliflcalion  algorilhm  by  determining  which  modes  change  Ihe  most  In  Ihis  research  work,  il  was 
found  lhal  as  many  modes  as  can  accurately  be  measured  should  be  used.  Modes  which  are  unchanged 
from  damage  conlain  informalion  on  elemenfs  which  are  undamaged,  and  Ihus  help  furlher  narrow 
down  Ihe  search  space. 

For  Ihe  case  where  slruclural  damage  is  confined  to  changes  in  Ihe  sliffness  of  slruclural  elemenfs, 
il  was  shown  lhal  Ihe  resulling  nalural  frequencies  from  damage  can  only  decrease.  Placing  an 
upper  limil  on  Ihe  largel  frequencies  corresponding  to  Ihe  nominal  frequency  of  Ihe  analylical  model 
minimized  Ihe  problem  of  mode  swilching  during  Ihe  ileralion  process  and  produced  resulls  which 
were  closer  to  Ihe  Irue  damage. 
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An  objective  function  was  introduced  which  represents  the  mismatch  between  the  measured  and 
analytical  partial  eigendata.  Two  methods  were  used  to  minimize  this  function:  ASTROS-ID,  and  APE. 
A  comparison  of  these  two  methods  follows.  Both  solution  methods  use  an  iterative  solution  technique 
to  minimize  the  same  objective  function  and  are  suitable  with  partial  measured  modal  data.  (An  iterative 
method  is  not  required  for  the  APE  solution  when  perfect  full  length  eigenvector  measurements  can  be 
obtained.)  Eor  ASTROS-ID,  design  variables  can  be  any  elemental  mass  or  stiffness  values,  whereas 
for  APE  the  design  variables  are  limited  to  changes  in  elemental  stiffness.  The  advantage  of  the  APE 
method  is  that  it  does  not  require  an  eigenanalysis  and  the  computation  of  the  sensitivity  values  to  each 
design  parameter  at  each  outer-loop  iteration  step.  APE  however  requires  the  decomposition  of  a  large, 
sparse,  possibly  singular  matrix,  which  requires  sparse  matrix  techniques  to  make  it  computationally 
competitive.  Both  techniques  benefit  greatly  from  an  initial  sensitivity  analysis  to  reduce  the  search 
space  prior  to  initiating  either  method.  Due  to  the  algorithm  percolation  process  used  in  the  APE 
method,  APE  typically  reduces  the  solution  space  to  fewer  elements  than  the  ASTROS-ID  method.  As 
currently  coded,  the  APE  method  is  computationally  an  order  of  magnitude  faster  than  ASTROS-ID. 
The  success  of  either  method  is  contingent  upon  accurate  measured  data  and  an  analytical  model  which 
correlates  well  with  the  measured  data  for  the  nominal,  as  well  as  the  damaged  configuration.  Neither 
method  is  guaranteed  to  converge  to  a  global  minimum. 

9.2  Recommendations  for  Additional  Research 

During  the  course  of  any  research  investigation,  additional  understanding  of  the  problem  is  always 
accompanied  with  additional  questions.  During  the  experimental  portion  of  modal  identification  on 
the  ETE,  considerable  difficulty  was  associated  with  extracting  partial  modal  data  for  modes  which  are 
spaced  very  closely  in  frequency.  The  use  of  additional  input/output  relationships  helped  to  minimize 
this  problem,  but  may  not  be  possible  on  orbit.  Thus,  it  is  recommended  that  additional  research  be 
conducted  on  extracting  mode  shapes  for  closely  spaced  modes  using  non-ground  base  actuators.  Eor 
model  tuning,  several  methods  can  be  adopted  to  speed  up  the  ASTROS-ID  algorithm.  Currently,  the 
design  variable  sensitivity  is  performed  through  a  finite-difference  method.  Since  these  sensitivities  are 
known  explicitly  from  the  finite  element  formulation,  the  finite-difference  method  can  be  replaced  with 
either  an  analytical  method  in  the  current  version  of  ASTROS  (12)  or  the  use  of  the  matrix  operator  P 
as  used  in  APE.  Additional  research  can  be  conducted  in  the  percolation  of  either  algorithm  for  damage 
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identification.  For  either  method,  an  adaptive  method  to  narrow  the  search  space  can  he  employed  to 
speed  the  processing  time.  The  ASTROS-ID  search  space  currently  is  not  reduced  during  the  solution 
process,  and  hence  requires  the  sensitivity  calculation  during  each  outer-loop  iteration  step  for  all 
design  variables  included  in  the  initial  search  space.  For  the  APE  method,  the  search  space  is  reduced 
during  the  iteration  process  hy  fixed  initial  tolerances  on  the  damage  fractions.  A  method  of  gradually 
tightening  these  tolerances  as  the  solution  progresses  is  desirable,  and  requires  further  investigation. 
With  respect  to  damage  localization  and  sensor  prioritization,  the  results  are  dependent  on  the  chosen 
thresholds  for  both  detectability  and  colinearity.  Although  the  threshold  values  are  problem  dependent, 
additional  studies  would  enhance  future  selection  of  these  design  parameters.  Both  methods  can  be 
expanded  to  include  structural  damping  in  the  formulation. 

Lastly,  the  damage  identification  method  presented  in  this  study  should  be  applied  to  other 
structural  applications.  This  method  is  suitable  to  any  modeled  structure  where  practicality  dictates 
that  only  a  small  portion  of  the  eigenstructure  can  be  measured.  The  method  is  ideally  suited  for  remote 
monitoring  where  more  conventional  non-destructive  testing  methods,  such  as  x-raying  and  acoustic 
emissions,  cannot  be  employed. 
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Appendix  A.  ASTROS-ID  Software  Modules 


A  full  description  of  the  initial  version  of  ASTROS-ID  is  contained  in  the  work  hy  Gihson. 
However,  during  the  course  of  this  research  effort,  several  shortcomings  of  the  software  were  discovered, 
requiring  either  modification  or  replacement  of  the  original  software  code.  The  overall  intent  has  not 
been  altered  from  that  originally  proposed  hy  Gihson.  This  section  is  intended  only  to  provide  a 
brief  description  of  the  software  component  modules  used  in  ASTROS-ID  that  were  modified  or 
alfered  affer  fhe  original  work.  Errors  encounfered  using  fhe  original  code  included  runfime  errors 
when  using  non-consecufively  numbered  modes,  as  well  as  when  using  a  differenl  number  of  mode 
shapes  fhan  modal  frequencies.  Ofher  problems  encountered  included  gradienf  calculafions  fhaf  did 
nol  properly  accounf  for  mode  normalizafion  and  mode  swifches  which  were  nol  idenlified,  as  well  as 
errors  in  fhe  reporfing  of  resulfs.  Second-order  eigenvalue  gradienf  informafion  was  also  included  in 
fhe  sensifivify  calculafions.  Changes  in  fhe  soffware  are  listed  below  which  corrected  fhese  problems. 
The  mafhemafical  foundalion  was  presented  in  Chapter  V.  This  lisf  should  be  used  fo  supplemenf  fhe 
original  work.  Informafion  on  fhe  use  of  ASTROS  can  be  obfained  in  Reference  55. 

A.  1  Modified  Software  Modules 

1.  MAPOL-ID.SEQThe,  mapol  sequence  was  modified  fo  incorporate  fhe  non-consecufive 
eigendafa,  correcf  for  poinf  normalizafion,  and  incorporate  fhe  use  of  second-order  gradi¬ 
enf  s. 

2.  TUNE. FOR  This  module  was  modified  fo  incorporate  fhe  non-consecufive  eigendafa  and 
incorporate  fhe  use  of  second-order  gradienfs.  Also  included  was  fhe  ability  fo  sef  fhe 
opfimizafion  parameters  from  fhe  inpuf  dafa  deck,  and  incorporate  fhe  numerical  opfimizer 
Design  Opfimizafion  Tools  (DOT)  for  fhe  inner-loop  opfimizafion. 

3.  REPORT.FOR  This  module  was  updated  fo  correcf  numerous  formaffing  errors  when 
reporfing  resulfs. 

4.  Misc.  The  original  code  included  an  unsfrucfured  dafa  base  enfry  called  ‘MTRACE’ 
which  kepf  frack  of  fhe  mapping  befween  mode  numbers  during  fhe  iferafion  process.  A 
simpler  mefhod  of  renumbering  fhe  stored  inpuf  dafa  whenever  a  mode  swifch  occurred 
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was  incorporated.  As  a  result,  most  subroutines  required  a  minor  modification  to  remove 
the  MTRACE  reference. 


A.  2  New  Software  Modules 

1 .  ORTHTESTFOR  This  routine  checks  for  mode  switching  by  performing  a  matrix  search  on 
the  modal  correlation  coefficient  matrix,  and  updates  the  TSHAPE  and  TEREQ  relational 
entities  corresponding  to  switched  modes.  This  routine  replaced  the  previously  used 
orthogonalization  test. 

2.  MTCHINDX.FOR  This  routine  matches  mode  numbers  and  index  numbers  between  fre¬ 
quencies  and  shapes  required  when  non-consecutive  data  is  input. 

3.  GETDLAM2.FOR  This  routine  retrieves  entities  required  for  the  second-order  gradient 
term  calculations. 

4.  DLAM2.FOR  This  routine  computes  the  second-order  eigenvalue  gradient  terms. 

5.  UTMCOPYFOR  Matrix  utility  routine  to  copy  a  matrix  entity  in  the  data  base. 

6.  PUTSCAL.FOR  This  routine  allows  multiple  calls  to  put  scalars  into  a  matrix  in  the  data 
base  from  within  the  same  loop. 

7.  NORMINDX.FOR  This  routine  determines  the  internal  degree-of-freedom  corresponding 
to  the  normalization  point. 

8.  NORMAL.FOR  This  routine  point  normalizes  the  analytical  eigenvectors  such  that  the 
degree-of-freedom  corresponding  to  the  max  measured  degree-of-freedom  from  the  input 
(TSHAPE)  data  is  unity.  A  warning  message  is  displayed  (and  the  proper  normalization 
performed)  if  the  input  mode  shapes  were  not  normalized  to  maximum  entry  equals  one. 
Eigenvectors  not  used  in  the  tuning  process  are  mass  normalized. 
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Appendix  B.  APEWARE:  Damage  Identification  Integrated  Software 


Integrated  Damage  Identification  Toolbox 

Structural  Damage 

IdentificatioD 

using 

P  artial 
Pigenstractnre 

Toolbox  for  Matlab 

Version:  0,1  Beta,  January,  1996 

Richard  G.  Cobb,  Capt,  USAF 
Air  Force  Institute  of  Technology 

Start  Help 


Figure  B.l  Software  opening  menu. 

The  four  tasks  associated  with  damage  identification  were  combined  into  an  integrated  software 
tool  which  was  programmed  in  MATLAB®.  This  software  package  provides  the  design  engineer  an 
easy  to  use  tool  to  aid  in  determining  system  identification  requirements,  sensor  placement,  damage 
localization  studies,  and  damage  identification  using  the  assigned  partial  eigenstructure  method.  Also 
included  is  pre  and  post-processing  routines  for  use  with  ASTROS-ID.  This  appendix  is  intended  only 
to  provide  the  reader  with  a  brief  introduction  to  the  functionality  and  capabilities  of  the  software 
package.  The  algorithms  for  each  task  are  based  on  the  work  presented  in  the  main  text,  and  were 
briefly  described  under  the  software  implementation  section  for  each  task.  The  specific  algorithm 
descriptions  are  not  provided.  Familiarity  with  MATLAB®is  assumed.  Information  on  the  use  of 
MATLAB®is  contained  in  Reference  52.  The  software  was  written  and  implemented  on  a  Sun  Sparc- 
10  workstation  running  under  Unix/Sun4.  While  most  routines  are  ascii  ‘.m’  format  which  can  be 
run  on  any  platform,  routines  related  to  the  sparse  matrix  operations  were  written  in  FORTRAN  and 
compiled  as  ‘.mex4’  code.  These  routines  must  be  re-compiled  when  changing  to  a  different  platform. 
The  routines  are  bspars.mex4(bspars.f),  fillmat.mex4(fillmat.f),  and  spqrmex.mex4(available  from 
MATLAB® technical  support). 
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APE  Control 

Setup 


LOAD 

ERA 

Display 

AnIm. 

Analysis 

Opt.  sens, 

Man 

Vis 

Sens.  Anal.  Tol. 

Solve 

APE 

Close  Els 


Figure  B.2  APE  control  menu. 

The  program  is  initiated  from  the  MATLAB® command  line  with  the  command  ‘apeware’. 
Program  control  is  facilitated  through  the  main  control  panel  as  depicted  in  Figure  B.2.  Model  files  are 
loaded  into  the  workspace  hy  depressing  the  [LOAD]  button.  The  model  files  consist  of  information 
on  the  finite  element  model  matrices,  grid  data,  and  element  connectivity  data.  A  complete  description 
of  the  required  data  in  the  model  file  along  with  a  list  of  utilities  which  can  he  used  to  generate  the  data 
is  included  in  Section  B.5.  The  remaining  buttons  are  described  below. 

[ERA.]  Initiates  the  ERA  submenu  used  for  system  identification  as  further  defined  in  section  B.  1 . 

[Display.]  Produces  a  wire  frame  drawing  of  the  model,  showing  all  instrumented  degrees-of- 
freedom.  A  typical  display  is  shown  in  Figure  B.3. 

[Anim.]  Initiates  the  animation  submenu  used  to  display  and  animate  deformed  geometry  for 
any  mode  contained  in  the  model  file.  Modes  may  be  produced  either  from  the  measured  data  using 
ERA,  or  from  the  results  of  an  eigenanalysis  on  the  FEM  model.  This  feature  is  further  described  in 
Section  B.2. 
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Model  View 


Figure  B.3  Geometry  display  screen. 

[Opt.  sens.]  Performs  an  analysis  to  prioritize  the  degrees-of-freedom  to  instrument  based  on 
the  modes  selected  (using  the  [Man]  selection)  and  the  tolerances  set  in  [Tol].  This  feature  is  further 
discussed  in  Section  B.3. 

[Man.]  Opens  a  dialogue  box  as  shown  in  Figure  B.4,  allowing  the  user  to  input  the  measured 
modes,  the  instrumented  degrees-of-freedom  and  the  element  numbers  of  the  elements  used  for  the  APE 
method.  The  instrumented  degrees-of-freedom  and  the  element  numbers  are  updated  automatically 
when  using  the  [Vis]  and  [Sens.  Anal.]  selections  respectively,  as  discussed  below. 


Modes: 
Sensor  DOFs: 
Elements: 


Case  Set 

OK 


1:5 

1  4  7  10  28  3T  ^ 
[32,  47,  52:60,  78] 


Figure  B.4  Manual  input  dialogue  box. 
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Figure  B.5  Sample  graphical  sensor  selection  window. 


[Vis.]  Produces  a  wire-frame  drawing  of  the  model,  and  includes  additional  buttons  which  allow 
the  user  to  select,  via  the  pointing  device,  the  degrees-of-freedom  to  instrument.  The  [Direction] 
button  shows  the  current  degree-of-freedom  direction  (X,  Y,  Z,  Rx,  Ry,  Rz)  corresponding  only  to 
the  degrees-of-freedom  included  in  the  model.  Repeatedly  depressing  this  button  toggles  through  the 
selections.  Selecting  the  [Input]  button  allows  the  user  to  select  degrees-of-freedom  by  clicking  (right 
mouse  button)  on  the  node  to  instrument.  Nodes  which  are  selected  can  be  deselected  in  a  similar 
fashion.  All  selected  nodes  for  the  current  degree-of-freedom  direction  are  highlighted.  The  selection 
process  is  terminated  by  making  the  final  selection  with  the  left  mouse  button.  The  [View]  button 
initiates  the  3-D  viewer,  used  to  rotate  the  wire  frame  for  ease  of  node  selection.  The  [DONE]  button 
terminates  the  [Vis]  command.  A  typical  screen  display  is  depicted  in  Figure  B.5.  The  3-D  viewer  is 
shown  in  Figure  B.6. 

[Sens.  Anal.]  Initiates  the  sensitivity  analysis  for  damage  detectability  and  damage  localization. 
This  feature  is  described  in  Section  B.3. 
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Figure  B.6  Sample  3-D  view  selection  control. 


[Tol.]  Initiates  the  tolerance  menu  containing  the  tolerance  values  for  sensor  prioritization, 
damage  detectability,  damage  localization,  and  the  APE  method.  The  tolerance  menu  is  depicted  in 
Figure  B.7. 

[APE.]  Performs  the  APE  method  for  the  search  space  defined  in  [Man]  using  the  current 
tolerance  values.  Graphical  results  obtained  using  APE  are  described  in  section  B.4. 

[Els.]  Opens  a  dialogue  box  which  allows  the  user  to  key  in  element  numbers.  Depressing  the 
[SHOW]  button  displays  the  keyed  in  elements  on  the  wire-frame  model  in  red. 

[Close.]  Terminates  the  software  package. 

B.l  System  Identification 

Identification  of  modal  parameters  from  measured  frequency  response  functions  is  fully  auto¬ 
mated  using  the  apeware  package.  The  process  is  initiated  by  depressing  the  [ERA]  button  from  the 
main  control  panel,  which  displays  the  ERA  control  panel  as  shown  in  Figure  B.8.  If  the  measured 
transfer  functions  are  not  included  in  the  model  file,  fhey  can  be  loaded  info  fhe  workspace  eifher  from 
fhe  command  line  or  fhrough  fhe  [LOAD]  bullon  on  fhe  main  confrol  panel.  Informalion  on  how  fhe 
dafa  was  recorded,  (sample  rale,  number  of  sample  poinls,  elc.)  should  be  keyed  in  using  fhe  [ERA 
specs]  bullon,  which  iniliales  fhe  parameters  confrol  panel  as  shown  in  Figure  B.9.  An  explanalion  on 
each  parameter  is  oblained  by  depressing  fhe  [HELP]  bullon.  After  selling  Ihe  ERA  parameters,  Ihe 
ERA  melhod,  including  Ihe  inverse  Fourier  Iransforms  lo  oblain  Ihe  lime  domain  dala,  is  performed 
by  depressing  Ihe  [ERA]  bullon.  The  slate-space  idenlified  quadruple  is  conlained  in  Ihe  workspace  in 
Ihe  variables  (aera,bera,cera,dera),  in  eilher  conlinuous  or  discrete  form  as  determined  by  Ihe  selling 
in  Ihe  ERA  parameters  panel.  If  desired,  a  display  of  Ihe  singular  values  of  Ihe  block  Hankel  malrix  is 
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Figure  B.7  Tolerance  menu. 


provided  to  allow  the  user  to  select  the  number  of  singular  values  to  retain,  using  the  pointing  device. 
Simulated  frequency  response  functions  of  the  identified  state-space  quadruple  can  he  computed  and 
displayed  according  to  the  value  entered  in  the  [plot  fit]  parameter  hox.  A  comparison  of  a  simulated 
frequency  response  function(s)  of  the  identified  model  can  fhen  he  ohfained  and  displayed,  as  shown  in 
Figure  B.  1 1 .  An  addifional  opfion  fo  include  a  low -pass  filler  on  Ihe  measured  dala  is  also  incorporated. 

The  exlraclion  of  Ihe  modal  parameters  (frequencies  and  shapes)  from  Ihe  measured  dala  is 
performed  using  Ihe  [Extract]  hutton.  The  resulting  control  screen  is  shown  in  Figure  B.12.  Modes 
are  selected  hy  using  the  pointing  device  to  click  on  a  modal  peak  in  the  displayed  frequency  response 
function(s).  The  eigenvalue  of  the  model  corresponding  to  the  closest  mouse  pick  is  automatically 
determined,  along  with  the  partial  eigenvector  for  this  mode.  The  selected  partial  mode  shape  is 
displayed  along  with  a  polar  plot  indicating  the  phase  of  the  identified  eigenvector.  When  used  on 
measured  data  of  structures  with  minimal  structural  damping,  the  phase  of  the  eigenvectors  should 
lie  along  the  abscissa  (corresponding  to  real  eigenvectors).  The  polar  plot  is  used  as  a  check  on  the 
identified  eigenvector.  If  the  phase  is  scattered  off  the  abscissa,  the  mode  was  poorly  identified.  Each 
identified  eigenvalue  and  eigenvector  must  be  assigned  a  mode  number.  The  software  automatically 
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compares  the  identified  partial  eigenvector  with  the  stored  eigenvectors  of  the  analytical  model.  A 
figure-of-merit  (FOM)  is  listed  on  the  display  screen  along  with  the  mode  number  of  the  analytical 
model  with  the  highest  FOM.  The  FOM  is  a  measure  of  the  colinearity,  with  a  maximum  value  for 
perfect  colinearity  at  unity.  Depressing  the  [Store]  button  saves  the  identified  eigenvalue/eigenvecfor 
pair  for  fufure  use.  The  mode  numbers  can  be  changed  by  changing  fhe  enfry  in  fhe  displayed  mode 
number  box.  This  pairing  of  fhe  idenlified  and  analyfical  modes  is  performed  fo  accounf  for  mode 
swaps  occuring  in  fhe  damaged  sfrucfure.  Note  fhaf  when  using  only  parfial  eigenvectors,  no  guarantee 
is  made  fhaf  fhe  highesf  FOM  corresponds  to  fhe  correcf  mode  number,  and  fhus  care  musf  be  exercised 
when  assigning  mode  numbers.  The  [Clear  inodes]  buflon  can  be  used  to  clear  fhe  previously  sfored 
modes  from  fhe  workspace.  The  currenfly  stored  modes  are  listed  on  fhe  display  screen,  and  is  updafed 
after  each  store  operafion.  The  mode  selecfion  process  is  repeated  by  depressing  fhe  [Select  again] 
button  for  each  mode  desired.  The  [Zoom]  button  can  be  used  to  help  the  selection  of  closely  spaced 
frequency  peaks.  A  green  circle  is  displayed  on  each  peak  selected  to  aid  in  the  identification  process. 
Modal  peaks  can  be  selected  from  any  frequency  response  function.  Because  the  frequency  functions 
are  all  based  on  the  same  ERA  identified  model,  they  all  contain  the  same  modal  information.  However, 
some  modes  may  not  be  visible  from  the  selected  transfer  function.  Different  transfer  functions  are 
selected  by  changing  the  [Input]  and  [Output]  numbers  on  the  ERA  control  menu.  Multiple  transfer 
functions  can  be  displayed  by  including  several  combinations  in  the  Input/Output  boxes.  The  [clear] 
button  is  used  to  clear  the  plotting  window. 

B.2  Model  Tuning 

Model  tuning  is  performed  using  ASTROS-ID,  a  software  package  independent  of  the  apeware 
package.  ASTROS-ID  however  lacks  a  graphical  pre  and  post-processor.  Apeware  partially  accomo¬ 
dates  this  deficiency  by  providing  an  easy  to  use  graphical  interface  to  display  deformed  geometries  and 
animate  mode  shapes.  This  information  is  useful  in  classifying  mode  shapes  and  correlating  analytical 
and  measured  mode  shapes.  This  feature  is  initiated  using  the  [Auim]  button  (animation)  on  the  main 
control  panel.  The  mode  view  control  panel  is  depicted  in  Eigure  B.13.  Several  types  of  displayed 
results  are  possible.  Deformed  geometries  using  the  wire-frame  model  can  be  displayed  for  any  stored 
mode,  from  either  the  analytical  data  as  in  Eigure  B.14,  measured  data,  or  both.  Additionally,  mode 
shapes  can  be  animated  at  fifteen  frames  per  second  to  aid  in  the  classification  process.  Alternatively, 
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ERA  Control 

INPUT#; 

OUTPUT  #: 

P__  Measured  | 

J  ERA  FIT  J 
ERA  Specs  | 

ERA  I 

Plot  FRF  I 

Extract  Shapes  | 

V-  Zoom 


DONE 


Figure  B.8  ERA  control  menu. 


clear 
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zT 


APE  Info  Window 


ERA  Parameters 


OUTPUTS  =  #  of  sensors 
INPUTS  =  #  of  actuators 

LH  FACTOR  is  a  factor  2  <  Ih  <  length(FRF)/ns  to 
determine  the  size  of  hankel  matrix  typicaiiy 
3-8.  The  resuiting  hankel  matrix  is  nxn  where  n 
is  approximately  ns*nmeas*ninput*lh/2 
STATES  =  #  number  of  states  for  model 
INTEGRATIONS  =  #  of  integrations  on  data 
(set  to  1  for  accei  to  veiocity,  etc.) 

FiLTER  ORDER  =  butterworth  fiiter  order,  ignored 
if  break  is  negative 

BREAK  =  cutoff  for  low  pass  filter  of  impulse 
response  in  Hz.  set  break  <  0  for  no  filtering 
ZEROS  =  #  of  zeros  at  low  freq  end  of  spectrum 
PLOT  =  figure  number  for  first  plot  (<0  to  surpress 
0  or  D  =  1  for  continuous  model  or  2  for  discrete 
CUT  =  1  to  prompt  for  singular  value  cutoff 
NDPTS  =  number  of  original  sample  points 
(ie  2"'N,  512,1024,2048)  set  to  zero 


Figure  B.IO  ERA  help  menu  dialogue  box. 


Model  View 


Figure  B.ll  Sample  frequency  response  showing  measured  data  and  simulated  data  from  an  ERA 
identified  model. 
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APE 

Mode  #: 

8 

FOM:  0.9741 

Select  again 

Save  Vector  | 

Clear  modes 

Done 

stored  modes:  [1:8] 


35.886Z  Hz  mode 


Figure  B.12  Control  panel  to  extract  measured  modes. 


the  mode  shapes  can  he  viewed  using  a  line  plot  format,  using  either  the  full  length  eigenvectors  as 
in  Figure  B.16,  or  the  partial  eigenvectors  corresponding  to  the  instrumented  degrees-of-freedom  as  in 
Figure  B.  15.  Vector  normalization  and  sign  convention  are  automatically  accounted  for  in  the  displayed 
results. 

Interactive  compiled  FORTRAN  routines  (db2mate,  db2xyz,  db2dvs)  can  be  used  to  extract 
the  eigenvalues  and  eigenvectors,  the  finite  element  model,  and  the  tuned  design  variables  from 
the  ASTROS  data  base.  The  results  are  stored  in  MATLAB®binary  format  for  use  with  apeware. 
Additionally,  the  MATLAB®routines  wrtshape  and  wrtfreq  can  be  used  to  generate  the  TSHAPE  and 
TFREQ  data  cards  from  the  identified  dafa,  for  use  in  ASTROS-ID. 


B.3  Sensitivity  Analysis 

The  sensifivify  analysis  is  performed  using  fhe  stored  analyfical  model.  The  analysis  is  performed 
on  fhe  modes  keyed  in  using  fhe  [Man]  buffon  on  fhe  main  confrol  panel.  Parameters  for  bofh  fhe 
damage  localizafion  and  fhe  sensor  priorifizafion  are  defermined  from  fhe  keyed  in  enfries  on  fhe 
tolerance  menu  obfained  using  fhe  [Tol]  buffon.  The  sensor  priorifizafion  is  performed  by  depressing 
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Figure  B.14  Sample  display  of  deformed  geometry. 
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the  [Opt.  Sens.]  button.  The  results  of  the  sensor  prioritization  are  shown  in  two  separate  display 
windows.  The  first  is  the  gradient  information,  depicting  the  values  of  for  each  sensor  location  as 
shown  in  Figure  B.  17.  Also  listed  are  the  degrees-of-freedom  with  values  of  below  the  detectability 
threshold  level,  and  degrees-of-freedom  within  the  colinearity  threshold  level.  A  second  figure  lisfs 
fhe  priorifized  locafions  in  order,  sfarfing  from  fhe  upper  leff,  as  shown  in  Figure  B.18. 

Damage  localizafion  is  performed  by  depressing  fhe  [Sens.  Anal.]  buffon  on  fhe  main  con- 
frol  panel.  The  resulfs  are  displayed  in  fwo  separate  figures.  Eigenvalue  and  eigenvector  gradienf 
informalion  confained  in  Dx  and  is  displayed  as  shown  in  Figure  B.19.  The  colinearify  analysis 
resulfs  are  displayed  as  shown  in  Figure  B.20,  which  lisfs  fhe  symmefric  elemenfs  as  well  as  fhe  un- 
defecfable  elemenfs.  Elemenf  numbers  are  color  coded,  wifh  each  symmefric  sef  assigned  to  a  disfincf 
color.  Eor  reference,  a  wire-frame  model  is  also  displayed  wifh  fhe  elemenfs  color  coded  according 
fo  fhe  displayed  informalion.  Brighl  green  elemenfs  are  classified  as  undeleclable  and  white  elemenfs 
correspond  lo  elemenfs  which  can  be  uniquely  idenlified. 
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APE 


APE  rank  ordered  selected  sensor  locations: 

12  28  22  16  18  36  d  30 

8  1  3  35  33  9  11  27 

25  13  15  23  21  17  19 


Figure  B.18  Sample  sensor  location  prioritized  list. 


Figure  B.19  Sample  damage  localization  gradient  results. 
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APE  undetectable  elements;  APE  symmetric  elements: 

1  2  5  6  11  16 
21  26  31  36  37  dO 

18  19 
41 

Modes  used  in  APE:  4  5  6  7  8 

Figure  B.20  Sample  damage  localization  analysis  results. 

B.4  Damage  Identification 

Damage  identification  using  APE  is  performed  using  the  [APE]  button  on  the  main  control 
panel.  The  APE  iteration  proceeds  according  to  the  parameters  set  in  the  tolerance  menu  as  shown 
in  Eigure  B.7.  The  iteration  tolerance  values  are  explained  using  the  [HELP]  button.  The  numerical 
results  of  the  APE  method  are  displayed  as  shown  in  Eigure  B.21.  If  the  [show  graphics]  box  is  set  to 
1,  the  results  are  also  shown  in  graphical  format  as  depicted  in  Eigure  B.22,  with  the  identified  damaged 
elemenfs  color  coded  according  fo  fhe  percenf  damage  idenlified.  If  [show  graphics]  is  sef  fo  2,  fhe 
intermediate  resulfs  for  each  iferafion  step  are  displayed. 

B.5  Model  File  Requirements 

The  following  is  a  lisf  of  variables  and  fheir  descripfions,  which  musf  be  confained  in  a 
MATEAB® binary  file.  This  file  musf  be  loaded  info  fhe  workspace  after  inifiafing  fhe  apeware  pro¬ 
gram.  The  exisfing  workspace  is  purged  when  fhe  program  is  inifiafed.  Note  also  fhaf  MATEAB®is 
case  sensifive. 
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Ape  Results 


r ' 


Figure  B.22  Sample  damage  identification  graphical  results. 
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1 .  M  Finite  element  mass  matrix. 

2.  K  Finite  element  stiffness  matrix. 

3.  C  Identity  matrix  the  size  of  M. 

4.  eol  Sorted  eigenvalues  of  M,  K. 

5.  vol  Sorted  eigenvectors  of  M,  K. 

6.  B  Matrix  as  defined  in  Equation  (6.10)  or  (6.13)  such  that  Ki  =  BiBf . 

7.  nb  Number  of  columns  in  B  corresponding  to  a  single  element  (scalar). 

8.  DOFTYPE  Number  of  degrees  of  freedom  per  node  (scalar). 

9.  xyz  Matrix  containing  the  grid  point  coordinates,  one  (x,y,z)  triplet  for  each  node.  Matrix 
size  is  (#  of  nodes  x  3). 

10.  con  Matrix  whose  row  points  to  the  row  of  grid  points  in  xyz  for  the  element. 

Matrix  size  is  (#  of  elements  x  2). 

1 1 .  CONINDX  Vector  which  points  to  rows  in  con  that  have  columns  in  B  (#  of  design 
variables). 

12.  connodes  Vector  listing  nodes  which  are  constrained. 

13.  DOFSTR  String  vector  of  degree-of-freedom  labels,  (‘X’,  ‘Y’,  ‘Z’,  ‘Rx’,  ‘Ry’,  ‘Rz’), 
used  in  the  model. 

B.5.1  Utilities.  To  aid  in  construction  of  some  of  the  model  file  requiremenfs,  fhe  following 
ufilifies  are  available. 

1 .  db2niate  Execufable  roufine  fo  exfracf  fhe  mass,  sliffness,  circular  frequencies,  and  eigen¬ 
vectors  from  an  ASTROS  dafa  base. 

2.  db2xyz  Execufable  roufine  fo  exfracf  fhe  xyz  and  con  mafrices  from  an  ASTROS  dafa 
base. 

3.  db2dvs  Execufable  roufine  fo  exfracf  fhe  funed  design  variables  and  elemenfal  properfies 
from  an  ASTROS  dafa  base. 
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4.  bildbeam  MATLAB®.m  file  to  construct  the  B  matrix  for  beam  elements  from  the 
elemental  properties. 

5.  bildrod  MATLAB®.m  file  to  construct  the  B  matrix  for  rod  elements  from  the  elemental 
properties. 

6.  eigsrt  MATLAB®.m  file  to  construct  the  sorted  eigenvalue  (eol)  and  eigenvectors  (vol) 
matrices  from  M  and  K . 
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